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SECTION  I 

MATERIALS  PROPERTIES 


1.  DENSITY  OF  COMPOSITES 


Mass  Fractions 

Based  on  conservation  of  mass,  the  rule- of - 


mixtures  mass- 

M. 

. + M,+  • • • M 

2>, 

M 


EMi  = M O) 

E™.  = 1 (2) 


Volume  Fractions 

Based  on  existence  of  definable  constituent 
volumes,  the  rule- of-mixtures  volume- 
fraction  relation  is  assumed  to  be  valid: 


c.  Mass  and  Volume  Fractions  for  Density 


lM! 

m2 

. . , 

:M«i 

mi 

m2 

* * * 

JjlJ 
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Figure  1 Mass  fractions  must  al- 
ways add  up  in  a compos- 
ite with  n phases . 
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Figure  2 Volume  fractions  need 
not  always  add  up  for 
many  reasons;  phase 
boundary  not  defined, 
phase  interactions  due  to 
mixing  and  binding,  voids, 
etc. 
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All  equations  are  equally  applicable  and  are  limited 
by  the  basic  assumption  that 


M. 


pi  Vi 


holds  within  each  constituent  phase. 

Equation  7 can  be  rewritten  as  follows: 

m. 


p . 


(8) 


(9) 


Figure  3 Rule- of-mixtures  relations 

imply  non-interacting  phases, 
within  each  of  which  Equa- 
tion 8 holds . 


d.  Density  of  Two- phase  Composites 


From  Equation  5 


P = 


<Vl+P2V2 


Fiber  Fraction 
v2  °r  m2 


From  Equation  7 

1 _ ^1  + ^2 

P “ P,  P, 


v2  or  m2 


TABLE  1 TYPICAL  CONSTITUENT 
DENSITIES 


MATERIAL 

SPECIFIC 

DENSITY 

GRAVITY 

kg  m3 

Kevlar 

1.45 

1450 

Graphite 

1.7 

1700 

Glass 

2.  6 

2600 

Boron 

2.6 

2600 

Steel 

7.8 

7800 

W 

19.  3 

19300 

Nylon 

1.1 

1100 

Epoxy 

1.2 

1200 

Polyester 

1.4 

1400 

Be 

1.8 

2.8 

1800 

A1 

2800 

Ti 

4.5 

4500 

Figure  4 Linear  plots  are  the  easiest  provided  proper 
variables,  p or  1/p  is  taken.  Relationships 


between  mass  and  volume 

seen. 

CO 

nte 

nts  are  read 

lily 

I 

I- — "t- i — ’ y“  t — — •* f-  - -* t- ■ - 

' 
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Void  Content 


Void  content  in  a two-phase  composite  can  be  defined  by- 


calculated 


This  relationship  is  not  accurate  because 


(11  The  assumed  non-interacting  phases  ignore  curing  stress  which  can  induce  up  to 
1 percent  strain. 


(2)  Absorbed  moisture  which  can  induce  swelling  of  several  percent 


(3)  Actual  voids  and  cracks  may  be  closed,  thus,  cannot  influence  the  gross 
density  beyond  the  detectable  level. 


Alternative  methods  for  determination  of  voids  will  be  covered  later 


2.  ONE-DIMENSIONAL  ELASTICITY  PROPERTY 


a.  ModuluB  vfl  Compliance 


E = Modulus  = — 


S = Compliance  = 


Stress 

Strain 

Strain 

Stress 


(12) 


(13) 


E • S = 1 


(14) 


Figure  6 Modulus  and  compliance  are  re- 
ciprocal of  each  other.  Either 
one  can  be  used  but  one  is  usually 
preferred  for  a given  situation. 
This  will  be  illustrated  on  sever- 
al occasions  later. 
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b.  Elasticity 

Property 

of  2 - Phase  Composites 

Figure  7 

For 

a c 

:onstant  strain  model 

with 

Figure 

8 

For  a 

constant 

stres 

s n 

nodel  with 

parallel  phases,  modulus  E is  the 
preferred  property  because  simple 

rule  of  mixture  eauation  annlies- 

in- 

is 

■series  phases,  compliance  S 
the  preferred  property  because: 
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TABLE  2 TYPICAL  CONSTITUENT  STIFFNESSES 


MATERIAL 

YOUNG'S  MODULUS 

E(GPa) 

ACOUSTIC 

VELOCITY 

C(km  sec 

Kevlar 

131 

9.50 

Graphite 

207 

11.0 

Glass 

87 

5.78 

Boron 

414 

12.6 

Steel 

207 

5.15 

W 

407 

4. 59 

Nylon 

4.8 

2.  1 

Epoxy 

3.4 

1.7 

Polyester 

3.4 

1.6 

Be 

241 

11.6 

A1 

69 

5.0 

Ti 


103 


4.8 


3.  TWO-DIMENSIONAL  ELASTICITY  PROPERTIES 


a.  Strains  in  Terms  of  Stresses  for  Isotropic  Bodies  Under  Plane  Stress 


(1)  Compliance  matrix  in  conventional  form 
1 v 

6x  = E °x'  E °y 


v 1 

TT  <T  + —a 
E x E y 


(3)  Compliance  matrix  in  matrix  form 


xy 
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cr 

xy 

^ xyy 

L gJ 

(19) 


(2)  Compliance  matrix  in  tabular  form 
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(4)  Compliance  matrix  in  index  form 
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S11  S12  0 


S21  S22 


(5)  Contracted  Notation: 
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t>.  Stresses  in  Terms  of  Strains  for  Isotropic  Bodies  Under  Plane  3tre 


(I)  Modulus  matrix  in  conventional  form 


E 

vE 

f > 

"x=  1- 

2 e 
V 

X 2 y 

1-  V 
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vE 
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y i. 
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+ , e (22)  < 

* , 2 y ' 
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J € 
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cr 

n.  xy^ 

(2)  Modulus  matrix  in  tabular  form 
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c.  Components  of  S. . and  Q... 


(3)  Modulus  matrix  in  matrix  form 

• 1 P4  ol  r. 


E 

vE 

l-v2 

l-v2 

vE 

E 

l-v2 

l-v2 

0 

0 

(4)  Modulus  matrix  in  index  form 

"ll  fQll  °12  0 1 


°21  °22 


Summation  convention 


°66j  le 6 


EQ..  e . 
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i.  j = 1,  2,  6 
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Bodies 
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e.  Kronecker  Delta 

Analogous  to  Equation  14, 
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SECTION  II 
STRESS  AND  STRAIN 
1.  STRESS  AS  COORDI  NATES  ROTATE 


a.  Introduction 

Stress  is  NOT  defined  as  (You've  flunked) 

It  is  defined  by  how  it  changes  as  its  reference  coordinates  change.  That's  why  we  need  to 
know  the  ground  rules.  Rigid  rules  are  particularly  important  for  composites  because 
composite  properties  also  change  with  coordinates,  by  a different  set  of  rules.  All  of  these 
rules  are  called  the  transformation  relations.  In  this  section,  we  will  discuss  only  those 
related  to  stress  and  strain. 


When  we  apply  a stress  to  an  isotropic  material,  we  can  analyze  the  response  of  the  material 
in  the  same  coordinate  system.  There  is  no  reason  to  look  into  any  other  coordinates  with 
the  possible  exception  of  the  plane  where  shear  stress  is  maximum. 

For  composites,  the  response  is  highly  dependent  on  the  orientation  of  the  material.  It  is 
therefore  important  to  know  how  an  applied  stress  can  be  transformed  to  the  material  axes. 

( See  Figure  1 1 ) 


' r 

i r fVi  m ■ f *•  « • a * ■ 

r U1C  oTTuSSCS  III 

A% 

y,  we  warn  to 

* 

know  what  apply  stresses  must  be  for 
any  coordinate  system,  we  simply 
rotate  the  0 backward,  or  apply  inverse 
transformation.  The  same  rule  applies 
except  now  we  use  negative  0 in  the 
formula. 

It  is  therefore  an  important  rule  of 
transformation  to  know  the  positive 


Figure  11  Transformation  of  applied  stress 
in  coordinate  x-y  to  material 
coordinate  x'-y'  by  a positive  rota- 
tion of  0. 


from  the  negative.  We  are  all  aware  of 
the  difference  between  tension  and  com- 
pression, but  we  rarely  pay  attention  to  positive  from  negative  shears  because  they  are  not 
important  for  isotropic  materials.  For  composites,  we  must  make  sure  the  proper  signs 
for  all  angles  of  rotation,  stresses,  etc.,  are  used.  For  rotation,  we  will  use  the  right 
handed  system. 
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b.  Formula  for  Stress  Transformation 

TABLE  3 STRESS  TRANSFORMATION  RELATIONS 


I 

R 

cr' 

X 

1 

cos  2(9* 

*6) 

tr' 

y 

1 

■ cos  2(9*6) 

cr' 

xy 

0 

-sin  2(9^6) 

.X" 


I = 


<r  + a 
1 


xy 


2 or 


tan  2 6|  = 


a - <r 
x y 


(28) 

(29) 

(30) 


Figure  12  New  coordinates  x'-  y'  in 
terms  of  old  coordinates 
x-y,  or  given  x-y  and  6,  find 
x'-y1.  Arrow  of  rotation  is 
pointing  up.  If  0 is  negative, 
new  coordinates  are  x"-y". 


6 = 6.  + 90  if  <r  < a ; 6 = 6.  if  a > a 
1 x y 1 x y 

Both  I and  R are  invariants,  as  shown  in  the  Mohr's  Circle. 
Special  orientations; 

(1)  When  9 = 6,  i.e..  Principal  Directions 


06 


<r'  , a' 
x y 


*y 


= maximum  or  minimum 
= 0 


(2)  When  9 - 6=±  »/ 4,  i.e.,  max.  shear 
orientation 


<7f  = 

X 

_ i _ 

0*' 

y 

= I 

u 

*y 

max , = 

R 

1 r ! — f — 

TABLE  4 

' ] 1 1 I T r 

ALTERNATIVE  FORMULA  FOR 

STRESS  TRANSFORMATION 
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nc 

Figure  13  Mohr’s  Circle  is  defined  by 
invariants  I and  R.  Phase 
angle  6 is  defined  by  a spe- 
cific combination  of  stress 
components  ax,  <s  and  sxy. 
As  reference  coordinates 
change  by  49,  the  rotation  in 
Mohr's  Circle  is  29. 

Special  attention  should  be 
given  to  value  of  6; 


6 = 


6,  if  o-  > tr 
1 x y 

6,  4-90  if  cr  < cr 
* 1 x y 
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c.  Graphical  Illustrations  of  Trigonometric  Relations 


Figure  14  Relevant  trigonometric  functions  for  stress  transformations.  Phase  angle  6 
displaces  the  function  to  the  right  (vertical  axis  to  the  left). 


4 
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2.  MAXIMUM  STRESS  FAILURE  CRITERION  FOR  UNI  01 RECTIONAL  COMPOSITES 


Off-Axis  Tensile  Failure  Stresses  Under 
Uniaxial  Stress  <r  x(cry  =0^=0).  The  stresses 
in  the  material  symmetry  axis  x'  (along 
the  fibers)  can  be  obtained  directly  from 
Table  4. 


Longitudinal  stress  = m <r 


Figure  16  Coordinates  for  off-axis  uni- 
axial tests.  Shear  stress  will 
be  positive  if  fibers  oriented 


mn<r 


Maximum  stress  criterion  assumes  that 
failure  will  occur  when  lowest  of  the 
following  3 <r  's  is  reached. 


' Jew.  f/tfLOgg 

fe*Jf  fAUMttr 

,5h€A*  ft/LMe 

/ <S*YJ 


mn 


e 17  Maximum  stress  criterion  is 
typified  by  separate  branches 


which  may  be  related  to  fiber 


and  matrix  (transverse  or 


shear)  failures 
interaction. 


There  is  no 


ssive  Uniaxial  Strength 


/ / SfKM 

<\  ' THAW.  Cerrfft 

• V 

' \ Lo*4<r.  4f/V. 


ear  strength  S is  assumed  to  be 
positive  and  negative  shear  stre 


Maximum  stress  criterion  as 
applied  to  off-axiB  compres- 
sive tests.  There  is  no  inter 
action  between  tensile  and 


Where 


compressive  strength 


Fi 

. (- 4 

jure  18 



— 

it 



Where 

11 

n a 

s cos  0 

r l 

n 

3 

= a 

iin  0 

1 ( 

c.  Off-Axis  Shear  Strength 

T 

y 


Substituting  <rx=  <rxr  = 0 into  Table 


4 (for  + 6), 


cr  = 2mn  o = ( sin  20)o 

x xy  xy 

o = -2mn  o = -(sin  20)o  > (341 

y xy  xy  ' 

o = (m^-n^)o  = (cos  20)o 

xy  xy  xy- 


Depending  on  the  signs  of  the  stress 
components  o^,  Oy,  o^y,  which  are 
functions  of  the  signs  of  0 and  oXy, 
the  following  strength  criteria  shall 
be  used.  Note:  Each  quadrant  in  the 
table  has  distinct  combinations. 


- 0 

T 0 

^xy 

X'/sin  20 

Y/sin  20 

X/sin  20 

Y'/sin  20 

- cr 

xy 

X / sin  20 

Y1  / sin  20 

X'  / sin  20 

Y/sin  20 

5 i 

'v 

\ / , 

i x 

' > ‘ /v' 

* / t / 2 

w V ' 

rs- 

7/&*{  Ten s. 

l^y,  't* 

\\  , 

' V . , '' 

\ \ > 

1 v i - ' > / Co*fM 

\ ! 

h' — '/\sFbi 

T/+  Till  6 

Y 

- . MW- 

7/?4+is 

W 

/ /t 

1 v ShCA*  i V 

X Xx,t 

/ V , \ ' 

( \ ' 

V ' ^ i 

V 

/ > } \ 

/>  / 1 « 

>,  A 

/ i » 

- 0-^-3 

4 ~ 4 

3 1 

4 4 

+ cr 

S/cos  20 

S'  /cos  20 

- cr 

xy 

S'/cos  20 

S/cos  20 

Figure  19  Off-axis  shear  strength  based  on 
maximum  stress  criterion.  Six 
possible  failure  modes,  including 
positive  and  negative  shears  are 
labeled.  In  this  case,  longitudin- 
al tensile  and  compressive  are 
never  limiting  cases.  Matrix 
tensile,  compressive  and  shear 
are  the  controlling  modes. 


Numerous  failure  criteria  other  than  the 
maximum  stress  will  be  discussed  later. 


d.  The  intersection  of  any  2-failure  curves  in 
Figure  17  are  as  follows: 

,-l  X 


1 S' 

0 - ,,  -1  Y 

(36) 

°2  = tan  S' 

0^  = cot 


■JT 


(37) 
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3.  STRENGTH  ESTIMATE  OF  RANDOM  FIBER  COMPOSITES 


Tensile  Strength 

It  may  be  assumed  that  the  strength  of  a random  fiber  composite  is  equal  to  the  average 
strength  of  the  area  under  the  off-axis  tensile  strength  of  the  area  under  the  off-axis 

— I-  „ r _ ' J • . • * . 


2 

rn 

= * \ 

f 

) 

xe  de 

2 

N 

fei  de 

i ~ /82-  de 

r12  de 

ji 

3 oos^e 

J sin  6 cos  e 
el 

“ + Y 

( 

*2  sin^e 

After  integration 


X 
Y 

Where  a 


iW-  -w 


s 

Y 


(38) 

(39) 


(40) 


(41) 


X/Y 

Figure  20  The  analytic  estimate  (based  on  high  shear  strength)  falls  below  most  available  data 
by  a factor  of  1.5;  but  the  trend  appears  to  be  in  general  agreement.  The  rule-of- 
mixtures  equation,  with  correction  factor  p,  is  shown  in  dashed  lines,  [j] 
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b. 

Compre 

issive 

Strength 

c.  Shear  Strength 


d.  Strengths  of  Particulate  Composites 


4.  STRAI N AS  COORDI  NATES  ROTATE 


a.  Formula  for  Strain  Transformati 


Only  the  shear  component  is  the  difference  between  the  strain  and  stress  transformation 
relations,  i.e.. 


TABLE  7 STRAIN  TRANSFORMATION  RELATIONS 


Figure  21  Positive  and  negative  rotations  for 
transformation. 


s Circle  representation  of 


strain  transformation 


Only  the  one-half  factor 


for  the  shear  strain  Is  the 


difference  between  the  strain  representation  and  that 


for  stress.  Special  attention  must  be 

6 = 6 , if  e _>  e ; 6=6,  + 90.  f < 


i 

1 

X 

e' 

i 

y 

€ * 
xy 

0 

TABLE  8 TRANSFORMED  AND  PRINCIPAL  STRAINS 


6-45 

6-30 

1 

6-15 

6 

6 + 15. 

_6+30 

6+45 

0 

2.000 

3.346 

4.  332 

4.693 

4. 332 

3.  346 

2.000 

1 

2.000 

.654 

-.332 

-.693 

-.332 

.654 

2.000 

3 

5.  385 

4.664 

2.693 

0 

-2.693 

-4.664 

- 57  385 

5 

- 

. 

- — 

( 

i 

1 

I 

(2) 

Give! 

n 

« 

r.  = 

2 

1 

E 

3 

(mm/m)  or  ( 10  | 

pim/m) 

i 

1 

3 

J 

f 

V 

% 

7 

i = 

1. 

500 

1 

- 

-i 

.581 

6 

= 54.217 

I 
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Since  there  are  3 strain  components  at  each 
point,  3- element  rosettes  are  in  general 
needed  to  solve  for  3 unknowns.  Assuming 
3 elements  are  mounted  at  3 different  angles 
from  some  reference  coordinates,  each  rosette 
must  satisfy 


Figure  25  Strain  rosette  orienta- 
tion with  positive 
angles.  If  an  angle  is 
negative,  it  shall  be  so 
entered  into  Equations 
(45)  and  (46). 


Other  Rosettes 


Four- element  rosette  is  an  over- determined  system  when  a 4th  equation  is  added: 


Methods  of  solution  for  an  over-determined  (4  or  more  equations  for  3 unknowns)  are 


available.  The  additional  strain  gage  also  serves  as  a redundant  gage,  in  case  of  a 


defective 


Two- element  rosette  is  adequate  if  additional  information  on  the  strain  components  is 


given 


For  example,  if  from  symmetry  considerations,  « is  known  to  be  zero 
two-element  rosette  can  be  used.  The  angle  betwe’Si  the  elements  (l.e., 


value,  although  normally  it  ia  90° 


0,  e ^0;  or  only  uniaxial 


One- element  rosette  is  adequate  if  it  is  known  that  e 
strain  is  needed. 


STRAIN  RnSFTTFS 

I i i — rn — i — i — — i — i "t ■ i t — i — 

Q 

/ 

a.  Three 

-Element  Rosette 

s 

/•  < 

In  i 

matrix  form 
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c.  Sample  Problems. 


(1)  From  a uniaxial  test, 

A = 0,  B = 60°,  C = -60°,  and  e^,  e^.  = 10,  -3,  -3  mm/m. 


find  Poisson's  ratio  along  a-axis. 


Figure  26.  Rosette  for  uniaxial 
test. 

(2)  Find  off-axis  mounted  2-element  orthogonal  rosette  in  order  to  achieve  full  15  mm/m 
-ange,  Assume  ultimate  longitudinal  strain  of  3 percent  and  Poisson's  ratio  of  0.  3. 
Is  there  an  optimum  orientation? 

t 


A 
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SECTION  III 

STRESS  - STRAIN  RELATIONS 
1.  COMPLIANCE  AND  MODULUS  MATRICES 


a.  Stress- Strain  Relations  in  Longhand  Form  (Plane  Stress) 


S..CT.  , 
1J  J 

i,  j = 1,2,6 

(48) 

Slfj 

= Sllal+S12a2+Sl6a6 

S2jaj 

= S2lVS2lVS26a6 

(49) 

*6  = Vj  = S61 V S62a2 + S66a6 


<r.  = Q..e. 

1 !J  J 


(50) 


a = Q e.  = Q e +Q  e +Q.,e, 

1 111  111  122166 


aZ  = °2i*i  " Q216l + Q22*2+ Q26*6 


(51) 


°6i*i  = Q6l*l+  Q62®2+  °66*6 


b.  In  Matrix  Form 


€1 

S11  S12  S16 

V 

< 

*2 

= 

S21  S22  S26 

< 

a2 

> 

S/s  S-  S,, 

0-/ 

L.  Ol  KJC  DO 

^ vj 

( 52) 


S..  is  symmetric,  i.e.,  = S^.  or  S^  = S^,  SJ6  = S^,  . S^ 


I I I I 
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r n 
1 


= 


f N 

Q 

Q , 

Q . 

e , 

11 

12 

16 

1 

°21 

°22 

°26 

i 

e2 

► 

.°61 

°62 

°66. 

<e6> 

(53) 


Q..  is  also  symmetric,  i.e.,  Q. . = Q..,  or  Q = Q , Q,  , - Q,.,  Q-,,  = Q,, 
ij  ljji  1221166126  62 

c.  Elastic  Symmetries 

TABLE  10  COMPLIANCE  AND  MODULUS  IN  TERMS  OF  ELASTIC  SYMMETRIES 


Symmetry 
(No  Indep.  Const) 


Anisotropic 

(6) 

or 

Generally 

Orthotropic 

(4) 


Specially 

Orthotropic 

(4) 


Isotropic 

(2) 


Compliance  Matrix  S.. 

.1  13 

(TPa) 


11 


12 

E11 


16 


21 

E22 


61 


22 


62 


26 


12 


_1_ 

E, 


LT 


TL 

'E„ 


LT 


_1 

E 

_v 
‘ E 


v_ 
' E 

1 

E 


_i 

G 


Moduli  '"trix  Q. . 
(MPa) 


Q11  Q1 2 °16 


°21  °22  °26 


°61  °62  Q66 


mE 


vmE 


vmE 


mE 


mEL  VTLmEL  ° 


vLTmET  mET  0 


LT 
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For  orthotropic  material:  m 


For  isotropic  material 


can  be  expressed  in  terms  of  engineering  constants  only  for  orthotropic  and 


Note 


isotropic  materials,  not  for  anisotropic  materials.  S. . does  not  have  such 


limitations 


TABLE  11  ENGINEERING  CONSTANTS  IN  TERMS  OF  S, 


Q 

Q17 

12 

: 

n 

1 

°12 

°22 

Q. 

r% 

1C. 

Q. . 

* 

i 

1 

, 

Q, 

6 

6 

. 

, — j 

— t 

E 

J 

□ 

E 

L 

^ | 

1 

r ! 

1 | 

f 

(■ 

1 1 

1 

— f 

_ i 

d.  Elastic  Constants 


TABLE  12  ENGINEERING  CONSTANTS  FOR  UNIDIRECTIONAL  COMPOSITES 


, 

Type 

Material 

Fiber 

Vol. 

vf 

Specific 

Gravity 

y 

el 

GPa 

et 

GPa 

^LT 

glt 

GPa 

B<  41 

5505 

Mod  II 

5206 

HMS 

B/Ep 

Gr  /Ep 

Gr /Ep 

0.5 

0.55 

0.48 

2.0 

1.5 

1.58 

204 

55 

185 

18.5 

8.83 

6.76 

0.23 

0.  30 

0.20 

5.  79 

5.24 

5.  86 

3002M 

T300 

5208 

Mod  I 

ERLA  4289 

Mod  I 

ERLA  4617 

Gr/Ep 

Gr/Ep 

Gr/Ep 

0.70 

0.51 

0.45 

1.60 

1.56 

1.54 

181 

188 

190 

10.3 

4.  14 

7.  10 

0.28 

0.20 

0.  10 

7.  17 

4.  83 

6.2 

AS 

3501 

B(4) 

WRD  9 371 

Mod  I 

WRD  9371 

Gr/Ep 

B/PI 

Gr  /PI 

0 . 66 

0.49 

0.45 

1.60 

2.0 

1.54 

138 

222 

216 

8.96 

14.5 

4.97 

0.  30 

0.  16 

0.25 

7.  1 

7.7 

4.  5 

S Glass 
1009-26-5901 

Gl/Ep 

0.72 

2.13 

60.7 

24.8 

0.23 

12.0 

Scotchply 

1002 

Gl/Ep 

0.45 

1.8 

38.  6 

8.27 

0.26 

4.  14 

T300 

SP  313 

Gr/Ep 

0.65 

1.55 

140 

9. 7 

i 

U«  jt 

Kevlar  49 

Kev/Ep 

0.60 

1.38 

76 

5.5 

0.  34 

2.  3 

Epoxy 

TABLE  13  COMPLIANCE  S. . IN  TERMS  OF  ENGINEERING  CONSTANTS  ( 1012Paf  1 


Composites 

s . i 

S - 1 

O 1 

11  el 

22  ‘ Et 

S12  'VLTS11 

66  ‘ G,_ 
LT 



(TPaf‘ 

(TPa)'1 

= 'VTLS22 

(TPaf  1 

T - 300/5208 

5.525 

97.087 

-1.547 

139.476 

B / 5505 

4.902 

54.054 

-1.128 

172.712 

Scotchply  1002 

25.907 

120.919 

-6.744 

241 . 546 

Kevlar  49 

13. 158 

181.818 

-4.474 

434.783 

S Glass 
1009-26-5901 

16.474 

40. 323 

-3.789 

83. 333 

HMS  300277 

5.4054 

147.929 

-1.081 

170.649 

Mod  I 

WRD  9371 

4.6  30 

201.207 

-1.157 

222.222 

TABLE  14  MODULUS  IN  TERMS  OF  ENGINEERING  CONSTANTS  (106Pa) 


Composites 

m 

Qir  mEL 
(GPa) 

°22  = mET 
(GPa) 

°12=  vLTQ22 

r v Q 

Tl  1 1 

°66=  CLT 
(GPa) 

S..Q.W=6.. 

ij  jk  lk 

T-300/5208 

1.0045 

181.811 

10. 346 

2.897 

7.  170 

B/5505 

1.004S 

204.985 

18.589 

4.275 

5.790 

Scotchply 

1002 

1.0147 

39. 167 

8.  392 

2.182 

4.  140 

Kevlar  49 

1 .0084 

76.641 

5.546 

1.886 

2. 300 

1009-26-590 

1.0221 

62.041 

25.348 

5.830 

12.000 

HMS  3002M 

1.0011 

185.271 

6.770 

1. 354 

5.860 

1_vltet/el 
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e.  Determination  cf  Strains  from  Stresses  for  a 


<r.  in  MPa 


S..  for  Gr/Ep  in  (TPai 
ij 


.06  > TPa 


5.525  -1.547 

1.547  97.087 


1 39.470 


Figure  27  Coordinates  for  cr^  and 
Sjj  are  coincident. 


e.  in  mm/m 

l 


Figure  28.  For  a given  set  of  Imposed  Stresses,  shown  are  the  resulting  strains  for 
T-300/5208  unidirectional  composites  for  all  orientations. 


f.  Determination  of  Stresses  from  Strains  for  a Given  Material 


Q . 
ij 


for  Gr/Ep  in  GPa 


181.811  2.897  0 

2.897  10.346  0 

0 0 7.170 


1 

» 

1 

{ 

J' 

X 

, tr'  <r' 

y xv 

1 J 

' 

— 1 

Figure  2 

1 1 1 

9 

.1 

F 

r 

or  a given  set  of  Imposed  Strains,  shown  are  the  resulting  stresses  for 
-300/5208  unidirectional  composites  for  *11  orientations. 
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2.  MAXIMUM  STRAIN  CRITERION 

a.  Determination  of  Ultimate  Strains 


b.  Typical  Ultimate  Strains  of  Unidirectional  Composites 


TABLE  15.  ULTIMATE  STRAINS  OF  UNIDIRECTIONAL  COMPOSITES  (mm/m) 


— 

Material 

eL 

eV 

. 

ei 

€ ' 

T 

CS 

v 

IT 

iM| 

mmm 

T-300/5208 

8.30 

8. 

30 

3. 

91 

23.  8 

9.42 

.28 

.0159 

B/5505 

6.  18 

12. 

2 

3. 

30 

1C 

.9 

11.6 

.23 

.0208 

27.  5 

15. 

8 

3. 

80 

14.2 

17.4 

.26 

.0557 

. 

37 


c.  Maximum  Strains  Criterion  in  Strain  Space 


Figure  31.  Three  views  of  maximum 
strain  envelope  in  strain 
space. 


Figure  32  Longitudinal  and  traneverse 

strain  in  uniaxial  axial  tension 
and  compression  tests. 


H vLT*L^  eT>  strains  envelope  may 
not  include  all  failure  points.  Point 
A,  for  example,  could  be  outside  fail- 
ure surface  as  shown  in  Figure  31  . 
This  can  be  an  inconsistency  of  this 
failure  criterion. 


e.  Maximum  Strain  Criteria  in  Stress  Space 

From  Stress-Strain  relation  S..<r.  = e. 

lJ  J 1 


ai  ' VLT<r2  - ELeL 


X 


-o'  + v o'  ^ X* 

1 LT  2 — 


■vTL<ri  + a2  ~ Y 


VTL<ri  ' *2  ^ Y' 


a6  ~ S 


(56) 

(57) 

(58) 

(59) 

(60) 


& 


Figure  33  Strain  criterion  is  limited  to  the 
material  symmetry  axes  shown 
here. 
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f.  Off-axis  Uniaxial  Tests 


(1)  Maximum  Strain  Criterion 
Substitute 

, 2 , 2 

<r'^  = n or ^ = m 

eL-  (n2’  VLTm2)  ai/EL 
2 2 

eT~  (m  ' vTLn  ' VET 


o'*  = mnir.  into  SSR 

b 1 

(61) 

(62) 


(2)  Maximum  Stress  Criterion 

A special  case  of  (1)  when  v = v = 0. 


3.  PLANE  ELASTICITY  SOLUTIONS 


b.  Stress  Concentrations  at  Various  Boundary  Conditions 
(1)  Elliptic  Hole  Under  Tension 


(5)  Circular  Hole  Subjected  to  Twist 


k- D +n(cos  6+nsin  6) 


(2)  Elliptic  Hole  Under  Bending 


Under  Tension 


Circular  Hole  Under  Internal  Pressure 


At  Point  B 


Point  A 


Point  B 


Circular  Hole  Under  Uniform  Tension  General  Reference 


S.  G.  LEKHNITSKII,  Anisotropic  Plates  , 
English  translation  by  Gordon  & Breach 
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c.  Numerical  Example  of  Elliptic  Hole  Under  Tension 


At  Point  B 


Figure  36 


Elliptical  opening  fransvers 
to  the  fibe  s. 


composite 


tss  concentration  factor 
various  aspect  ratio. 
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4.  TRANSFORMATION  RELATIONS  FOR  COMPLIANCE 


a.  Formulas 


TABLE  17  TRANSFORMATION  OF  COMPLIANCE 


h ‘ i(Sll  4 S22-  ZS124  S66> 


R1  - r JT 


SH  * S22|Z  * (S!64  S26'2 


*2 

8 ^<S11+  S22  ‘ ZS12  " S66)  + 4(S26  ' Sl6’ 

e c 

tan  26, 

16  26 

1 

S11  * S22 

2(S16-  S26> 

tan  46^ 

S11+S22'2S12'S66 

For  orthotropic  material 
61  " 62 


X1 

SH+  S22+  2S12>  (82) 

(83) 

(84) 

(85) 

(86) 

(87) 

(88) 


For  anisotropic  material 

{ 

‘l  * 62 

in?) 

rrri 
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b.  Typed  Numerical  Data 


TABLE  18  INVARIANTS  FOR  COMPLIANCE  MATRIX  FOR  VARIOUS  COMPOSITES 
(TPa)’  1 


~ 

Material 

!is 

T- 

300/5208 

24.88 

B/550 

5 

14.  18 

Scotchply  1002 

33.  3 

Kevlar  49 

46.51 

1009-26-5901 

12.  30 

uus. 

37. 79 

30.65 

29.24 

50.2 

79.84 

18.41 

40.77 


R1S 

R2S 

45.78 

4.22 

24.55 

13.94 

47.5 

84.33 

10.2 

28.86 

11.92 

2.  37 

71.26 

1 . 89 

; 

(3)  Variation  in  Engineering  Constants  as  Coordinates  Rotate  (GPai 


I 


i 


1 1 


'll 


'22 


22 


12 


'21  * 1 S 


22 


1 


TABLE  19  TRANSFORMED  COMPLIANCE  & ENGINEERING  CONSTANTS  FOR  VARIOUS 
COMPOSITES  (TPal*1  AND  (GPa) 


T- 

300  / 5208 

B(  41  / 

5505 

S 1 

Engineering 

S.  ' 

Engineering 

IJ 

Constants 

Constants 

5.52 

-1.55 

0 

Eu 

= 

181 

4.90 

-1.13 

0 

E11 

= 

204 

E22 

= 

10.  3 

-1.13 

54.05 

0 

E22 

_ 

18.5 

0° 

-1.55 

97.09 

0 

1 12 

= 

.28 

1 1 2 

= 

.23 

1 21 

= 

.016 

0 

0 

172.71 

‘■21 ' 

.02 

0 

0 

139.47 

C12 

7.2 

G 1 2 

5.79 

13.77 

-3.66 

30.20 

El, 

= 

72.63 

15.  16 

-8.10 

36.43 

El. 

= 

65.95 

E22 

= 

10.75 

E22 

_ 

17.  32 

15° 

-3.66 

93.06 

15.58 

V12 

= 

.27 

-8.  10 

57.73 

-11.85 

v12 

= 

.53 

vz\ 

= 

.04 

‘ 2 1 

= 

. 14 

30.20 

15.58 

131.03 

C12 

7.63 

36.43 

-11.85 

144.84 

°12 

6.90 

WM 

34.75 

-7.88 

46.96 

E11 

= 

28.78 

38.  10 

-22.03 

45.42 

E ’ 

= 

12.42 

OJ 

o 

o 

-7.88 

80.53 

32.  34 

.23 

-22.03 

62.67 

-2.86 

V 1 2 

— 

v22 

= 

. 10 

46.96 

32.34 

114.15 

G 

_ 

8.76 

45.42 

-2.86 

98.06 

I L 

59.75 

-9.99 

45.78 

Ell' 

= 

16.74 

E22’ 

= 

16.74 

45° 

-9.99 

59.75 

45.78 

V12 

= 

. 17 

V21 

= 

.17 

L 

45.73 

45.78 

105.71 

g|2' 

= 

9.46 
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5.  TRANSFORMATION  RELATIONS  FOR  MODULUS  MATRIX 


a.  Formulas 


TABLE  20  TRANSFORMATION  OF  MODULUS 


— 

T. 

X2 

R1 

R2 

Qil 

1 

1 

cos  2(6-  6 ) 

cos  4(6-  6^' 

°22 

1 

1 

-cos  2(6-  6j) 

cos  4(6-  6,1 

Qi2 

1 

-1 

0 

cos  4(6- 

°66 

0 

1 

0 

- cos  4(6-  62) 

°16 

0 

0 

sin  2(6-  6j ) 

- sin  4(6-  62) 

°26 

0 

0 

- j sin  2(6-  6^) 

sin  4(  6-  6.,' 

T1  = 4 <Q11  4 Q22+2Q12’ 


H ~ 8 <Q11  + °22  ' 2Q12  + 4Q66' 


(Q0) 

(911 


2 > ('Q11+  Q22)2+  4(Q16+  Q2b' 


,rQ22-2Q!2-  4Q66'2+16<Q16-  °26'2 


tan  26 


2(QI6+Q2^ 

Qll‘  °22 


, 4(Q16-  °26> 

tan  462  • + Q,  , + Q,,-  2Q,,-4Q 


11  22  12  66 
6^  = for  orthotropic  materials 


(92) 


(93) 


(94) 


(95) 


6j  ^ for  anisotropic  materials 


b.  Typical  Numerical  Data 


TABLE  21  INVARIANTS  FOR  MODULUS  MATRIX  FOR  VARIOUS  COMPOSITES  (GPa) 


Material 


IQ 


2Q 


R 


IQ 


R 


2Q 


T- 300/5208 

B/5505 

Scotchply 

1002 

Kevlar  49 

S-Glass 

1009-23-5901 

HMS  3002M 


49.49 
58.03 

13.0 

21.49 
24.76 

48.69 


26.88 

29.77 

7.47 

10.95 

15.47 

26.60 


85.73 

93.20 

15.4 
35.  55 
18.  35 

89.25 


19.71 

23.98 

3.  33 
8.65 
3.47 

20.74 
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d.  Typical  Data  for  Transformed  Modulus  for  Various  Composites 


TABLE  22  TRANSFORMED  MODULUS  FOR  VARIOUS  COMPOSITES  (GPat 


I 


T - 300  / 5208 

B / 3502 

GUI) 

ij 

181.81  2.90  0 

2.90  10.35  ~ 0 

0 0 7.17 

204.98  4.28  0 

4.28  19.59  0 

0 0 5.79 

q!15) 

1J 

160.47  12.75  -38.50 

12.75  11.98  -4.36 

-38.50  -4.36  -17.05 

180.51  16.27  -44.07 

10.27  19.08  -2.53 

-44.07  -2.53  17.78 

q(30) 

1J 

109.38  32.46  -54.19 

32.46  23.65  -20.05 

-54.19  -20.05  36.78 

122.41  40.25  -61.13 

40.25  29.21  -19.59 

-61.13  -19.59  41.78 

56.66  -42.32  -42.87 

63.82  52.24  -46.60 

„(  4 5) 

Q. . 

U 

-42.32  56.66  -42.87 

42.87  -42.87  46.59 

52.24  63.82  -46.60 

-46.60  -46.60  53.76 

9 

! 

• 
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6.  MODULI  FOR  RANDOM  COMPOSITES 

a.  Constant  Stress  or  Series  Model 


"72 

I 

- */2 


sij  de 


Il+12  VX2  0 


h+h  0 


The  invariants  are  those  for  S..  in  Table  18. 

lJ 


's  - I1+!2 


8 ' ri+I2 


s 4I„ 


b.  Constant  Strain  or  Parallel  Model 

n/2 


h+l2  h-h  ° 


■i  / °. 


. de  = 


h+h  0 


h J Ve6 


The  invariants  are  those  for  Q..  in  Table  21. 

yi2 

= Xl+I2 


Ep  = d-"  >w 

Gp  * >2 

c.  Isotropic  Constants  for  Various  Random  Composites 

TABLE  2 3 PREDICTED  ELASTIC  CONSTANTS  FOR  RANDOM  COMPOSITES  (GPa) 
Material  E G v E G v 


Material 

E 

8 

T- 300 / 5208 

18. 

B/5505 

23.0 

1002 

12.0 

E 

P 

G 

P 

V 

p 

69.7 

26. 

.295 

78.7 

29. 

B 

.320 

19.0 

15. 

i 

.270 
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SECTION  IV 


MICROMECHANICS 


1.  ELASTIC  PROPERTIES 


a.  Definitions 


e. . = 7 / e..dV  = -777  f (u.  .+ u.  )dV 

D V J ij  2V  J i,j  J,i 


= Tv  V 

u.n.+  u.n.)dS 

11 

(104) 

s 

j j 

Average  stresses  depend  on  the  boundary  tractions 


s = sl+s2 


V = v1+v2 


oniy. 

] 

Figure  38  . 

A composite  body,  Di 

1 

r 

1 r 

placement  continuity 
and  equilibrium  are 
maintained  at  the  inter 
face. 

O'.  . 

V 

1 

0\ 

xj 

.dV 

) 

V 

7 

V 

(0 

X.), 
Ik  J 

kdV 

1 

r,r 

v 4-  T v Iris 

2\ 

r . 

/ 

s 

1 i 

j 

n* 

Virtual  work 

1 

f 

1 r 

or.  .e. . 

IJ  IJ 

- 

V 

J 

V 

o',  .e  . .d\ 
ij  *J 

r 

- 

V 

4 

T 

.u.< 
1 1 

lb 

(106) 

O’  • 

Qta Hrallv  aHmi  sflihlp  t 

? • 

lei  npmaHralW 

admissible 

ij" 

7 

'ij* 

7- 

Divergence 

theorem 

c 

r 

4 * 

n.dS 

J 

= 

j 

v 

f. 

] 

.dV 
'< J 

(107) 

V 

, 

I 1 4. i 1 — r — • — ♦ — 

! 

1 

1 i_L  1 1_L 
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(1)  Sample  Problems 


I 


i 


I 


C11  = zV  PuinidS 

s 

'll  ■ “v  / |uri»ds  ■ 7a  6ia  ■ T - 6 

(b)  A composite  bar  is  subjected  to  the  » x 


pressure 

p at  both  ends 

. What  is 

IUC 

avci age 

oliCO  9 

^ ^ 

<rll  = 

"v  /V 

,ds  « P 

1 X 

s 1 1 
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b.  Homogeneous  Boundary  Conditions 


If  “i  = eijxj  on  S'  then 
Tj  = <rijnj  on  S,  then 
If  either  of  the  above,  then 


e. . 
V) 


cr. . = cr. 


ij 


<r..e. . 


<r.  .e. . 

ij  U 


The  above  reduction; 

s follow  from  the  definitions  and  the  divergence  theorem 

e. . a - 
ij  < 

n.  + u,n.)dS 

J J i 

> 

1 Co 

: 

= [€i?k/XknjdS+€J?k/XknidSl 

b 


o',  .e. 
!J 


V J T.u.ds 

V S 11 


_L  f ° 

- IT  J <r.;;n,u; 


dS 


<r? 

1 


H-  /»,», 


J i 


dS 


(108) 

(109) 

(110) 
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(I)  Sample  Problem 


(a)  Prescribe  the  homogeneous  boundary  conditions  on  a two-dimensional, 
rectangular  body. 


c.  Effective  Elastic  Constants 


Displacement  B.C.  u.  = e°.x.  on  S. 

1 D J 

u.(x)  = u.kf(x)«°(  inV 

eij(x)  = Wkf  ' «ijkf  = I(uikf,j+ujkf.i) 

Vx)  = cijkf(x,ew(x) 

7 = c*-k/*kf  ■ c*k  = V f C - (x)e  k!(x)dv 
1J  ljki  kf  ljkf  V J ljmn  mnki 

Traction  B.C.  T.  = a°.n. 

i lJ  J 

= S*  7 , S*  - Tr  l S--  (x)a  t ,(x)dV 

ij  ljkf  kf  ljkf  V Jy  ljrrm  mnki 


= a 

ijkiv 

,<rki 

Strain  energy 

w‘ 

1 - - c 

< 

0 0 1 r-* 

_ 

- _ 

Displacement 

B.C 

" 2 e 

ijki  ij  kf  2 ij 

kf € 

*ijekf 

wa 

' * 

t 

jkf 

O 

a. . 

° _ i c* 

kf  2 ijk 

tai 

jaw  : 

Traction  I 

s.c 

• 

Sta 

fi  fltipsIlM 

■ V.  <i 

trypr>  a An  4 Krt/itf 

Ola 

iisiicaiiy 

11 U 

igcncuuD  uuuy 

l 

c;.,  .sr. 

ijk/  kimn 

9 

: 1 

, i.. 

ljmn  ljmn 

= 

J (&im 
c.  inn 

&.  + & 
Jn 

• 6.  ) 
in  jm 

d.  Prediction  Based  on  Averages 


Determine:  e.  = S.  .a.  (6  equations) 

1 ij  j 


(120) 


Unknowns  : <r.,e,  a e , <r  ,e 

ti  fi  fi  mi  mi 


(36) 


(121) 


Available  Equations: 


Const.  Eqns.  e = S„  .a..  , e r S <r 

fi  11.1  fj  mi  mij  mj 


Definitions 


e.  r v . e + v e . , <r.  = v.<r_.  + v a 

l f fi  m mi  i f fi  m mi 


( 12  equations)  ( 122) 

(12  equations)  (123) 


No.  of  additional  equations  required  : 6 


To  be  provided  from  the  equilibrium  and  compatibility  conditions. 


Suppose 


J 


(124) 


Then 


VfSfikBfkj+ Smik(6kj"  VfBfkj' 


(125) 
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e.  Approximation  for  Unidirectional  Laminae 


For  equilibrium 


<r  ~ -n . V , i = 2-  6 . 

mi  'i  fi 


For  compatibility 
6fl  = 11  l^ml 


# 1 

E,  - (n  v.E,+  v E ) 

ri . v + v 'iff  mm 

'If  m 


* 1 

v_,  = ; (n.^.v.+  v v ) 

12  ' ti,v  +v  Iff  mm 

if  m 


(1281  Figure  39  Unidirectional  composite  and 
the  reference  coordinate  sys- 
tem. 


± JZi  , 1 

m*  E / vr+9 


-+T1  V 

f am 


(r|  ri  E v -E  v )(v  /E  -v  /E.) 
1 a f m m f m m f f 

f m (ri  v.E  + v E )(v.+ri  v ) 

Iff  m m f a m 


v = v E / E, 
al  la  a 1 


( Vf  ”cVm  \ 

1 

K Gm  ’ 

vf+TVv™ 

i cm 

K ( 

1 + V. 

3 ( V 

f 

V.+  Tl  V \ f 

E + 

f am 

f 

= 2,  3 c 

II 

4^ 

O' 

ti  v -rr* 

a m E 


Note  that 


* * 

V V 

23  32 

* ^ * 
E2  E3 


— — unless  n. 


f.  Transverse  Isotropy 


TT*  - IT* 

( l 351 

E2  ' E3 

^2  = 

”3 

* 

*2  3 

1 

1 

E* 

^2 

E2 

2G4% 

— ^2  = ^4 

( 1 36! 

A # 

C44  = G55  

^4  = 

^5 

( 1 371 

Independent 

T>  '6  ! 

v ”2' 

’ ^ 6 

Independent 

elastic  moduli  : 

el 

Et 

(1381 

et  = 

E2  = 

E3 

(1391 

II 

H 

v E 
TL  L 

/et  = 

* 

*12 

(1401 

glt  = 

r* 

G66 

(141) 

‘tt 

* 

*23 

* 

V 

32 

(142) 

gtt(=  « 

J44  r 

G*^)  is  determined  from 

(143) 

et 

G „ = 



(144) 

TT 

2(1+  Kj.t) 
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Figure  40  Longitudinal  modulus  can  be 
predicted  by  Eq.  (128)  with 


1 


1.  Constituents  are 


E-glass  fibers  (E^=  72.4  GPa, 

*^=  0.20)  and  epoxy  matrix 

(E  = 3.45  GPa,  v = 0.  35). 
m m 


10. 5msi 
4.  34 
.2 

0.  185msi 
.35 

0. 185msi 


Figure  41  Prediction  of  transverse  modulus 
by  Eq.  (130)  with  q ^ = 0.5,  h = 1, 
is  comparable  to  the  more  rigorous 
solutions.  Constituents  are 

E-glass  fibers  (E^=  73.1  GPa, 


v = 0.22)  and  epoxy  matrix 

(Em  « 

m 
[2.3] 


(E  = 3.45  GPa,  u = 0.35). 
m m 


10. 5msi 

4.38 

0.2 

0. 5msi 


II  II 

6.-> 

1 

G : 

m 

c 0.185msi 

0.6 

V 3 

m 

= 0. 35 

f 

1 

",  1 

[ . -05 

vf  ' 

! : J I J 

f i 0.6,  . _i 4 . . . 

J— 1 . 1 J L 1 J. 

L.  1 1 1 i 1 _ . 1 ! L . 
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MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BURFAb  OF  STANDARDS  I96i  A 


Long.  Modulus  E , 10  psi 


Figure  40  Longitudinal  modulus  can  be  Figure  41  Prediction  of  transverse  modulus 


predicted  by  Eq.  (128)  with  by  Eq.  (130)  with  r| 2 = 0.5,  q = 1, 

Tij  - 1.  Constituents  are  Is  comparable  to  the  more  rigorous 
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Long.  Poisson's  Ratio  v 


1.45  0.50  0.55  0.60  0.65  0.70  0.75  0.80 

Fiber  Volume  Fraction  v 


Longitudinal  shear  modulus  increases  with  decreasing 
by  Eq.  (132).  Good  agreement  is  seen  whenri,  = 0.4 


Fiber  Volume  Fraction  v 


Fiber  Volume  Fraction  v 


Figure  43  Longitudinal  Poisson's  ratio  ca 
be  predicted  by  Eq.  (129  ) with 


Figure  44  Transverse  Poisson's  ratio  pre 


dieted  by  Eq.  (133  ) withfi.,=  0.5 


1.  Constituent  properties 


is  higher  than  the  upper 


the  same  as  those  in  Figure 


bound  from  a more  rigorous 


solution 


g.  Approximations  for  Particulate  Composite 
Assume 

v . - Ti.  TT i -=  1-3;  ~5  . * n,  or,., 

mi  '1  fi  mi  '2  fi 


_1_ 

/ 

Vf 

) 1 

E* 

= l 

Ef 

+ E 

m 

/ vf+TlVm 

r „ k 

£ + r»v 

i/  E. 

* 

V 

= - 

f f 

v I 

m 1 m m f 

r 4-  ~ ir  IT 

f 

I I * 

m ' m 

f 

Isotropy  requires  that 

^1 

= 

^2  = ' 

n 

Masnmand  fruikmn-^y  / A* 
Upper  Bound  // 

/y  vzf' 


^u''i  Appronmaie  / 

5 Solution  / 

U s'  Ishai's  Approximate 

3 / Solution 

■“  / / 

■f  and  Shtnkmen 

| / lower  Bound 

£.  q =0.5 

..  i ■ 

0 0»  OKI  0)0  0.<0  0 50 

Volume  fraction  of  AI^Oj  dispersion  m glass  matri* 


ftul’s  Approumaie 
Solution  , 


Ishai's  Appronmate 
Solution 


6 0.20  0 40  0 60  0 80 

Volume  fraction  of  sand  in  epoxy  matrix 


Figure  45  Youngs  modulus  of  particulate  composite  can  be  predicted  by  Eq.  (146  ) with  n = 0.4- 
0,5  : (a)  Al^O^/glass  (Ej=s  396  GPa,  E^=79.3GPa);  (b)  Sand/epoxy 
(E,®1 124  GPa  , E *«  3. 45  GPa) . [5] 

i m 


a) 

Ef 

E 

_ 

57. 

,4  msi 
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E, 

18  msi 

_ 

11. 

, 5 msi 

I 

E 

0.5  msi 
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h.  Thermal  Expansion  and  Swelling  Coefficients 


Hygro- thermo- elastic  relations  for  homogeneous  constituent  materials 

e.  = S.  ,<r.  + a.  6 or  (148) 

j ij  j 1 

<r.  = C..e.+  p.  0 

i ij  3 1 

where  0 is  the  temperature  difference  (T)  or  the  moisture  concentration  c . If  both,  a.& 

T H 1 

should  be  replaced  by  (a.  T + a.  c). 

(j  Q 

Let  e.  and  e.  be  the  strains  resulting  respectively  from  the  homogeneous  stress  B.C.  with 
0=0  and  the  stress  free  B.C.  with  0 = 0°. 


— <r  o _*  _ 0 

e.  = S . . it.  = a. . <r.,  e. 

i V)  3 l3  3 1 


Therefore 

* _ * — 

e.  = S. . a.  t o.  6 . 

I ij  j i 


Similarly 

a.  = C*  e.  + p*  0 

i i]  3 1 


(1)  Sample  Problem  [2] 
(a)  Show  that 

* 


a .+(a,.-a  ,)S,  (S  - S ) 

mi  fj  my  fmjk'  ki  mki' 


where 


fmij  Sfij  ^mij 


Proof 

From  the  virtual  work  theorem 

~ 5"  — ~e  _ * — 

"i  ei  = Vi  = Vi6 


O' 

<j 

9 

e 

s. 

I 

c 

V 

r ® , 

(T  4- 

O' 

9 

(155) 

i 

i 

W . T 

J 

O'. 

1 

a. 

i 

= 

<r  € 

€.  c r. 
1 J 

1 , * 

+ o',  a. 

i i 

e 

= 

0 

(T 

ri  ai 

9 

(156) 

T 

’he 

reform 

o-.c 

i 

L, 

l 

= 

V 

r 

a 
n i 

a 

.<r 

tu  mi 

+ 

vf' 

O' 

a£i<rfi 

(157) 

Let 

— <r 

tr . = 

: ] 

B 

t j 

. T1 

fi 

fij 

j 

» 

(158) 

= 

g 

L + 

Vsn, 

-s 

B. 

O] 

p 

mil 

r 

mj 

2S 

* 
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= 

Sj 

-1 

k {S\ 
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/ 1 Am 

IlJ 

:mi 

‘j 

mK] 

\ I 
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O' 
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V, 

O' 
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we 
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mi 
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(161) 

oince 

L 

is 

arbitrar- 

y,  the  proof  follows. 
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EFFECTIVE  EXPANSION  COEFFICIENT, 


i.  Approximations  for  a. 


Following  the  same  procedure  as  in  Section  e,  we  obtain 

°L  = °1  = f+V^m^SWV^m* 


(162 


= Vf+  vmam  + (Vf  vt+  VmQm  ‘'m''  (t1  lvf  lf+  Vm  *'m)al 


(16: 


Figure  46  Thermal  expansion  coefficents 
predicted  by  Eqs.  (162-16  3)  with 


— 

s 1 arc  comparable  to  the 

i _ roi 

more  rigorous  solutions*  [£•] 

TABLE  24  TYPICAL  COMPOSITE  THERMAL  EXPANSION  COEFFICIENTS 


1 


r 

LONG 

TRANS. 

COMPOSITE 

CONSTITUENTS 

vf 

EXP. 

EXP. 

(im/m/,K 

(»m/m/  “K 

B/Ep 

B ./5505 

4 

0.  5 

4.  32 

22.  1 

B/PI 

B . /WRD  9371 

0.49 

4.90 

28.4 

Gr/Ep 

Mod  n/5^06 

0.55 

-0.  23 

34.0 

i 

or/ii^p 

HM5/ 3UUCIVL 

0.  48 

-0.  23 

33.  5 

Gr/Ep 

T300/5208 

0.70 

0.01 

12.  5 

! 

ur/n,p 

Mod  I/EJxvi-iA 

0.  51 

-1.  10 

31.  5 

Gr/Ep 

Gr /PI 

Mod  I/ERLA  4617 

Mod  I/WRD  9*^71 

0.45 

0.45 

-0.90 

0. 

33.  3 

25.3 

Gl/Ep 

S-Glass/1009 

0.72 

3.8 

16.7 

Gl/Ep 

Scotchply  1002 

0.45 

4.16 

15.5 

2.  STRESSES  IN  CONSTITUENT  PHASES 


i 

a.  Microstresses  Under  Longitudinal  Tension 


Figure  47  Hexagonal  array  of  fibers  and  the  coordinate  system  chosen. 
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b.  Stress  Concentration  Under  Transverse  Tension 


Figure  49  Stress  concentration  under  transverse  tension.  Solid  and  broken  lines 
represent  the  stress  concentration  factors  at  the  interface  and  in  the 


matrix,  respectively.  [7,  8] 

c.  Stress  Concentration  Under  Longitudinal  Shear 


Figure  50  Stress  concentration  under  longitudinal  shear.  . Solid  and  broken  lines  represent 
the  stress  concentration  factors  at  the  interface  and  in  the  matrix,  respectively. 

[7.61 
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ORIENTATION,  6,  DEG 


Residual  stresses  at  the  interface  in  a composite  with  hexagc 
Constituent  properties  are  the  same  as  those  in  Figure  48  , 
correspond  to  the  unconstrained  matrix  shrinkage  strain  of  1 


d. 

Residual 

Stresses  After  Curing 

TABLE  25  TOTAL  STRESSES  LNDER  LONGITUDINAL  TENSION 

<cf.  [9]) 


Fiber  Volume 

Content 

vf 

Modulus 

Ratio 

E,  / E 
f m 

Radial  Stress,  o 

Iksil r 

Hoop  Stress,  a 
(ksi)  9 

8 = 30 

Maximum  Shear 
Stress,  r . 

(x.i)  r8 

0 = 0° 

e = 3o° 

Curing  Stresses  (Resin  shrinkage  = 0.6%) 

0.  64 

150 

-1.2 

0.  3 

3.0 

0.6 

0.70 

150 

-1.2 

0.  6 

3.6 

0.6 

0.64 

26 

-1.5 

0.  3 

3.0 

0.6 

0.70 

26 

-1.8 

0.6 

3.6 

0.6 

Stresses 

Due  to  Longit 

udinal  Tensio 

n of  100  ksi 

0.64 

150 

-0.7 

0.4 

7.5 

1.5 

0.70 

150 

-1.2 

1.0 

9.0 

2.5 

0.64 

26 

-1.0 

1.0 

7.5 

2.0 

0.70 

26 

-1.3 

1.3 

9.0 

2.5 

Combine 

d Stresses 

0.64 

150 

-1.7 

0.7 

10.5 

2.1 

0.70 

150 

-2.4 

1.6 

12.6 

3.1 

0.64 

26 

-2.5 

1.3 

10.5 

2.6 

0.70 

26 

-3.1 

1.9 

12.6 

3.1 
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3.  STRENGTH 


i 

i 


a.  Longitudinal  Tensile  Strength 
Brittle  matrix 


Influencing  factors 

Fiber  damage  - fiber  surface  treatment  (Gr  fibers),  heat  treatment  (B/Al) 

Strength  of  B^'Al-6061  (v^  = 0.25) 

Heat  treatment  No  treatment 

98  MPa  104  MPa 


Figure  j 

52 

Longitudinal  tensile  strengths  of  B /Al-6061  and  BSC/Al-6061.  [10] 


TABLE  26  TYPICAL  FIBER  STRENGTHS 


MATERIAL 

Kevlar 

Graphite  (Gr) 

MODMOR 

II 

HTS 

MODMOR 

I 

HMS 

T300 

AS 

Glass  (Gl) 

E- Glass 

S-Glass 

Boron  ( B) 

Borsic  ( BSC) 

Steel 

Tungsten  (W 

') 

Beryllium  (Be) 

DIAMETER 

nm 


TENSILE  STRENGTH 
MPa 


.7 

2760 

.6 

27  60 

1850 

2240 

,_L ! --L4-H  i 
■ ■ ■ i l ■ 


«•  » 

^ 4 


TABLE  27  TYPICAL  MATRIX  STRENGTHS 

material 

TENSILE 

COMPRESSIVE 

SHEAR 

STRENGTH 

STRENGTH 

STRENGTH 

MPa 

MPa 

MPa 

Epoxy  (Ep) 

Narmco  2387 

29 

159 

10 

Narmco 

5505 

59 

128 

Epon 

1004 

B28 

72 

82 

150 

207 

83 

134 

ERLA 

4289 

34 

93 

ERLA 

4617 

132 

226 

Polyeste 

r 

72 

Aluminum  (Al) 

2024  - T3 

427 

221 

255 

6061 

T 

6 

290 

241 

186 

6061 

0 

131 

- 

83 

Titanium 

(T:) 

6A1-4V 

958 

951 

558 

Pure 

552 

483 

290 

TABLE  28  TYPICAL  COMPOSITE  STRENGTHS  [11,12] 


■ 

COMPOSITE 

CONSTITUENTS 

Vf 

LONG. 

TENS. 

MPa 

LONG. 

COMP. 

MPa 

TRANS. 

TENS. 

MPa 

TRANS. 

COMP. 

MPa 

SHEAR 

MPa 

B/Ep 

B . / 5505 

0.5 

1260 

2500 

61 

202 

67 

B/PI 

B . /WRD  9371 

0.49 

1040 

1090 

11 

63 

26 

Gr/Ep 

4 

Mod  11/5206 

0.  55 

1110 

970 

36 

170 

63 

Gr/Ep 

HMS/3002M 

0.48 

680 

690 

16 

186 

72 

Gr/Ep 

T300/5208 

0.70 

1500 

1500 

40 

246 

68 

Gr/Ep 

Mod  I/ERLA  4617 

0.45 

840 

880 

42 

197 

62 

Gr /PI 

Mod  I/WRD  9371 

0.45 

807 

652 

15 

71 

22 

Gl/Ep 

S-Glass/1009 

0.72 

1290 

822 

46 

174 

45 

B/Al 

B4/AI-606I-F 

0.5 

1110 

1480 

103 

159 

103 

BSC/Ti  BSC  / Ti-  6A 1 - 4 V 0.5  1310  - 441  - 483 

Gr/Ep  Mod  1/4289  0.5  1120  990  4.  2 - 34 


Gl/Ep  Scotchply  1002  0.  45  1062  610  31.4  118  72 


b.  Transverse  Tensile  Strength 


Failure  occurs  when 


K t = X , 
T m m 


K^,  : effective  stress  concentration  factor 


Y = 


1 


Tm 
1 


(v  + v,/n0)X  : matrix  fail  ire 

m i £ m m 


„ (v  .+  v./T1,)X.^X. 
K^,i  m f '2'  1 i 


interface  failure 


(166) 


(167) 


2.000 


1,000 


I 

■ 


! 

! 

1 

1 

I 


) 

i 


Figure  53  Transverse  strengths  of  B^/Al-6061  and  BSC/Al-606l  [10] 


Influencing 

factors 

Fiber  splitting  decreases  Y 

b4/ai 

Heat  treatment  increases  X and  in  turn  Y . 

m 

Voids  increase  K and  in  turn  decrease  Y . 

Fiber  surface 

treatment  increases  interface  strength  and  in  turn  Y . 

j I J 1 | 

1 1 1 1 * 1 

f j — f ‘ 11 

_ -j  1 — |— f f r — i j — 1 ; 

- J L I . J 

— 1 — f 1 — — — pi — r — — t — ' — f — 1 — : 

i.-l  L 1 11  1 11  11  M 1 1 
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Kg/cm‘ 


Figure  54  Effect  on  transverse  strength  of  fiber  diameter  and  matrix  condition  (B/Al) 
Failure  of  B /Al-6061  is  due  to  fiber  splitting.  [10] 
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Figure  57  Failure  mode  in  flexure  test  depends  on  the  shear  strength 

(Gr/Ep,  = 6.6):  U - fiber  surface  untreated:  T - fiber  surface 


treated. 

[14] 

- I j 1 . | J • | ..  - f 

d.  Longitudinal  Compressive  Strength 


(1)  Buckling  in  Shear  Mode  - , 

G xv  / d \ 

x'  - 77^7  + * -nr-  Ef  \-t)  ri5i 

L : gage  length 


(170) 


1 for  simply  supported  ends 
4 for  fixed  ends 

(2)  Actual  Composites 
G , 

— ■ -m  - *»  — GPa  : too  high'  (171) 

1-  v,  1-  v. 
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4.  SYNERGISTIC  EFFECTS  OF  MATRIX  AND  FIBER  STRENGTH  SCATTER 

a.  Role  of  Matrix 

Minimizes  interaction  among  fiber 
breakages  - beneficial. 

Induces  stress  concentration  and 
hence  possible  fracture  - deleterious. 


Figure  58  Fiber  breakages  in  composite. 


b.  Bundle  Strength 

Fiber  strength  in  stress 

Rf(Xf)  = exp  r-L(Xf/Xfo)a],  Xf  = XfQL' 1 /ar<  1 + 1 /a ) 


(175) 


L.inea.r  Rtrp.Rf 

1-st.rain  rolatinn 

O = Efc 

(176) 

Fiber  strength  in  strain 


R,( 

Y,) 

1 = exD  T-LiY./Y. 

*1. 

Y 

if 

/E 

/ 1 *7 

r 

r 

‘ ‘ ' f fo'  r 

fo 

fo 

f 

\ * « « 7 

Fraction  of  remaining  fibers  at  e 

= Rf(c). 

Gross  s 

i tress  of  bundle 

a 

= i 

ELeRJe) 

1 1 

f F 

da 

in  — 
dc 

Failure 

occurs  whe 

= 0. 

* * B 

YR 

= 

Yf  (La) 

rw“ 

(179) 

Bundle  strenffth 

x_ 

_ 

E.Y_RJY_) 

= x,  1 

(Lfl 

t e] 

r1J 

<0t 

_ 

X,( 

,-l 

ae) 

/a 

/rn  + i/ai 

(180) 

D 

l B r 

B 

fo 

r 

Bundle 

secant  modulus 

l i — * 

. 

a — 

ft  - 

eb 

s: 

7 - 

= Ef 

expr-1 

/Yl 

)o> 

j . 

(181) 
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Figure  59  Normalized  bundle  strength  vs.  the  inverse  of  shape  parameter.  [16] 

c.  Composite  Strength  [17] 

X - v ( strength  of  bundle  of  length  6) 

= VfXB(L/6)'1/a 
= vfXfo(6aef1/“ 

d.  Determination  of  6 


Elastic  shear  lag  analysis  [17 J 

r 

i 

°f  ^f  j cosh  7?(X-x 

a E " , — 

cosh  rj* 


T 1 1 f sinh 


cosh  rjt 

2 

sz 

/ 

G V' 

m 1 — 

n 

= 4 

, . < sr 

1/2  \ 

in  1 

J7  'X 

s x/d  , 1 = l/df 

1 

l w 

1V\ 

r ” * 

) 

7 

f 

afdx 

= | 

Jf|! 

[l-b/2)  ► C 

i/2  = 1/r;  as  1 -»  • . 

•'o 

K = Z 

l 

1 -VV 

1/2 

• 

e 

_ 

_6_ 

J 

df 

?.  G 

m V^T 

F 

’erfectly  plastic  shear  lag 

analysis 

Xffd,/4 

= 

t trd  6/2 

y f 

T : matrix 

yield  stress 

(189) 

7 

• 

• • 

6 = 

*f 

/(2Ty) 

(190) 
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Composite 


f 

mm 


Xfo 

Ef 

Gm 

MPa 

GPa 

GPa 

3323 

79.  3 

0.024 

4413 

86.2 

0.024 

0 1 87 

372 

3994 

0.024 

6(  elastic) 
mm 


> 

1,  MPa 

Theory 

Exp 

R.M.  * 

144 

71.7 

150 

755 

145 

681 

215 

265 

230 

1656 

1207 

1365 

140 

132 

161 



E-glass 
Epon  815 
Epon  815 
S- glass 
Epon  815 
Epon  828 
Boron 
Epon  815 


6.20 


7.68 


11.11 


0.  127 


0.  127 


0.  102 


0.095 
0. 442 

0.095 
0.  565 

0.061 


7.73 

3.67 

8.06 

1.11 

15.68 


*Rule  of  Mixtures:  X = X,[v , + (l-v  )E  /E„]  . Data  from  [18]. 

f f l m f 

Remarks 

(1)  Stress  concentration  due  to  fiber  breakage 

(2)  Dispersion  of  failure  sites 
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Figure  60  Normalized  composite  strength  vs.  bundle  strength.  [16] 
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SECTION  V 
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IN-PLANE  PROPERTIES  OF  SYMMETRIC  LAMINATES 


(2)  Constant  In-Plane  Strain 

o (1)  ■ (n/2) 

'l  = €1  = •"  'l 


Figure  61  Ply  orientations  and  modulus 
of  symmetric  laminates. 
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(4)  Stress  Resultant 


N]  = h^i 


h (r_  ) or  N.  = h <r. 

2 ( l l 


where  N.  = Distributed  load  per  unit  width  of  a plate  with  thickness  h;  in  Nm"  or  Pam 


(5)  In- Plane  Modulus  A.,  (same  as  N.,  or  Nm  or  Pam) 


Substitute 


a(.t}  = Q(t)  e(t> 
i n i 


Q(t)  e° 


N.  = 2h  £ q!*’  e°  t 

1 ° tTl  lJ  3 


a/2 

= 2h  e°  V Q(t)  t 

° J t=l  « 


= A..  e° 
»J  J 


where  A..  = 2h  V*  Q. . t 

lJ  0 tTl  *3 


Since  most  practical  laminates  have  upto  4 ply  orientations, 


h Y Q(o)n 

o Z— f « 


°!  ’ °2’’ 


na  = total  number  of  plies  with  t orientation. 


Q. . * Modulus  of  a orientation. 

‘J 


is  governed  by  simple  rule  of  mixtures;  stacking  sequence  is  of  no  consequence. 
Stacking  sequence  is  critical  for  other  laminate  properties,  such  as  flexural  rigidity. 


I 

. 


b.  In-Plane  Modulus  and  Compliance 


f 

N, 

1 

N2 

i = 

.N6> 

A11  A12  A16 


A21  A22  A26 


r o > 
ei 


or  N. 
1 


A..  e° 
l3  1 


Let  Compliance  = a..  = [A.l  or  a..  A = 6. 

H ij  1 ij  jk 


ik 


(203) 


(204) 


r o^i 

e. 

1 

fa.. 

3.  n ^ | Z 

! 

r > 
Ni 

1 

1 1 

Id  10 

o 

4 21 

a a 

M 

0 

= a M 

(205) 

e2  | 

' 

22  26 

2 

i 

ij  j 

J | * 

1 

Ve6> 

*61 

1 

a62  a 66 J 

N, 

^ 6 J 

where 


wnere 

i 

(A 

l66 

-a! 

- 

a. 

_ 

(A 

l66 

-A2.) 

11 

A 

22 

22 

A 

11 

10 

A 

o 

a 1 7 

= 

A 

<A,AA 

26* 

A 1 2A 

66^ 

, 

a 

66 

= 

(A 

i/ 

k22 

-a121 

(206) 

A 

A 

k 

2 

ai* 

> 

- 

A 

(A 

12A 

26* 

-A. 

L 

!2A 

16'’ 

a, 

26 

“ 

A 

(A 

k 

12^ 

16 

-AnA 

l26‘ 

1 

’a 

A 

A " 

11 

12 

16 

i 

A 

= 

A 

21 

A 

22 

A26 

A 

A 

A, , 

l — 

■ 

61 

; - - 

62 

66  J 

. -J 

1 1 

In- 

Plane 

Encrineerine  Constants: 

| 

< 

1 

_ O 

1 

, 

o 

i 

1 12 

..o 

li 

| 

1 

E 

11 

= 

ha.  , 

9 

1 

E22 

hi 

» 

V 

12 

= 

“ ' 

•„ 

9 

*12 

= “ 

a2 

9 

0 

1 

Ai 

> = 

h; 

1 66 

1 1 

dd 

C 

Units: 

(mm 

- 

1 

/ 

j 

1 

8 . 

1 

m 

or 

**r  i 

i.s; 

! 

. 

i 

1 

k J 

[ “1 

' 

N. 

i 

(Nm 

or  Pam 

) 

ij 

LN 

m c 

ir(  Parr 

5 1 

] 

A 

/ XT. 

_ ] 

l 

_ n( 

k 

□ 

E 

,o 

/ M* 

-2 

30  C 

H»  J 

ij 

1 

\ Nl 

m 

Ol 

r 

MTI 

) 

□ 

■E 

'll 

:: 

t 



1 J 

h~ 1 

■ 

1 ] 

□ 

L 

l 

1 1 1 

1 

(207) 
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2.  FORMULAS  FOR  IN-PLANE  MODULUS 


a.  General  Multi-directional  Laminates 

A 

. ,/h 

= - Z Q<"'n 

(208) 

n » a2  lJ  n 

A 

ll/h 

= Ij+  I2  + VjRj+  v2r2 

(209) 

A. 

22 /h 

= I!+I2-  V1R1+V2R2 

(210) 

Al2/h  = Ij-I2  -V2R2  (211) 

A66/h  = I2  -V2R2  (212) 

Al6/h  = ' 2 V3RrV4R2  <213> 

A26/h  = -IV3R1+V4R2  (214) 


Ply  orientation  a is  opposite  to 
angle  of  rotation  for  trans- 
formation (used  inTbble  20). 


where 

= — — Z n„  cos  2or 
n ® 

1 'O 

(215) 

V2 

= — Z na  COS  4 or 

(216) 

1 O . ^ 

V3 

= > n-  sin  2 a 

n * 

(217) 

i 3 1 Oi 

V4 

- ~ L.  ***  sin  4a 

18) 

Since  A . transforms  the  same  way  as  Q..,  some  invariants  must  exist: 
1J  lJ 


*ii 

\ ~ 4 (A11+ A22+ 2A1 

2'  = hIl 

(219) 

1 -- 

lZJ 

\ = 8 'll  22"  1 

2"  ™66>  = nl2 

;)2+  4<A1( 

,>2 

R1A 

2 >|("A11+A22 

6 A2( 

V,  V2  -L  O 

(221) 

“ 1 1 3 1 

R» 

-L  J (A  + A„„- 

• 2A.  4 

iA66)Z 

+ 1 

6<A16 

- A 

_.»2 

2 

A 

2 i'll  22 

12 

Zb' 

(7531 

. f 2 2 

~ 

h «1  V,+  V R. 
12  4 i 

> 
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i 

• 

M 

| 


i 

1 

f 


I 


J 

< 


i 
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b.  In-Plane  Modulus  for  TO  /90  / 45  /-45  1 

1 p q r sJ 

Simplification  of  ply  composite's  factors  because 


a 

cos  2 a 

cos 

sin  2 9 

sin  4a 

0 

1 

1 

0 

0 

90 

-1 

1 

0 

0 

45 

0 

-1 

1 

0 

-45 

0 

-1 

-1 

0 

~ Ivv 


— ( n + n„  - n - n ) 

n o 90  45  -45' 


V3  = ‘T  <n45-n-45> 

v4  = 0 

r 

1 2 

2(  A 

16+A26>  n45*n-45  V3 

tar 

L 4 

1 All'A22  "o  ‘n90  V1 

2 Au*A22-2A12.4*66  V2  - 62  - »• 

1-  4 -4—~  -4 

(223) 

(224) 

(225) 

(226) 

(227) 

±2  n (228) 
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TABLE  29  FORMULAS  FOR  IN-PLANE  MODULUS  FOR  fO  /90  /45  /- 45  1 

L r>  n r e 1 


1 

L 

b 

J 

L_ 

L 

i 

r 

9) 

2 

2 

i r 

V1 

+ V 

3 

= 

n 

i 

(nc 

rn< 

J+  ( 

n4! 

rr 

i 

-45' 

(22' 

i 

V* 

+ V 

.2 

-Li 

IZ'U 

ru 

2 

4 

n 

o 

< 

)0 

4 

5 -4 

15' 

- n 

tan 

1 2S 

'l 

= 

-- 

-45 

1 

(231) 

o 

“90 

Laminate  is  o 

rthotropic  if  : 

s n _ 

-45 

• 

0 

1 

1 

(23< 

>\ 

IN  Ol 

vi 

r2 

o 

At  lower 

bound  (equal  to  zero 

>), 

the 

laminate 

! is 

is  i 

otropn 

C. 

o 

At 

upper 

1 bound  (equal  to  ur-ti 

f). 

the  laminate  ia 

i homogeneous , 

o 

In 

genera], 

laminate 

is  a 

ire 

always  less 

anisotropic 

than  1 

the 

constituent 

ply 

L 

— 

n 

. 



1 

. 

"1 

. 

h j 

1 

. 

' 

' 

1 

_i 

J 

L 1 

.U 

c.  Numerical  Example  of  In-Plane  Properties  of  T- 300/5208  Laminates.  For  this  material, 
the  ply  properties  in  GPa  are:  Ij  = 49.  5,  12  = 26.9,  R j = 85.  7,  and  1*2  = 19. 7 

TABLE  30  IN- PLANE  PROPERTIES  OF  T- 300/5208  LAMINATES 
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3.  CALCULATION  OF  PLY  STRESSES 


a.  Calculation  of  In- Plane  Strains 


c 

€. 

5 

= 

a. . 

N. 

+ a 

. ~N_+  a. 

6Nf 

] 

l 

U 1 

16  6 

1 

V 

*2 

= 

a2,Nl 

+ a 

22V 

a2 

6N( 

C . 

_ 

a . _ 

. N, 

+ a 

,,N_  + 

a. 

.N, 

f 1 

1 

) 

6] 

i 1 

62  2 

- 1 - ! 

6 

. 

6 f 

or 

r — f 

.J 

■ 

S 

n 

f 

ai  l 

1 

a 

1 7 

a,> 

(233) 


‘21 

l 

61 


22 

a 

62 


"26 

a66. 


b.  Calculation  of  Stresses  of  the  t-th  ply 


»(*>. 


‘11 


21 


L°61 


12 

22 

62 


16 


26 


66  J 


(234) 


(235) 


For  the  a oriented  p|yy  | replace  jndfex  ('t"  byi"aM  iii  the  above  equation 
I TAjBLE.31  SAMPLE  CALCULATIONS  OF  PLY  STRESSES 
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4.  PLANE  ELASTICITY  SOLUTIONS  FOR  SYMMETRIC  LAMINATES 


c.  TABLE  32  SAMPLE  CALCULATION  OF  ELLIPTIC  HOLE  UNDER  TENSION  IN  A 


T- 300/5208  ORTHOTROPIC  LAMINATE 
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SECTION  VI 


1.  FLEXURE  - CURVATURE  RELATIONS 


a.  Definitions  and  Assumptions 


(1)  Symmetric  laminates,  as  before  is 
assumed. 


Instead 


laminate 


varying  strain  across 


thickness  is  assumed 


Figure  64  Linear  strain  variation.  The 


is  curvature 


where  k 


(3)  Moments 


(239)  T 


[ure65  Ply  stress,  is  piece- wise 


r and  discontinuous 


distributed. moments  oz 


moments  per-unit  length  (N) 


b.  Flexural  Rigidity  and  Compliance 


i 

| 

* 

I 


. 


I 


y 


In-  plane  modulus  or  rigidity 


(244) 


Flexural  rigidity 

d..  = \ Lo(t)  (h3-  h3  ) 

lj  3 ij  t t- 1 


(245) 


A follows  a linear  rule  of  mixtures,  and  independent  of  stacking  sequence. 

ij 

D follows  a weighted  or  cubic  mixtures  relation,  and  is  highly  dependent  on 

ij 

stacking  sequence. 

Inverse  of  rigidity  D..  Is  compliance  d. . such  that 


k'l 

'dll 

A 

d12 

dl6 

k2  / = 

d21 

d22 

d26 

J 

-d61 

d62 

d66- 

M, 


M„ 


M , 


or  k.  = d .M. 
‘ U J 


(246) 


Unit  of  d„  is  (Nm) 


-1 


c.  Engineering  Constants 

Similar  to  those  for  in-plane  properties,  engineering  constants  for  flexural  properties 
must  also  be  derived  from  the  compliance  d not  the  rigidity  D_  because  the  latter 

can  be  related  to  engineering  constants  only  when  D, . is  orthotropic 

U 


Therefore 


J 

"ll 


12 


Id 


1 1 


h3d 


(El), 


11 


11 


£ 

.£ 

12 

/ 1TT\  - L 

( 

247) 

'22 

h3d 

'2 

d 

22 

| 

d 

12 

I-  • 

12 

‘ d 

lll 

! 1 

1 1 1 

■ 

i 

•y  a i 

r>\ 

J 

1 

2 

c 

J12 

h3d 

. 



66 

i l-  -J 

• 1 
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2.  FORMULA  FOR  FLEXURAL  RIGIDITY 


l 

i 

i 


i 

i 

i 


i 


i 


a.  Derivation 

For  a general  symmetric  structure  consisting  of 

(1)  A substructure  of  modulus  Q°  and  half  depth 

z , or  equivalent  plies  n (=  z /h  ). 
o o o o 

(2)  Two  symmetric  facing  sheets  (with  respect 
to  the  mid-plane  of  the  total  structure)  with 
composite  laminates  of  (n  /2)-  n plies  each. 


Rearrange  of  formula  for  rigidity 


i2_n  . t 

\3  3 

°)  r^o  J 2 \ 

1 n/2 

(249) 

(250) 

h3  ij  ' 'h 

where  F = 

Express  in  terms 

-)  % *(,N  ) 

t3-(t-l)3  = 

of  invariant- f 

Lj  * 

t=no+1  ^ 1 

3t2-  3t  + 1 

orm  transformation 

liD  Ho)  , 

(251) 

h3  1 l" V h / 

Du  - (I1+I; 

j)  + VjRj  + 

V2R2 

12  /2z  \3 

— D -( — ) i 

1 22  V h 1 

Q°  = (!,+  !,)  - V.R.  + 

V2R2 

(252) 

h ' ' * 

? ”..-W 

Qi2  = (Ir 

V 

V2R 

2 

(253) 

5 ■.-($*)’ 

°66  = lZ 

- 

V2R 

2 

(254) 

if  D - 

Q°  = 

- 4 V-R. 

V.f 

(255) 

h3  16  V h / 

16 

Z 3 1 4 

£ 

12.  /2Zo\3 

n 

1 „ 

3 26” \ h / Q26 

- ^ V,R 
2 3 

1 + V. 

4R2 

(256) 

where 

= (2)3  Lf 

(2 

) Z-F  cos  4or 

(257) 

W t 

cos  2cr^,  V. 

z 

\n 

/ t t 

V, 

-(I)  ^Ft 

sin  2ar  . V 

_ _ 1 

(1 

\ 3 

J EF  sin  4« 

(258) 

3 

\n/  t 

t 4 
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Vn 
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TABLE  33  FORMULAS  FOR  FLEXURAL  RIGIDITY  WITH  ISOTROPIC  SUBSTRUCTURE 


4 “i  ,-fcr) 

B 

. 3 22  \ h / W11 

ii 

h3  12  A * / I2 

— Di  . (Q°  -Q°  j 

3 66  2\  h 7 'Wll  U12' 


r ; n — r 


,3  16 

T ■ 

4 °26 

IT — I 1 


co*2(6-  6 j)  cos4(0-62) 


-00*2(6-6^  cos4(9-62) 


0 -cos4(6-62) 


0 -co84(0-62) 


sin2(e-61)j  -sin4(9-62) 


-■2sin2(0-61)  sin4(6-62) 


« 

VF  C082aj  + ( VF  sin2ff  ) 

^ ' 

— i i 

, \ 3 

r— 

ill. 

, 

[ 1 

F 

t Ftco84<»t)2  + ( 1 

pF^sin4o^) 

111  M 

2(D,  i' 

f D, 

2 

— i — L J — 
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TABLE  34  RIGIDITY  OF  CROSS- PLY  LAMINATES 


D,  , not  affected  by  stack- 


6 = w/4  (This  is  the  same  as  the  ±45  laminate) 


TABLE  35  TRANSFORMED  RIGIDITY  OF  CROSS-PLY  LAMINATES 


(n© stacking  sequence  effect) 


— 

b.  Stacking  Sequence  Effect  on  Symmetric  Cross-ply  Composites 


Assume  no  substructure;  i.e.,  zq  = 0. 

TABLE  43  STACKING  SEQUENCE  EFFECT  OF  CROSS-PLY  LAMINATES 


I 

1 


* m = any  positive  integer,  and  has  no  effect  on  stacking  sequence,  but  strong  effect 

on  D. . values  . 
ij 


TABLE  44  FLEXURAL  RIGIDITY  OF  CROSS- PLY  LAMINATES 


*2 

T 2 1 

Vl  + V3  R1 

R2 

— D' 
h3  U 

i 

1 

cos  26 

cos  46 

ii  D- 
3 22 

n 

■4  D'- 

i 

1 

-cos  26 

cos  46 

i 

-1 

0 

-cos  46 

•j 

n It 

h 

— D' 

.3  66 

0 

1 

0 

-cos  46 

i 

h 

12 

h3  16 

11  Di 

0 

0 

- q;sin26 

sin  46  | 

h 

Laminates 

5 

HEfl 

6i 

■ 

EfEfl 

62 

[0  /90  1 

1 m m J s 

0.75 

0 

1 

0 

[0  /90 

1 m mJ  2s 

.75 

0.375(  = — ) 

0 

1 

0 

[0  /90  ]. 

1 m m J 4s 

0.  1 875(  = 

0 

1 

0 

[0,/90./0-] 

0. 1875 

0 

1 

0 

[0/902/02/902/0]s 

0.0468 

0 

1 

0 

[0/90  }x 

o 

II 

qKj 

0 

1 

0 
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c.  Flexural  Rigidity  of  T- 300/5208  Cross- Ply  Plates  and  Beams 


I 

Ij  = 49.5  I2  = 26.9  Rj  = 85.8  R2  = 19.7  (GPa) 

h = .005x25.4i.l25  mm. 
o 

TABLE  45  RIGIDITY  OF  T- 300/5208  CROSS- PLY  LAMINATES 


d.  Numerical  Examples 


( 1)  Three  Point  Bend  Test  of  T - 30 0/5208  : [0,,/90,,]  , h = 8mm,  L = 10cm,  b = 2cm 

lb  lbs 


6848  186 


1357  0 (Nm)  d..  = 


(a)  Beam  along  0 on  outside 

(El),  = = .137  (kNm2) 

1 dll 


£ _ 48EI 


= 6576  kNm 


.146  -.020  0 

.740  0 


(kNm) 


(i)  Load  at  First- Ply  Failure  (FPF) 
e^,  = 3.9  mm/m 


Figure  66  Three-point  bend 
test. 


.\k  = — = 1.95  m'1  (268) 

z 

M = Djjkj  = 6848  x 1.95  = 13.4kN/unit  width 


10.72  kN 


.-.  p 

2Mb  2 x 

13.4  x .02 

4 32 j>!*s 


6 fpf  = iitr  = K63mm  (269) 

(ii)  Load  at  Ultimate,  assuming  X = 1503  MPa  or 


Figure  67  1-direction  of 
[0/90J . 


e^t  = 8.3  mm/m 


kUlt  = “P  = 2.08  m'1 


— = 2.08  m (270) 

6848  x 2.08  = 14.2kN  (271) 


Pmt  = 11.36  kN,  6UU  = 1.73mm.  (272) 


mm 


Figure  68  2-direction  of  [0/90] 


(i)  Load  at  First- Ply  Failure  (FPF) 


Load  at  Ultimate 


Figure  69  Load  deflection  curves  of  0 


beams  of 


e.  Natural  Frequencies  of  Transverse  Vibrations  of  Beams 

l 

n 1 

w = — \ 

" L2 


4 


dn^ 


(273) 


TABLE  46  FREE  VIBRATIONS  OF  BEAMS  WITH  UNIFORM  CROSS-SECTIONS 


End  Conditions 

s 

K2 

K3 

X for  large  n 

Clamped-  Free 

1.875 

4.694 

7.855 

(2n-  1)  /2 

Hinged- Hinged 

3.  142 

6.283 

9.425 

nJi 

Clamped- Hinged 

3.927 

7.068 

10.210 

(4n+  1)  /4 

Free-Hinged 

" 

" 

« 

" 

Clamped-  C lamped 

4.730 

7.853 

10.996 

(2n+  1)  / 2 

Free- Free 

II 

If 

" 

II 

(274) 

(275) 


Example:  T- 300/5208  beams  from  C®i6^®]6^s  P^ate 

L = 10cm,  b = 2cm,  h = ,8cm,  density  = 1600  kgm 
p = mass /unit  length  = 0.256  kg/m 


-3 


dn  = .146  (kNnf  \ ~ = 1 


37  Nm 


11 


- - x2 

n n 


\l  137 


2313X 


-1" 


256 


For  hinged-hinged,  u>j  = A x 2313  = 22832  s 


-1 


By  treating  3-point  bend  test  as  a 1-degree  of  freedom  approximation 


HIP' 

6576000 

* l 

"l  i 

6M  ~ 1 

| .025612 

22666  s 

The  difference  by  the  2 methods  is  negligible  as  expected. 


(27  6) 


The  effective  mass  M is  one  half  of  total  mass  of  beam  for  a beam  under  3-point 


bending. 

\ l 1 - l .1 

* 
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f. 


Special  Isotropic  Homogeneous  Plates  Consisting  oi  Orthotropic  Plies 


I 


' 


t 


TABLE  47  SPECIAL  ISOTROPIC  HOMOGENEOUS  PLATE 
[- 60 /O/iO^/O/- 60/60/0/- 60/0/60] # 

For  T-300/5208  material,  Ij  = 49.5,  I2  = 26.9,  Rj  = 85.8,  R2  = 19.7 


(GPa) 


r!  = V2  “ U4'  V3  = V4 


0-  6i  *«2  =o 


*1 

h 

12 

n 

1 

l 

h3 

ii 

12 

h3 

°12 

I 

- 

i 

12 

h3 

°16 

0 

0 

144  R1 


144  R2 


cos  29 


1 • 

— sin  26 


cos  46 


-cos  46 


- sin  46 


o 

1 5 

\c\ 

A.  G 

7 C 

QO 

12 

76 

75 

76.4 

? D 

i 

7 7*i 

7 ( 

0 

76. 

6 

3 

7 0 

0 

.9 

75.9 

4 D 

22.5 

22 

5 

22 

7 

22 

7 

22.6 

h3 

12 

ii  D 
. 3 ( 

26.8 

26 

8 

27 

0 

?7 

0 

26.9 

h 

12 

r3 

0 

03 

14 

•3DI 

38 

27 

0 

0 

Since  stacking  sequence  is  not  important  for  in-plane  modulus  A..,  this  special  laminate 


is  equivalent  to  [ 0 / 60  / - 60 ] ^ s 

, from  which 

£cos2<»t  = 

£sin2«t  = 

£cos  4e( 

t = Esin4-,t 

-J 

0 

(277) 

i 
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' 

Then  A /h 

- 

= A22 

/h 

= I 

1 + b 

1 1 

A12/h 

a 

''  !l  ' 

h 

, A 

l66/h 

~ 

= 


16 


= A26  = 0 


Then 

A. . is 
13 

isotropic 

12 

— — D..  for  this  plate,  the  laminate  is  also  homogi 

Since 

A. 

1 

1 

' h J 

(278) 
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• 2 >.  I 2?.  0 * ?•?. 

T ^ * 

Figure  71  Specific  rigidity  as  function  of  percentage  of  core.  Note  optimum  combina- 


• 2*a  0 


tions  of 

density  and 

stiffness 

ratios . 

, 
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4.  DETERMI NATI ON  OF  STRESSES 


a.  Basic  Relations 

TABLE  50  FLEXURE-CURVATURE  RELATIONS 


— 

kl 

kz 

k6 

M1 

M2 

M6 

M, 

D,, 

D,_ 

D,  , 

k 

d 

d , 

d , 

1 

11 

12 

16 

1 

1 1 

12 

16 

m2 

°21 

°22 

°26 

k2 

d21 

d22 

d26 

M6 

°61 

°62 

°66 

k6 

d6l 

d62 

d66 

Since  e,  = zk,  , e0  = zk-  , e,  = zk 

lie  CO 


Finally,  from  stress- strain  relation  of  the  t-th  layer. 


Figure  73  Flow  diagram  for  ply  stress  calculation. 


(281) 
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(4)  Stress  a. 


(t) 


„(t)  (t) 
= Q. . e 
»J  J 


(4) 


St  r e s s 


a<‘> 


n(t)  (t) 
Q. . e . 

U J 


TABLE  54  PLY  STRESSES  AT  VARIOUS  LOCATIONS 


z ( mm) 
t 


Stresses  at  z 

Stresses  at  z. 

t-  1 

t 

Figure  75  Stresses  and  Strain  across 

laminate  thickness. 
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d.  Stress  Concentration  of  Orthotropic  Plates  with  a Filled  Hole 


! 

i 

i 

M 


1 

I 


* 


. 


. 


Assume  a cross-ply  T-300/5208  laminate 

'°16/,016l 


(1)  Elastic  Constants 

From  previous  calculations 


6848 

“ 

186 

0 

D. . 

1357 

0 

(Nm) 

307^ 

Figure  76  Bending  of  plate  with  a hole 
filled  with  rigid  core. 


11 


160,  E 


22 


= 31.6. 


12 


.136,  G 2 = 7.2  (MPa) 


(291) 


Complex  parameters  are  roots  of  equations 
D22^  + 2(D12+ 2D66)Ki2+ D1 1 = 0 <292> 
M-  = - . 590  ± 2 . 1 67^  = 2.246  Qcos(±l . 30)  + i sin(±l . 30)^ 


= d 

t 1. 

.50 

Ocos(±  .65>+  i sin 

*■  d 

: 1. 

. 50 

(_.  80  ± . 60 ij  = 

± (1.2±  09 i)  (293) 

n = = ± 2.4i,  dt  1.8  (294) 


k = 

-PjP2  = l-2  + -9  = 2-25 

(295) 

(2)  Bending  Moments 

(a)  At  point  A 

^ +n 

M = (1  + — 

r 1 - v - 

c 

k 

•) 

= 1.86  1 

A 

12  21 

(296) 

j = ,050M  , 

[ = 0 

r0 

(b)  At  point  B 

M °2 

M = - — ( 

r 1 * wl2 v2 1 D1 

V1Z 

+ V 

l2n+k)  = 

- . 52  3M 

(297) 

= u \A  - - A7  17U 

M = 

n 

6 12  r 

’ re 
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(3)  Ply  Stresses 

From  previous  calculations 


(a)  At  Point  A 


M. 


1.86  ^ 

.SO  / M , then  k 


At  z = 4mm  (outer  most  ply) 


e.  = 4k. 

l l 


M , <r.  = Q^’e. 

1 ij  J 


At  z = 2mm  (Outermost  90°ply) 


e.  = 2k. 

i i 


M. 


<r.  = Q(.90'e. 
1 ‘J  J 


M 


FPF 


(b)  At  Point  B 


M 


Ult 


(298) 


M (299) 


M (300) 


M (301) 


(302) 


M.  J 

-.07, 

1 1 

l 0 J 

M,  then  k 


M 


(303) 


At  z = 4mm,  e.  = 4k 
i i 


At  z = 2mm,  e.  = 2k.  = < 


M.  . 


}m.  . a[ ’°>v( 


M (304) 


M (305) 


M 


FPF 


“Ult  " 


(306) 
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SECTION  VII 


PROPERTIES  OF  UNSYMMETRI CAL  LAMINATES 

1.  BACKGROUND 

An  unsymmetrical  laminate  is  one  that  does  not  have  mid- 
plane symmetry  in  the  layup  of  plies.  This  type  of  laminate 
has  not  been  used  extensively  in  actual  structures,  out  may 
be  considered  if  special  constraints  or  effects,  such  as 
minimum  gage,  aeroelastic  tailoring,  and  bimetallic  set- 
up are  desired.  Unsymmetric  laminates  will  warp  after  curing  and  cooling.  The  degree  of 
warpage  will  change  by  temperature  and  moisture  absorption,  in  addition  to  applied  loads. 

The  properties  of  unsymmetric  laminates  are  a little  more  complicated  than  the  symmetric  ones. 
Because  of  the  lack  of  symmetry,  the  laminate  will  bend  or  twist  when  an  in-plane  load  is  applied, 
or  it  will  stretch  when  a moment  is  applied.  This  coupled  response  is  unique  and  has  no  counter- 
part in  symmetric  laminates.  Unsymmetrical  structures,  however,  are  not  uncommon  in  real 
fife.  Floor  slabs,  fuselage,  and  many  other  built-up  structures  are  usually  unsymmetrical. 

Thus,  the  theory  of  unsymmetrical  laminates  will  be  discussed.  They  may  uniquely  fulfill 
requirements  not  possible  with  symmetric  laminates. 


Figure  77  Symmetric  versus  unsymmetric  laminates.  Symmetry  is  based  on 
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2.  COUPLING  MODULUS 

a.  Definition:  In  symmetric  laminates,  in-plane  and  flexural  behavior  are  independent  of  each 
other;  i.e.,  they  are  not  coupled.  But  for  unsyrametric  laminates,  they  are  coupled;  i.e., 
in-plane  extension  requires  both  in-plane  stresses  as  well  as  moments  (to  keep  the  laminate 
from  warping),  and  conversely, the  bending  of  a laminate  induces  both  in- plane  stresses  and 
moments.  The  coupling  modulus  can  be  defined  as  follows 


(307) 


h/2 

N.  , f ... 

L/z 

I 


Q..e.dz 
V)  J 


If  e.  = e.  + zk. 
i l l 


(308) 

(309) 


N.  = / Q.  ,(e°+ 2k.  )dz  = e.  / Q dz+ k.  / C 

1 J 1 J J 1 Jj  ij  J J 


e.  |U,,dz+  k,  /Q..zdz  (310) 


~ - 

A..N. 
lJ  J 

+ B..k. 
ij  J 

Alternatively, 

. h/2 

\A 

. J _ 

i 

cr . za  z - 
1 

■ W.  4 C . ZUZ 

J J 

(31 1) 

1 

1 

/* 

M.  = 

1 

J 

QU( 

O 

Bii 

,€  . 

1 J 

• = / Q..(e.  +zk.)dz  = e?  /Q..zdz  + k.  /Q..z2dz 

1 J U J 3 J J U J / lJ 


(312) 


Thus,  the  c< 

aupling  modulus  is  the  s 

ame  between  N and  k as  that  between  M and  e°.  Needless 

to  1 

say,  B. . i 
VJ 

s identically  zero  fo; 

r 8 

ymmetric  laminates. 

• For  [0/90]t,  -Bn  = B22,  B,2  = B66  = B^  = B26  = 0. 

* F°rf®/-e]f  B16'  B26*  °'  B11  = B22  = B12  ' ®66  = °* 
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b.  Formulas  for  Coupling  Modulus 


For  N plies  with  uniform  thickness  h 

r i 

(N  even) 


78  Integration  of  unsymmetric 
laminates  must  be  performed 
for  the  entire  thickness  from 
-h/2  toh/2.  (For  symmetric 
laminates,  integration  can  be 
limited  to  from  0 to  h/2). 
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TABLE  55  FORMULA  FOR  COUPLING  MODULUS 


*4  Ri 

_J_  R 
n2  R2 

h2  Bu 

^Ftcos2rr.t 

T,Ftcos4'lt 

7 B“ 

- ^f'tcos2'tt 

^Ftcos4rit 

2 

h2  B12 

-lFtc°S4,t 

4 B66 

n 

- 21^00840,. 

4 b,6 

h 

4 E^tsin2«t 

- £Ftsin4at 

_L  B 
,2  B26 
h 

4 HFtSinZnt 

2]Ftsin4rtt 

^ orientation  of  t- ply,  -n/2;t<n/2 


*3  B = 4 (B11+  B22+  2B12)  = 0 


*26  = 8 (BH+  B22-  2B12+  4B66>  = 0 


'IB 


= 1 ^-B11+B22)2+4(B16+B26)2 


rv  _ r"1  ■ 

= — \ (^Ftco»2rt)2+(X;FtBin2a)2  = Rj  ^ vJ+V 


R2B  = 8 ^ ^Bll4  B22* 

21 

316-  4B6( 

■/ 

+ 1 6 ( B. , 

- E 

,2 

= — \ (2^Ftco»4n 

)2h 

• (2ZFtsin4n)z 

= R2 

j 

v: 

< 

2+v2 
> 4 

N 
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0082(6-6^ 


cos4<e-&2) 


- cos2(6-  6 j) 

0 

0 

--IsinEfe-ej) 

- ^sinZte-Sj) 


cos4(9-  6.,) 


-co84(0-62) 


-cos4(6-  62) 
- sin4(9-  &2) 
sin4(9-62) 
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TABLE  58  CALCULATION  OF  REDUCTION  FACTORS  (V's)  FOR  COUPLING  .vIODULUS 
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TABLE  59  CALCULATION  OF  REDUCTION  FACTORS  (V's)  FOR  COUPLING  MODULUS 


I ! I 


T IT 


TABLE  62  CALCULATION  OF  REDUCTION  FACTORS  (V's)  FOR  COUPLING  MODULUS 
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(3)  [0./90  /0./90  .]  . B..  for  this  laminate  is  exactly  one  half  that  of  [0./90  ]_. 

4 4 4 4 1 lj  8 8 T 

So  for  N - even,  B../N  = constant. 

11 


3.  COUPLED  BENDI NG  AND  EXTENS I ON  OF  LAMI  NATES 


a.  Constitutive  Relations 


Between  2 generalized  stresses  (N.  & M.),  and  2 generalized  deformations  <e°&  k ).  there 

11  i i 

are  6 possible  relations,  all  of  which  are  listed  below  in  terms  of  the  original  modulus 


matrices:  Inplane  A , Coupling  B.,  and  Flexural  D 
l3  i)  ij 


0 • c 


D 


} 


a 

b = 
d = 


= A 


-1 


.-1 


(327) 


-aB 


. Ba  D-BaB 


B-AbD  Ab 


-bD  b 


r „ *\ 

f 

N 1 

A-BdB 

Bd 

i 

- 

.k  J 

_ - dB 

d . 

Db 


-bA 


B-DbA 


(328) 


(329) 


(A-BdB)' 


-(D-BaB  )_1Ba 


-(A-BdB)  Bd 
(D-BaB)'1 


(33d 


All  submatrices  are  3x3,  which  are  numerically  simple  to  calculate.  The  unpartitioned, 
original  matrices  are  6x6,  which  normally  require  bigger  computer  to  evaluate.  A sample 
problem  will  be  solved  in  the  next  section  using  the  3 x 3 submatrices  only. 


b.  Numerical  Example 
T- 300/5208 


QU  = 


182 


2.9 

10.3 


/ f*ir\ 


r °g /9°gJ  T»  h = 2mm 


0 

0 

7.2 


GPa  (331) 


-/ 


wet 


'8 


j 1=0 


Figure  80  2-layer  unsymmetric 
laminate. 
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192.3 


5.8 

192.3 


0 

0 

14.4 


MNm'1  (332) 


J 


’-85.85 

0 

0 

64.  1 

1.93 

0 

85.85 

0 

kN  (333)  D..  = 

»J 

64.  1 

0 

Nm 

0 

0 

0_ 

. 

4.8_ 

(334) 

| A | = 532017 


(GNm' V1 


(335) 


|b|  = 0,  B..  is  singular,  therefore,  there  is  no  inverse. 


1 

d| 

= 

19740 

di 

j 

= 

15.58 

15 

.469  0 

.58  0 

208 

(kNm)'1 

(336) 

r.44^. 

I 1 

4ft 

0 
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±3a 

1 3. 

0 

4o 

0 

0 

0 
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(337) 

*- 

J 

146 

-15 

1.48 

0 

aB 
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i: 

3.48 
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0 
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_1) 

-1 
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-1337 
40.26 
. 0 


-40.26  0 

1337  0 
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- 1 


(342) 


BdB  = 


114 


3.46 

114 


MNm 


( 343) 


A- BdB  = 


78.3 


2.  34 
78.3 


0 

0 

14.4 


MNm 


- 1 


(344) 


(D-BaB)_1  = 


38.8 


-1.16 

38.8 


0 

0 

208  J 


(kNm) 


-1 


(345) 


(A- BdB) 


-1 


12.8 


( D-  BaB)"  * Ba  - 


-.382  0 

12.8  0 

69.4 


(CNm'V1 


*1 

0 

0 

17.29 

0 

(MN) 

0 

0_ 

(346) 


(347) 


(A-  BdB) ' 1 Bd  = 


17.13  0 0 

0 17.13  0 

0 0 0 


(MN) 


-1 


(348) 


151 


c.  Summary  of  Constitutive  Equations  for  T-300/5208  [0  /90  ] , h = 

8 8 T 


Units:  - 


N is  MNm 
M is  N 

e°  is  lO'^ro/rn  or  mm/m 


-85.85 

0 

. 0 


0 

i 

0 I 
14.4  | 


-85.85 


0 

85.85 

0 


| 64.  1 


0 

85.85 

0 


1.93 
64. 1 
0 


Units: 


12.8 

1 

u* 

00 

rsi 

0 

17 

- .382 

12.8 

0 

i 

• n 

o 

f) 

Ao  a 

o 

/ , -j 

. J 

17 

0 

0 

38.8 

0 

-17 

0 

-1.16 

. o 

0 

0 

0 

-1.16 

38.8 

0 


Units  ; 


(GNm'V1  I (MN)"1 


(MN)-1  J (kNm)~ 1 


(1)  Determine  under  uniaxial  tension. 

1(  FPF) 


= 12.8  Nj  . 


1 7Nj  , 


•••  ej  = (12.8+  1 x 17)Nj  = 29.8  Nj 

z=  1 mm 
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Let  e 


FPF 


3.91  mm/m 


(354) 


N1(FPF)  29. 8 = -131  MNm 


ffl(  FPF)  = Nl(FPF)/h  = 65*6MPa 


(2)  Determine  if  the  same  plies  were  symmetrical; 

e.g.,  instead  of  [0^90^,  use  [ C>4/904]  g . 

For  the  latter  laminate. 


A remains  same,  but  B.  =0. 
‘J  ij 


ej  = 12. 8Nj  , kj  = 0. 


N1(FPF)  = 12.8  ‘3°5  MNm 


“’ll  FPF)  = 153  MPa 


(About  I 

1.3 

times 

hight 

:r  than  the  unsymmetric  laminate) 

(355) 

(356) 


(357) 

(358) 

(359) 
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d.  Alternative  Constitutive  Equations  of  T-300/520^  [0  /90-  J Laminate  ( From  Equation  328) 

o o J T 
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0 0 14.4  | 0 


1337 
-40.26 
■ 0 
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(1)  Determine  ply  stresses  in  uniaxial  extension  (N^  only)  of  a tubular  specimen 


of  [0/90]  ply  orientation. 

O O 

For  long  cylindrical  tubes;  k = k =k,  = 0 

1 Zb 

(364) 

o 

ei  = 

5. 20Nj 

(365) 

o 

e2  = 

-.  157Nj 

(366) 

M!  = 

44  6N 

(367) 

m2  = 

1 3. 48Nj 

(368) 

Moments  are  induced  because  curvature  is  prevented. 


At  FPF  - 


N 

l(FPF) 


3.91 

5.20 


.7  52  MNm 


(369) 


° 1(  FPF) 


375  MPa 


(370) 


The  tube  has  higher  FPF- stress  than  both  unsymmetric  (65.6)  and  symmetric  (153) 
laminates.  The  induced  moments  have  beneficial  effects. 


(2)  Are  they  still  beneficial  if  the  tube  is  under  axial  compression? 


i 


r " T ** 

(3) 

What  is  the  difference  between  shear  response  of  flat  laminates  and  tubes? 

1 

| > 1 l 

.. 
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(1)  Sample  Problem 


i 

Maximum  water  absorption  in  a typical  epoxy  is  6f . 

What  is  the  maximum  amount  of  water  which  can  be  absorbed  in  a Gr/Ep  (v^.  = 0.65i? 
Specific  gravities  for  the  epoxy  and  composite  are  1.25  and  1.6,  respectively,  and  the 
void  content  is  0.4ff'. 

Solution 

Since  c,  = 0,  c = (c  v s +v)/s=  1 .87% 
f m m m v 

6%  epoxy 

= .65 

S = 1.25 
m 


S{  = 1.6 

V =0.4% 


SECTION  VIII 
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STRESS  AND  DEFORMATION  DUE  TO  CURING  AND  SWELLING 
1.  UNI  DIRECTIONAL  LAMINAE 

a.  Moisture  Concentration 

Dry  State 

M - M + M, 
m f 

V = V +V,.+  V 
m f v 


Wet  State 


M'  = M+M  4-  M,  + M 

mv/  fw  vw 


Moisture  concentration 


M'-M 

M 


= c m + c,m,+  M /M 
mm  f f vw 


,cmvmPm+c(vfPf+v„Pj  / P 

m m m iit  vw 


(c  V S + C V.S  + v ) / s 

m m m f f f v 


V : volume  of  voids 


v 


M : 

mas 

s of  water  in  matrix 

mw 

M.  : 

fw 

ii 

" " fiber 

M 

n 

" " voids 

cr , . p Mr , , / M.  r , 

[f,m]  [f.mjw  [f.rn] 

s ; specific  gravity 


In  many  cas 
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« 0 

c = 

( c V 

m 

s + v ) / s 

ID  id  V 

’ " | 

11  -I  - i- 

L 1 1 LI 

i — r- 

(371) 

(372) 


( 373) 


(374) 


(375) 


(376) 
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b.  Approximations  for  Curing  and  Swelling  Strains 


v E e + v E a 
f f f m m m 

r| . v E + v E 
Iff  mm 


: longitudinal  strain 


e„  = ve+v  e +v.w,e,+  v v e 
ff  mm  fff  mmm 


transverse  strain 


N = T (curing),  or  H (swelling),  or  T + H (both) 


Cff,  m]  6[f,  m]  + eff  m]  ! 8tra*ns  resulting  from  a change  of  temperature  and 

moisture  concentration,  measured  from  the  initial, 
e stress-free  state  to  a final  state. 


e 

y Final 

/ 

Initial  Q 

Elastic  constants  at  the  final  state  of  interest. 


Curing  strains 


V-m]  " “[f.m]  AT 

AT  : (temperature  of  interest)  - (curing  temperature) 


strains 

0 
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0 
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re 
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(383) 


XH 

m 


H 

CL 


H 


1 

3 Sm 


« 0 


1+  ^ 
m 


1 + v 


volume  additivity 


v s c 
m m m 


->  o 

In  general 

tjj 

u 

e = 

T 

a (c  - c ) H (c  ) 

1 O O 

e T 

H 

ftT  : 

Swelling  coefficient 

(384) 


(385) 


(386) 


H ( ) : Heaviside  step  function 

: = v^/s  '•  *n  general,  threshhold  value  of  c below  which  no  appreciable 

swelling  occurs. 


(1)  Sample  Problem 

What  is  the  expected  swelling  coefficient  of  the  Gr/Ep  in  sample  problem 
l.a.(l)?  Poisson's  ratio  of  a typical  epoxy  is  0.35. 


c = v /s  0.25% 

o V 

XJ  1 + V 

i 

H m _ 

Ot  = 3 s = 0.72 

— 1 — 

f - 

1 

+ — j f-  -) ‘ * - 

- * _ * , 1 1 i -1  - 1 ' |-  -J 

. 

! I T 

1 

1 1 i — 

J LJ J J.  . 

JlJZ 

jj 

1 11  1 i 1 1_L 
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c.  Approximations  for  Residual  Stresses  in  Unidirectional  Laminae 
Refer  to  Section  IV. 


Matrix 


’VfW^m 


m — R 

<r 

v,  mL 


v E 
m m 


E^l*m‘V 


r\ , v E • f v E 
'Iff  mm 


Residual  stresses  manifest  themselves  in  the  residual  fringes  in  photoelastic  composites. 

d.  Change  of  Densities  Due  to  Curing 
Volume  changes 


m ‘If  m 


, "TfVm 


' (e  - e ,/n , 

mm  i ] 


v E 
m m 


Av 

f 

v. 
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*n 
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L f ’VfV'^m 


In-situ  densities  (dry)  after  curing 


P = p 
m mo 
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m 
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For  Gr/Ep,  v = 0.35, 
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T T 

e = 0.  *e x as  0.779. 
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(394) 


1 + 3e. 


P,  = P 


f°  1+  ZeJ  (1  + k£)+(1-  2i/f)e^ 


fo 


p t p : densities  of  constituents  themselves  at  the  temperature  of  interest 
mo  to 


Composite  density 

M = 

M + M, 
m £ 

V = 

M /p  + Mf/Pf+  v 
mm  t t 

V 

P = 

P^  V + P f vf 
mm  t t 

e. 

Change  of  Densities 

Due  to  Moisture  Absorption 
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= 
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f.  Remarks  on  Thermal  Expansion  Coefficients  of  Composites 
Assume 


c[f.m] 

T 

«ff.m] 


a..  ,(T-  T i 

[f,  m]  o 


independent  of  T 

initial  reference  temperature 


(403) 


From  Eq.  (377) 
T 

rtT  (T)  = 

1 J 


^VfEf(T'+VmEm(T' 


(404) 


wnere 

T 

aL(Ti  = 

T 

eL(T) 

T - T 

o 

Comments 


Define 


ttL(T2’Tl' 


T T 

eL(T2)-  eL(Tl) 

T2-Tl 


(405) 


1.  T 

c 

is 

I 

not  arbitrary  and  e must  be  measured  from  T . 

9 T . 

2 . n.  _ c 

T T 

lepends  on  T while  a and  ft 

i m 

do  not. 

i 

* 

(406) 


r r 

urn  n,q.  \ 

3 1 ( 

) 

n 

n i 

1 rn 

T 

rEf<T2>"f  + V 

E 

m 

(TZ)  a. 

T 

m 

J 

VfEf(T2)o.f 

f V 

] 

E (T,)flr 
m m c m 

"l 

(T 

2’1 

rl' 

= 

+ 

V 

VEiT; 

,)+  V 

7 rr 

E (T_) 

T|  lVf^/  ^2' 

+ V 

E ( 
m m 

T2' 

b 11 

\ 

T 

1 

’\ 

t rTl'*f 

+ V t 

m 

m 

Tl>"rr 

l] 

T1 

- 1 

0 

11 1 

VfEf(1 

V' 

■f  v E 
m 

(■ 

m 

"l' 

7 

T, 

i 

ri 

(40  7) 

r 

T 

*T 

(T 

a.I 

r.) 

becomes 

independent  of 

1 temperatur 

e ii 

r sc 

> are  the 

moduli . 

c 

1 

. 

- 4 
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2.  multidirectional  laminates 

jj!  u M Hrn  ! j | 

a.  Curing  and  Swelling  Strains 

Constitutive  relations  of  a constituent  ply 


<r.  = C (e. 

1 lj 

j' 

+ CIA 

.(eA 

V 


in- plane  stresses 


A ^Aj^j  ej^  + Gab^b'  V : ou*>of-plane  stresses 
Reduced  stiffnesses  Q. . 


(408) 


(409) 


Q..  = C 

lJ 

(410) 

o'.  = Q..(e.-  e.)  + C C_1 

1 lJ  J J iA  AB 

"b 

(411) 

Ir 

>- plane  strains  and  curvatures 

o 

e;  = e.  + zk. 

1 J } 

(412) 

Classical  laminated  plate  theory 

— - 1 . 

\—r  t — j 

O , 
Or.  «Jt  k 

^ independent  of  z 

/a.dz  , 

/V.  zdz  = 

N = A e°  + B..k.  . 
1 *•)  J 'j  j 

M = B.  ,e°  + D.,k  - 
1 }i  } ij  j 

■ NN  + T <r 
i iB  B 

MN  + R,  <r 
i iB  B 

(413) 

(414) 

h/eAdz 

where 

€A  " CAB  11 B ‘ 

CABCBj  *j  ‘ 

c;1*  z 
AB  Bj 

kj  + GABeBjej  + eA  <415) 

r.rN  . -N, 

/•h/2 

1 

in.  , m. 

' i i J 

J Q.,e  [1,  zjdz 

-ti/2  1 ■> 

1*10) 

j 

fTiB-  RiBJ 

rh,z  -i 

- i,z  Vab 

[ 1 , zjdz 

(417) 

-■i 

rhfZ 

J ^ >dz 

~ j — | — j — j — 4— 

-h/2 

(418) 

— J_J LL 

Ml 

1 11  1 [ j j 

i 1 1 1 1 

-j"*"  ii 

ITT 
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1 


Nonmechanical 

laminate 

s strains  a 

N. 

l 

ii 

H* 

= 

O' 

A 

= 0 

0 -IN  -IN 

e.  = F..  (N  - B.,D  M ) 

1 lj  j jk  Kn  n 


kN 

X 

II 

o~}  (MN  - : 
11  1 

Bikek> 

e . = + C"’  Cn.  e.  - C'*  C„.  e°  - C**  C z kN 

A A AB  Bj  j AB  Bj  j AB  Bj  j 


(419) 

(420) 

(421) 

(422) 


where 

F 

A 

B. 

i 

I 

1 R 

ij 

ij 

Ik 

kn  nj 

ot  rain-  displacement 

relations 

o 

= 

5 u/dx 

o 

- 

Sv/r'x-, 

e 

o 

✓ = C 

iu/Jx,  + 5v/ 

dx. 

1 

(424) 

kl 

= 

- 

<^w/c 

2 

,X1‘ 

k2 

= 

l ; 

"w/cW. 

l 

l * 

k6  = 

2 

2d  w/dx  dx, 
1 < 

> 

(425) 

All  the  elastic 

constant! 

» are  taken  at  the 

final  state  ( 

T.I 

4). 

In  the  material 

. symmetry  axe 

IS  c 

>f  transversely  isotropic  laminae 

el 

/U 

1 

E 

't  1 

( 

LT 

/U 

l 

0 

0 

= 

E 

: /u. 

0 

(426) 

u 

T 

1 

- 

Gi 

-T  . 

J 

1/1 

Z 

0 

0 

T 

l 

— 

! 

IB 

l l + 

*T' 

r)/ 

Et 

: { ) 

1 

- 

* 

1/gl: 

r - 

h 1 



I 

'i 

= 

Q. 

> 

ei> 

1 

1 

1 

1 

. 1 

J 

— 

1 

_J 

1 i 

i 

.11 

1 

J 

1 

_l 

j 

n 

1 
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0 


T 

c = c 

iA  Ai 


1LT(1  + ‘tt'ET/,U2  0 


( u + v i/  ) E /U_  0 

' TT  LT  TL'  T 2 


(428) 


U1  = 1_  ‘lt'tL 


U2  ~ ^1+  VTT^1*  ^TT*  2*LT*TL^ 


(429) 

(430) 


b.  Curing  Strains  and  Residual  Stresses  of  Symmetric  Laminates 


N 

B..  = 0 , M.  = 0 , 

i)  1 


, N 

k.  = 0 

t 


(431) 


TABLE  63  LAMINA  CURING  STRAINS  FROM  CURING  TEMPERATURE  [12] 


Vf 

8 

T 

CL 

% 

T 

eT 

% 

T 

o 

°K 

Boron/Epoxy  (B/Ep) 

0.50 

2.03 

0.118 

0.451 

450 

Boron/Polyimide 
( B/PI) 

0.49 

2.00 

C.081 

0.443 

450 

Graphite/Epoxy 

(Gr/Ep) 

0.45 

1.54 

-0.009 

0.  318 

m 

Graphite /Polyimide 
(Gr /PI) 

0.45 

1.54 

0. 

0.391 

450 

Glass /Epoxy 

(Gl/Ep) 

0.72 

2.13 

0.072 

0.530 

435 

j 

j 

.J 

1 

l _J 

1 1 1 

, . , 

■ 4.  1 — J 

. 

1 1 

( 

r — r 

. . 

1 - 

LLlj 
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TABLE  64  LAMINA  MECHANICAL  PROPERTIES  [12] 

el 

et 

glt 

1/ 

X 

Y 

S 

GNm  2 

GNm  2 

GNm'2 

LT 

MNm  2 

MNm  2 

MNm  2 

B/Ep 

201 

21.7 

m 

0.17 

1375 

56.0 

62.3 

B/PI 

222 

14.5 

7.7 

0.16 

1040 

10.8 

25.9 

Ci/  Ep 

190 

7.10 

n 

0.10 

1115 

41.9 

61. 5 

m 

216 

4.97 

n 

0.25 

841 

14.9 

21.7 

Gl/Ep 

60.7 

24.8 

12.0 

0.23 

807 

46.0 

45.0 

t 

. 

TABLE  65  CURING  STRAINS  OF  [Q2/±45]  LAMINATES  AND  CORRESPONDING 
STRESSES  WITHIN  45°  PLY  [12, 19] 


■B9 

T 

'*90 

. % 

•i 

m 

EXP. 

CAL. 

Y 

s 1 

B/Ep 

0.123 

0.109 

0.234 

0.259 

0.171 

0.184 

0.99 

0.13 

B/PI 

0. 100 

0.075 

0.214 

0.186 

0.163 

0.131 

4.09 

0.  33 

m 

-0.018 

(-0.009) 

-0.019 

0.104 

(0.068) 

0.122 

0.036 

(0.026) 

0.052 

1.09 

0.11 

m 

0.002 

-0.004 

0.031 

0.051 

0.017 

0.028 

1.20 

0.11 

Gl/Ep 

0. 105 

0.117 

0.260 

0.361 

0.182 

0.239 

1.39 

0.65 

Remarks:  Numbers  inside  the  parentheses  for  Gr/Ep  are  measured  during  the  cooling  stage  of 
curing  and  the  heating  stage  of  postcuring. 


c.  Warping  of  Asymmetric  Laminates  Alter  Curing 
Input  mechanical  properties  at  room  temu.-rature 


Material 


Gr/Ep 

Gl/Ep 


et 

glt 

lLT 

GPa 

11.2 

GPa 

7.10 

0.33 

13.1 

4.70 

0. 30 

Input  curing  strains  for  lamina 


Material 


e , mm/m 
L 


eT  , mm/m 


e , mm/m 
L 


e^,  , mm/m 


Gr/Ep 

Gl/Ep 


0. 304 

1.  345 


- 3.503 

- 4.077 


0.267 

0.743 


- 2.167 

- 2.235 


Amount  of  \ 

varping 

1 

i, 

kT 

2 T 2 T 

W 

‘ 2 ( 

kl 

x + k y + k.  xy  ) + c x + c.y+  c 
Co  12 

.1 

T 

w 

= 

0 at 

(0,0)  , (a,  0)  , (0,  b) 

For  f + 0; 

/-  e 

l]  laminates 

T kj 

= 

kT  - 
*2 

: 0 

T 6 

, w (a, b)  = - ab 

max  2 

2 

i 

ci 

<t  - 

2 kl 

a 

- o Cl  - 

2 

1 T 

2 

"2 

b - 

— k 1 
2 2 

b 

° C2  = 2 

DEFLECTION 


ANGLE,  DEG 


Figure  81  Maximum  deflection  of [9,/- 9 ] 

r ~ *•  4 4.  J 


laminates  after  curing,  a = b = 6.35cm 


Swelling  Strains  of  Symmetric  Laminate 


Absorption 


Desorption 


Properties  of  Gr/Ep  (AS/3501)  lamina 


EL'  GPa  ET’  GPa  GLT'GPa  "LT  ‘'tT  eL’  mm^rn  eP,mm/m 

153  11.2  7.10  0.33  0.33  0.127  -4.228 


0.  2 0.4  0.6  0.8  1.0  1.2  1.4  1.6 

Moisture  Content,  % 


Figure  82  Transverse  swelling  strain  in  [0]  Gr/Ep  (AS/3501)  subjected  to  absorption  and 
desorption  tests.  Note  that 


aT  = 0.5' 

c 

'o 

= 0.47o. 

0 0.4  0.8  1.2  1.6 

Moisture  Content,  % 


■ - • f-4  j j — j- — | — | — I — }-  j"  ■ f I { f ' 4 — 

Figure  83  Swelling  strain  in  thickness  direction  pf  Gr/Ep  laminates. 


e.  Warping 

of  Asymmetric  Laminates 

Due  to  Swelling 

_ , N 

Total  w 

T H 

= w + w 

For  square  [0/90]  laminate  of  side 

. N 

a,  k = C 

1 . k1 

M , N 

, = k?  • 

O 

. 2 _ 

N 1 / a \ T H 

w = -7  7'  <k,  + k,  ) 

N 

w = 

0 at  (0,  0) 

inaA  t-  ^ s-  j i l 

} 


Figure  84  Maximum  deflection  of  [0^/90^]^  Gr/Ep  (AS/3501)  laminate  due  to  curing 
and  moisture  absorption,  a = 7.58cm. 


(441) 


(442) 
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SECTION  IX 

r pi — 

FAILURE  THEORIES 

1.  INTRODUCTION 

Definition  of  failure: 

Any  consistently  definable  point  on 
the  stress-strain  relation 

Simplifying  Hypothesis: 


Po,»r 

Only  c 

>ne 

point  of  : 

failure  associated 

with  one 

state  of  stress 

Stress 


R«yt*n! 


On  Stft  ef 

nonl.neu.nty 

Figure  85  Definition  of  failure. 


Geometric  interpretation  of  failure  (in  stress  space) 


Dimension  of 
Physical  Space 


Dimension  of 
Stress  Space 


Failure  Contour 


One  (Xj) 


One  (ffj) 


4— 


Points 


-1  V ° X1  *1 


Two  (x  , x2  ) 


Three  a^) 


3-D  Surface 


H- 


Three  (xj.x^Xj) 


six 


6-D  hypersurface 


Failure  Criterion: 
Representation: 
Application: 


Definition  of  the  loci  (or  contour)  of  failure  surface. 

Material  is  safe  for  any  state  of  stress  within  the  failure  surface. 
Predictive  design. 


Strength  Anisotropy  - Failure  surface  dependent  on  material  orientation 


1 -D  example 

; • ; • • 


: i -) I I I L_.  - 1 1-  i i 


1 J 


X,'  X, 

8 A 


T*n»;un  strenftli 
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2.  GEOMETRIC  INTERPRETATION  OF  EMPLOYING  FAILURE  CRITERION  TO  INTERROGATE 
POTENTIAL  FAILURE 

A given  state  of  stress  Oj>  O 2-D)  can  be  represented 


same  direction  as  the  stress  vector.  The  strength  in  this  direction 
can  be  characterized  by  a strength  vector 

171  ■ I,,*2  * „2*2  * „6*V/2  (-45, 

* * * 

where  (j  j ’ <*2  ' P°*n*  on  ^e  failure  surface. 


Figure  86  Failure 

surface. 


Safety  factor  can  be  considered  as 

= 1 imminent  failure 


1 


-jL  = _ 

^ Safety  factor 


<1  no  failure 


(446) 


Utility  of  Failure  Criterion 

A,  For  a given  state  of  stress  Oj.  interrogate  whether  failure  is  imminent 

B,  For  a given  state  of  strain  e^,  interrogate  whether  failure  is  imminent 

C,  For  a given  stress  ratio -Si- , 1 - etc,;what  are  the  failure  stresses? 

D,  For  a given  state  of  stress,  what  is  the  margin  of  safety? 


Applications  for  different  failure  criteria  will  be  illustrated  by  examples.  All  examples 


will  use 

! the  following  mechanical  i 

properties  of  a 

0" 

graphite/epoxy  lamina 

17. 

6 

0.  581 

0 

- 

1 J 

.3  . 

1 

\i 

Q 

ij 

— 

0. 

58 

1 

1.  21 

0 

> 

t 1( 

> 1 

tsi, 

fin /in 

(447 

0 

_i  0 

0.760 

[1 

J 

1 < 

1 

i 

1 

. 

1 
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t 


r 

0.057 

-0.  027 

0 

= 

-0.027 

0.  838 

0 

x 10  in/in/ksi 

0 0 1.32 


(448) 


X1  = 

« 149  ksi 

X2  = 6.3  ksi 

X6 

■ 10.  5 ksi 

x'  = 

Xi 

= 103  ksi 

X'  = 18.2  ksi 

X' 

, 

2 

6 

= 10.  5 ksi 

(<7l  = 

177  ksi 

Biaxial  Strength 

f * 

\ 

2 “ 

- 12.  8 ksi 

t 


(449) 


(450) 




L— 

' 

— 

— . 

i 

i 

"T 

■"  J 

J 

id 
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3.  FAILURE  CRITERIA  FOR  ONE  GIVEN  ORIENTATION 


NON  - Interacting  Failure  Criteria:  biaxial 


not  change  uniaxial  strength 


Maximum  Stress  Criterion 


Maximum  stress  criter 


ion  in  stress  space 


Maximum 


in  Strain  Space 


For  linearity- 


maximum  stress 


criterion  in 


strain  space 


xl“ 

C11 

xi 

cu_ 

X2 

c,- 

22 

X' 

_2 

r* 

The 

strain  allowables 

are 

e. 

* 

12 

- e„ 

°16 

. 

i 

C 

11 

L 

\ 

cn 

C 

12 

C16 

1. 

€ . 

i = 

' C 

11 

e2 

C11 

C 

'13 

C26 

e< 

i 

i = 

* b 

' c.„ 

61 

- 

c„ 

r 

C.C 

r 

e. 

12 

e. 

* 

26 

6 

V. 

1 

^22 

A>  L 

_ °16 

C26 

1 1 



) 

C 

66 

*1 

C66 



c 

16 

C26 

1 

e. 

3 

'c 

el 

C66  ' 

; 

□ 

DO 

1 

DO 

(2)  Example 

Given  stress  <7.  interrogate  failure  by  maximum  stress  criterion 
CTj  = 100  ksi 
= 8 ksi 

cr^~  9 ksi 

Substitute  into  Eq.  ( 451  ),  if  any  one  equation  not  satisfied,  failure  occurs 


100  <X1  = 

149 

no  failure 

8>X2  = 

6.  3 

failure 

9<X6  = 

10.  3 

no  failure 

(3)  Example 


Given  strain  interrogate  failure 
_ „ , ..  ,„-3 

3 


1 

e2  1 

= 6 

: 

s 2 > 

x 10 


x 10 


-3 


Failure  for  lamina 


i __J 

~ - | 

. 

— 

J 

i 

□ 

. 
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(4)  Ex-imple 

For  the  following  state  of  stress  (ratios),  using  maximum  stress  criterion,  find 


the  failure  stresses. 

CTj  = 70 

*2  3 3 

a,  3 5.5 


We  need  to  find,  along  the  stress  ratios 
what  is  the  intercept  of  the  failure  surface. 


°_Z  °_6 


■Sftrwies 


Since  the  maximum  stress  failure  criterion  (a  rectangular  parallelpiped)  is 
a surface  bounded  6 sides,  the  stress  vector  J may  penetrate  any  one  of 
the  sides.  Accounting  for  the  signs,  we  may  interrogate  the  positive  surfaces 
only. 
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(5)  Example 

For  the  state  of  stress  given  in  example  4,  what  is  the  stress  vector, 

strength  vector  and  the  factor  of  safety  predicted  by  the  maximum  stress  criterion. 


b.  Maximum  Strain  Criterion 
E!=S11X1^  °1 


El’ 

= Sllxl> 

- 61 

E2 

= S22X2- 

°2 

E J 

= S„.,X'> 

- 

2 

22  2 - 

2 

E6 

= S66X6  ^ 

e6 

(453) 


E6,=  W**  6 6 

E^  Positive  Uniaxial  Maximum  Strain  in  i component 


E.1  Negative  Uniaxial  Maximum  Strain  in 

i component 

Geometry:  Rectangular  Parallelepiped  in  Strain  Space 


~ 2 

— 

E; 

l ~ 

e, 

P 

/ 

L 

‘ 

P 

;2  ~ 

i 

| 

1 

' 

1 

t 

■ 

l 

“1 

1 

P . 

rl 

E'  _ 

Figure  8? 

S Maximum  strain  criterion  in  s 

train  space. 

1 i — 1 

i i _L 

1 
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Maximum  Strain  Criterion  in  Stress  Space 
For  linearity  e.  = S.^cr,  combine  with  Eq.  (453). 
The  allowable  stresses  are: 


12 

' 1 ' ’ sn  1,2 
S12 


— ^CT.  + X, 
S11  6 1 

S16  . 


“ii 

S12 

= - - — 

(7, 

S 

11 

^26 

<7, 

+ 

X. 

c 

*2 

s 

22 

S12 

S22 

1 

al 

- 

s 

s 

s 

’22 

26 

22 

6 

*6 

- 

X2' 

d 

>* 

maximum  strain 

criterion  in 

stress  space 

(454) 

S,A 

s 

26 

° b = 

c 

hb 

- 

s 

— i 

L L 

*2 

+ 

X6 

e 

< 

DO 

5 . . 

UD 

s_ . 

°b 

= 

16 

c 

al 

- 

c 

26 
— I 

or2 

- 

y ' 

f 

1. 

'66 

o 

66 

(1)  Example 

For  graphite  epoxy  lamia, 

find  the 

maximum  s 

itrain 

failure  criterion  using 

me  given  data 

# 

1 

' 



i 

1 

L J 1 

, 

[ 

' 

J 

. 

J 

i i 

1 

, 

1 



l 

— 

1 

1 

J 

. 

, 

J 

1 _ 

1 

| 

_J 

J 

| 

1 

. 

I 

’ 

1 
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(2)  Example 

Giver,  stress  0.,  interrogate  failure  by  maximum  strain  criterion. 
0^  = 100  ksi 

^2  = 8 ksi 

0^  = 9 ksi 


(3)  Example 


Given  strain 

e ^interrogate  failure  by 

r maximum  « 

strain  criterioi  . 

e , 

= - 

•5. 

5 x 10'3 

1 

_ 1 

e2 

= 

6 x 10 

*>.  . = 

2 10'3 

6 

. . 

- - - l _ 

> 

* 

— | * 

— — -L  — .i. — 

_LJ 

L _ 1 - _ 1 1 

: 

— H 
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4.  lit 

tfTERACTING  FAILURE  CRITERION 

a. 

j 

formulation 

f(0rj 

= F 

’iffi 

+ 

F.  .o.o.  + F 
iJ  i J ij 

ikai 

0 0 

■ J 

k + 

Important  features 

(1)  Each  term  Scalar  product;  independent  of 
material  coordinate. 

(2)  Can  be  expanded  to  include  sufficient  (but  not  excessive)  number 
of  terms  to  describe  a given  composite 

For  a 2nd  order  polynominal  i.e.  , 

F.O.  + F..O.O.  = 1 expanded  for  two  dimension,  i.e..  i,  j = 1 , 2,  6. 

l l ij  l j 

The  failure  tensors  F.,  F.j  are  determined  by  9 experiments 


■ fcr > = 1 


(456) 


I 


“•-liWr  *TlAaj  X1 


Reference  x. 

(i) 

X. 


♦ 

x: 

(2) 


(3  r (4) 


K — 
(5) 


L 

(6)  (7J 


\\u 


T"  01  ' 

a, 

18)  (91 


Figure  89  Reference  coordinate  definition  and  9 guiding  experiments  for  measuring 
failure  tensors  in  a quadratic  polynominal  failure  criterion. 


For  tests  1 through  6 
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For  tests  7 through  9 (or  any  multiaxial  experiments): 
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The  optimal  biaxial  ratios  can  be  solved  from 


-(F.B  + F.)  ± [(F.B  + F.)2  + 4(F..B2 
1 1 1 J 11 

2(F..B2  + 2F..B  + F..) 

>J  JJ 


+ 2F.  .B  + F..)l1/2 
LI JJ 


(460) 


Similarly  the  approximate  optional  ratios  for  Tests  8,  9 are  shown  geometrically 


Geometric  meanings  of  the  failure  tensors  (individual  effects) 
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Computation  uaing  the  TENSOR  POLYNOMIAL  CRITERION  (of  2nd  order) 


E.  .0.0-  + F.CT.  < 1 

ij  » J 11 

F.  .a  .a.  + F.cr.  =1 

ij  i J ii 


no  failure 
failure  occurs 


i,  j = 1,2,6 


Plotting  the  equation  F.  .<7.0 . + F.a.  = 1 in  the  a,,  <7-,.  <7.  space,  one  will  obtain  a 
ij  1 J 11  12  6 

failure  envelope  in  the  stress  space. 


Let: 


'1 


= R cos 

a 

cos  P 

a 

R sin  a cos  g 


= R sin  p 


where: 

, 2 2 

R = [tTj  + a2  + 

? 1/7 

a6  ] 

(T 

-1  2 

a = tan  (— ) 

b 

g = tan 


-1 


2 2 
+ (r2 


(463) 


(464) 


(465) 
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Thus,  ’ 

we  can  rewrite  F. 

f,  + Fi 

(T- 

= i 

as 

follov 

781 

1 

J i 

l 1 

— i 

[ -i  -h— t~  — 

AR  2 + RR  = 1 

1 

where 

K = (F 

2 X , 
jj  cos  a + I 

F'12sin2o( 

+ F22 

sin 

2 a)' 

cos 

2e 

+ F66: 

. 2 

sin  3 + (Fj^cosa 

(466) 


26s 


B = (FjC08a  + F2sina)cosfl  + Fusing 

Solving  Eq.  (466)  for  R,  taking  positive  root,  we  obtain 

£ 


(468) 


R*  = 


-B  + \ B +4 A 
2A 


(469) 


For  an  applied  stress  a.,  the  stress  and  strength  vectors  can  be  computed  from  Eqs. 
(465a)  and  (469)  , respectively,  i.  e.  : 

- R (from  Eq.  (465)) 

3F  - R*  (from  Bq.  (469)  where  a and  g are  from  Eqs.  (465  b,  c)  ) 

If:  d < 7 ^ no  failure 

If:  A 7 ^ failur  e occurs 

CT^ 

To  compute  the,?  vector  for  a given  stress  ratios,  e.  g.  — , — , a and  g are  computed 

CT1  CT1 

from  Eqs.  (465b,c)by  setting  ffj  = 1,  and  f can  be  obtained  from  Eq.  (469). 

Caution:  The  signs  of  the  angles  a and  g should  be  set  according  to  Figure  90. 

b.  Computation  of  the  Failure  Tensors 

Using  the  strength  given  in  Section  IX.  2 and  Eqs.  (457-458), 

Fi  = x]  ' xy  = 74?  " To?  = '°‘ 003  TT^T 
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F rom  Eq.  (459  ) 

** 

B .i— !i 

7 

1 ft? 

*2  -12 

. 8 

1 J • O Ct 

2 

^ 2 

(B  Fn  + 2BF 

12 

+ F22)ff2  +(BF!  +F2)?2  = 1 

2 

■V  « / -»  w t e notn  . -3... 

12 


+ 

(13.  82 ( - < 

).  003) 

+ 0.  105) 

- 12.8  = 1 

Solve  for  F 


12 


=0.20  x 10 


3 1 


(ksi)2 


c.  Application  Examples 

(1)  Given  stress  a.,  interrogate  failure  condition  by  tensor  polynominal 


CTj  = 50  ksi 
<T ^ =-15  ksi 
= 5 ksi 


■ 
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(2)  Given  strain  e.,  interrogate  failure  condition  by  tensor  polyminal 
e = -5.5  x 10  ^ 

, -3 

e =6  x 10 

-3 

= 2 x 10 


ff2  a6 

(3)  Given  stress  ratios  — = 0.  5,  — = -0.2,  CT,>0  compute  the  strength  in  this  direction 
ffl  °l  1 
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5.  TRANSFORMATION  OF  FAILURE  CRITERION  FROM  ONE  ORIENTATION  TO  ANOTHER 
ORIENTATION. 


f(«T>  iW) 


Method  1 

For  maximum  stress  or  maximum  strain  criterion,  transform  stress  (or  strain) 
to  reference  direction  where  failure  criterion  is  known. 

Ol' 


* =>! 


Apply  failure  criterion 
Eq.  (451)  or  (455). 

Only  one  state  of  stress 
can  be  evaluated  per 
operation. 


Method  2 

For  tensor  polynominal  failure  criterion,  either  use  method  1 or  transform  failure 

tensors  from  reference  direction  to  desired  direction, 
i 2 


■ 


F 

F.‘ 

1 

1 

F.'a.'  + F .'cr.'cr  ' s 

F.. 

ij 


F..' 


iJ  1 J 


(47C) 


AJ 

1 

1 

Reference  direction  Arbitrary  direction 

Any  state  of  stress  in 

X 

► 1 

i 

orientation  X.'  can  be 

evaluated 

- -f- — -4  -4-  — t — -4- - -4 *- --  j ■— — J 

L J L.i.  1 1 1 - 1 1 _L J 1 

. 1 
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LAMINATE  STRENGTH 

1 1 

nr 

i 

i _ j 

a. 

\ 

Fir 

8t  ] 

Ply 

Failure 

: 

i — \ 

To  assess  laminate  strength  two  approaches  can  be  used 

(1)  Assume  laminate  as  homogeneous,  following  method  outline  in  Sections  3 and  4, 
but  replace  the  lamina  failure  criterion  by  a laminate  failure  criterion:  i.  e.  , by 
appropriate  numerical  value  for  F^,  F.^  and  higher  order  terms  when  appropriate 

(2)  Use  laminate  analysis  as  discussed  in  Section  V to  obtain  stress  in  each  layer 


and  interrogate  layer  by  layer  using  pro 

cedure  outlined  in  Section  IX.  5. 

In 

i thj 

is  c 

:ase  stress  vector 

' and  strength  vector  for  the  kth  laver  are  resDectivelv 

i 

— j 

: 

> 

f(k 

j. 

_ <k)  2 

(k),2. 

1/2 

! ] 

J 

ff2  ) + 

r 1 

/2 


(471) 

(472) 


where 

a.: 

i 

are  the  failure  stresses  of  the  kth 

F.(k)a 

1 

*<k>  + F (k)o  *(kl 
i ij  ui 

1 — — r 4 — }— • 

'a.*™  = 1 

J 

and  F 


(k) 


layer  from  the  roots  of 


F.^  ^ are  the  failure  tensors  for  the  direction  of  the  ktk  layer 


(473) 




l 

J 

[ 

determined  by  transformation  as  described  in  Section  III.  4. 

Repeating  this  for  all  the  layers  in  the  laminate,  the  potential  failure  layer  can  be 
seen  graphically  in  representation  such  as  in  Figure  91. 

In  Figure  91a,  the  layer  0^  ^ is  closest  to  failure  and  layer  0^  has  the  greatest 
margin  of  safety.  Further  loading  would  lead  to  first-ply  failure  in  ff  ' , 

Note  that  margins  of  safety  are  in  general  not  equal  for  all  laminates.  The  layer 
thickness  and/or  orientation  may  be  varied  to  achieve  a uniform  failure  condition 
as  shown  in  Figure  91b  which  is  the  optimal  design.  The  methods  tor  varying 
these  parameters  are  sizing  and  mathematical  optimization. 

The  procedure  of  determining  the  strength  of  a laminate  configuration,  i.  e.  , 0^ 
(k) 

b , undergoing  applied  loads  Nj,  N2>  and  is  following: 

(l)  Compute  Q..(k),  F..(k),  and  F.(k) 

ij  ij  i 

(ii)  Compute  A..  = Vq.  /k^h^ 
ij  ^ ij 

(iii)  Compute  A..  * 

o f-.L  j -1 — I . 

(iy)  Compute  a.  = 
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(v)  Compute  0.^  = ° 

o (k) 

(vi)  Compute  = R (from  Eq.  (465a)) 

(vii)  Compute  JC  ^ = R*  (from  Eq.  (469)) 

(k)  (k) 

(viii)  Comparing^'  to  to  determine  whether  failure  has  occurred  at 

the  0^  layer 
1.  Example 

Given  [457-45/0/0^:  h(°  * + h(45)  + h(_45>  = 0.  200  inch 


Nj  = 10. 

i 

N2=  »• 

i (klbf/in) 

N,  = 4.  1 

6 j 

1 
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b.  Behavior  After  First  Ply  Failure 

"Failure"  o i a ply  (any  mode)  in  a laminate  "degrades"  the  laminate  but  may  not  produce 
ultimate  failure.  Several  schemes  have  been  proposed  to  account  for  this  degradation 
including: 

(a)  Total  ply  discount  - assign  zero  stiffness  and  strength  to  the  failed  ply,  all  modes. 

(b)  Mode  limited  discount  - assign  zero  stiffness  and  strength  to  transverse  and  shear 
modes  if  ply  failure  is  in  the  matrix  phase:  if  fiber  phase,  discount  all  modes  as  (a) 
above. 

(c)  Assign  residual  properties. 

(1)  Example 

Consider  a laminate  [90,  0 , + 45,  -45]  subject  to  N.  = [3380,  1320,  0]  lbs/in. 

and  t = 0.  005  in.  all  layers. 


with  Q. . = 

VJ 

20 

0. 

0.  3 

1 

0 

0 

0 

0 

0.  5 


x 10  psi 


[e  } = { e ]°  = [A]  1 [N]  all  layers 

x,y  x,  y 1 x,  y 


Assume  a membrane  state  of  stress  exists  such  that 

rs 

Use  a maximum  strain  failure  criterion  with 
= 0.  010,  -0.  010 
E2  = 0.005,  -0.007 
Efe  = 0.015.  -0.015 


LAYER 
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6 2 

6 12 

M.S.  (1) 

31 

0.006 

0.  002 

o 

0 67 

1.5 

mmmm 

inT/vZ 

0.  006 
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4.  0 

-0.  17 
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*iP 

0.004 

0.004 

-0.  004 

1.  5 

0.  25 

2.75 

B3H 

0.004 

0.004 

0.004 

1.  5 

0.  25 

2.75 

A transverse  tension  failure  is  predicted  in  the  90'  ply. 

Following  degradation  scheme  (b) 


Q 90  -Q  90  = Q 90 
22  U12  U66 


n . 90  _ 90  „ 

0 . e2  - e6  =0 


Recalculate  [A]  and  { e }* t 2>  evaluate  M.  S 
Problems 


Confirm  and  complete  the  above  analysis 
Repeat  using  the  first  scheme 
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{ « 3k  = { e°)  + *k{  k } 

2.  Composite  laminates  are  'statically  indeterminate1  systems  in  one  or  both  of 
two  senses: 

(a)  structural  - a structural  member  has  boundary  conditions  or  geometry  such 

that  the  constitutive  equations  are  necessary  in  determining  stress  resultants 
at  any  location. 

(b)  lamination  - the  determination  of  stresses  in  any  given  layer  requires  a 

knowledge  of  the  laminate  constitutive  equations  and  each  layer  stiffness, 
ex.  homogeneous  materials  N's,  M's  <r's 

layered  materials  N's,  M's  e 's  -“^tr's 

3.  Each  cycle  of  a stress- strength  analysis  involving  degraded  laminate  properties 

will,  in  general,  require  a new  stress  analysis  for  { N,  M}. 
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7.  SIZING  FOR  STRENGTH 


1 


In  Section  6 , we  outlined  the  procedure  for  checking  the  strengths  when  the  lamination 
configuration  is  given.  In  this  section,  we  discuss  the  methods  for  determining  the 
lamination  configuration.  The  lamination  configuration  consists  of  the  orientation  and 
thickness  of  each  lamina. 


J th  j(k) 

Using  to  represent  the  failure  condition  of  the  * layer  lamina,  i.  e. , ■ ^ „ , =1  is  the 

^ sJf1  k>  (k) 

state  of  imminent  failure,  whereas  ■*  . <1  is  a no-failure  condition.  When  the  lamination 

configuration  is  not  optimized,  as  in  Figure  91a  , each  layer  would  be  at  different  degree 

from  failure.  For  example  layer  8^  is  close  to  failure  and  has  the  greater  margin  of 

safety.  It  can  be  seen  that  the  orientation  of  the  individual  layers  or  their  thickness  or  both 

can  be  varied  such  that  they  will  have  the  same  factor  of  safety.  This  is  illustrated  in 

Figure  91b  and  is  the  optimum  configuration.  This  optimum  configuration  cannot  be 

obtained  explicitly  for  the  reasons  to  be  explored;  it  can,  however,  be  determined  by 

formal  optimization  and  nonlinear  programing.  We  present  several  direct  sizing  methods 

for  estimating  the  optimal  configurations. 


h h 


Figure  91  Stress  RatiOj^k'/jC^  versus  the  Thickness  h^k'. 
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Computation  for  Laminate  Strength 


h3 

• 

Obtain  laminate  average 

Q..(k)  = Q.^tQ.  °,8  ) 

H tj  ij  k' 

strain  by: 

(474) 

• 

• 

• 

(II) 

• 

or  (III) 

■V.  = A..(Q..(k),h(k') 

H ij 

(475) 

Apply  failure 

or  Apply  failure 

(476) 

criterion  in 

criterion  in 

strain 

stress 

N. 


iJ 


(477) 


g(k)(e.°)<  >1 


f((T.k)<  > 1 
1 


o In  order  to  compute  a ; h(k\  must  be  known 

ij 

o However,  in  design  where  only  load  carrying  capacity  N.  is  given,  a = A (h^ 

(6)  . i ij  ij'  ’ ’ 

Qjj  ) is  generally  not  explicit,  the  determination  of  thickness  requires  either 

mathematical  optimization  or  direct  sizing. 

Sizing:  Strength  at  Minimum  Weight  (no  pre  existing  Haws) 


Given  a generalized  loading  configuration,  determine  directly  the  composite  lamination 
configuration  for  minimum  weight. 

Definition  of  Problem 

Given:  , N^,,  , find  h^  and  0^  (thickness  and  orientation  of  k*1 11  layer) 

(1)  Method  1:  Assume  total  decoupling  of  layers  and  fiber  direction  carry  the 

load  completely.  I 


I 


I 


Limitation  for  Efiber»  Ematrix  (example,  fiber  reinforced  elastomer). 
Method  1A:  Use  0°  and  90°  layers  exclusively* 


Applicability  useful  for  complex  but  homogeneous  state  of  stress 


xi  1 


= -!,  if  Nj  <0 


b(,0,=^ifK2>0 


= Xj'  lf  N2  < 0 


° For  E2  > El 

^1*  are  longitudinal  tension  and  compression  strength  of  0°  lamina 


For  > E2 


X1  = E2  °11  ' X1  = Q11 


Method  IB:  Same  assumptions  and  limitations  as  method  1A. 


use  0°,  90°,  ± 45  layers 

Applicability  useful  for  complex  and  nonhomogeneous  state  of  stress. 


Figure  93  Method  IB. 


199 


— 

, c 

N1 

n 

□ 

n 

I 

90 

N2 

' — 

n 

X1 

if 

N1 

> 

[) 

1 

* 

X1 

,f  n2 

>0 

M 

= 

N1 

if 

N, 

< 

0 

— 

__2 

if  N 

< 0 

X1 

1 

X1 

2 

N, 

h(- 

•'  ? 

w 

c 

X, 

, if  N 

6 s 

0 

* 7 = 

- ifN6>C 

1 

N£ 

. 

N6 

"X1 

11  IN 

6 < 

0 

K 1 6 * 

Foi 

E_ 

> 

E, 

..  .1  J 

z 

1 

xr 

X 

j'  are 

longitudinal 

tension  and 

compression  strength  of 

o 

For  F. 

> 

E 

i 

2 

X1 

= E2 

D11 

, > 

i 

1 

= 

E2  Q11 

Discu 

isions 

(479) 


(480) 


° Equivalent  to  netting  analysis,  matrix  contribution  in  stiffness  and 
strength  are  ignored. 


o Geometric  interpretation. 

I1  - I ■ ■ 1 

,90  fadune 


t — r 


-r 


envelope 


— 0 "failure 

•n  VtUP 


4 


727&Z/?\ 


So  -fa.  dure 
en*Uj>d 


1 


0°  failure 


S>*/1 

Fieure 

€ 

k 

B < 

1 

6^ 

wot 

CorvServmTi  ✓«- 

4 

L 

— 

o 

For  composite  which  does  not  satisfy  criterion  foi 

the  actual  tensile  and  compressive  strength  must 

being  conservative, 

he  replaced  by 

T 

1 ult  1 1 

‘ 

xi ' •**,%! 





Explicit  £ 

Sizing 
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listing  flaws 

) 

( 1 ) Method  2 

This 

method  accounts  for 

the  interaction 

betw 

sen 

layers  by  using  line 

ftr 

laminated 

plate  analysis.  Explicit  relationship  are  presented  if  lamination  orientation  are  pre- 
selected . 
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n total  number  of  layer 

Q.k  stiffness  of  k**1  layer 
ij 

h thickness  of  k layer 
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Stress  to  be  carried 
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strain  criterion) 
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n x 1 matrix  Thickness  of  each  layer 


To  solve  for  [H],  let  [K]  = [E]  [Q],  combining 
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Inputs  to  Eq.  (491) 

D. . Stiffness  of  0 lamina 
N.  Principal  stress  from  given  N. 

e i are  determined  from  Figure  94. 

for  Nj,  > 0 use  point  A,  for  Nj  < 0,  > 0 use  point  B 

N^,  ^2  < ® use  P°*nt  f°r  ^2<  > 0 use  point  0 
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SECTION  X 


FRACTURE  TOUGHNESS 


1.  ELASTIC  STRESS  ANALYSIS 

a.  Stress  Intensity  Factors 

Kj  = <r  \J  ira  , 


K 


II 


"a 


(499) 

(500) 


b.  Crack  Tip  Stresses 
K, 


I 


■vfijrr 


Re 


/ ^2 

1 l 


Figure  95  Reference  coordinates. 
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+ — " Re 
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0j  = cos  9 + |o,j  sin  9 , 0^  = cos  9 + (jl^  Sin  9 

p.j  and  are  the  roots  of  the  characteristic  equation 

S11M’  ‘ 2S16M'  + (2S12+ S66^  " 2S26M’+ S22  = °- 


(501) 


(502) 


(503) 


(504) 


(505) 


Note  that  u3  = R.,  n4 
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d.  Energy  Release  Rate 
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e.  Reductions  for  Orthotropic  Materials 
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Figure  96 


isotropic  and  orthotropic  finite  width  correction  factors 


rvw  * 

= 0 

xy 

2.  FRACTURE  TOUGHNESS  OF  UNI  Dl  RECTI  ONAL  LAMI NAE 


a.  Crack  Parallel  to  Fibers 
(1)  Mode  I and  II  Loadings 


o 


= Kuc 


(5.191 


Crack  Half  Length,  in.  Crack  Half  Length,  in. 


Figure  97  Mode  I loading,  Scotchply.  Figure  9 8 Mode  IX  loading,  Scotchply  . 
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(2)  Mixed- 
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Kj  i / xlO  , psirZn 


Figure  99  Interaction  between  K and  K , Scotchply  (E  = 34.  5 GPa, 
Et=  11.5  GPa).  [23]  1 u L 

b.  Crack  Normal  to  Fibers 
(1)  Crack  Tip  Damage 

Composites  with  brittle  matrix  - longitudinal  cracking 


f f 1 
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Y 

K|C  (90) 
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Figure  10  0 Comparison  of  crack  tip  stress  ratios  with  strength  and  fracture 

toughness  ratios  for  unidirectional  Gr/Ep.  Elastic  moduli  are  ET  = 145  GPa, 
Et=  11.7  GPa.  Glj=  4.48  GPa,  * = 0.  21.  [24] 


Composites  with  ductile  matrix  - longitudinal  plastic  deformation 
Elastic  crack  opening  displacement  (COD)  at  the  center  of  crack 
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Plastic  crack  tip  opening  displacement  (CTOD)  [25] 
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Figure  101  Crack  opening  displacement  in  [0]Q„  B/Ai  -6061,  E = 245  GPa 
>=5.2.  [24]  8T  L 


(2)  Notch  Sensitivity 

Longitudinal  cracking  reduces  the  notch  sensitivity, 


AS- CUT  • 
SHARPENED  O 


(a)  X = 690  MPa 


Figure  102  Effect  of  interfacial  bond  on  notch  sensitivity  of  carbon/epoxy 
(a)  surface  treated:  (b)  surface  untreated.  [26] 
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(3)  Work  of  Fracture 

W s work  accompanying  fracture  / unit  crack  extension 
Matrix  fracture 


W1  * Wm(1-V 


Fiber  fracture 


W2  = Wfvf 


Fiber  / matrix  debond 
W3  = 2Wd(L/df)vf 


Work  after  debond  [27] 
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Figure  103  Typical  fracture  mode  in 
unidirectional  composite. 
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3.  FRACTURE  TOUGHNESS  OF  MULTI  01  RECTI  ONAL  LAMI  NATES 


a.  Crack  Tip  Damage 

(1)  Nondestructive  Examinations  (NDE) 


Transmitted  light  - translucent  composites  (Gl/Ep) 
Radiograph  with  tetrabromoethane  - Gr/Ep,  B/Ep 
Dye  penetrant:  regular,  fluorescent 
Ultrasonic  C-scan 
Reflected  light  - smooth  surface 


Frgure  104  NDE  of  crack  tip  damage  in  0)/±45^2a  Gr/Ep. 


(2)  Typical  Crack  Tip  Damage 

Subcracks  along  fiber  directions 
Delamination  between  subcracks 


215 


T f 


-f- 


R R -k  4 -©  I 
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[O/  + 30I  |0/*l)5l  [0/?60]  [0/90]  [O/tl.5/90] 

Figstr.  105  Typical  subcrackB  in  Qo/±e3anid  Q>/±45/9(0  laminates. 


(3)  Size  of  Crack  Tip  Damage  Zone  (c) 

, . 1 I ....  —4 — . — 4 )- — -j 1 : — 1 — * 

Direct  measurement  - projected  length  along  the  crack  plane 
Compliance  match  - compliance  based  on  COD 
Plastic  zone  analogy  - analytical 
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Compliance  biased  on  the  diaplacement  between  two  points  3. '8  mm 
and  below  crack  surface.  The  material' is  Gr/Ep.  [29] 
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Figure  107  Comparison  between  analysis  and  method  of  compliance 

match:  .(a)  £q/±*5)b  Gr/Ep;  (b)  [o/90/±45]]s  Gr/Gp.  [24.  29] 


b.  Resistance  Curves 


*eff  ' aeff  = a + c 


kr  = Y'V"%ff  • 

c is  determined  by  the  method  of  compliance  match. 


c.  Prediction  of  Notched  Strength 
(1)  Final  Fracture 

I ( 1 t 1 | I . , | ‘ , 

K = YX  V*(a+c  ) 


-4-4 


c cab  be  interpreted  as  an  effective  Critical  damage  zone  (ECDZ)  size. 
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(538) 


(5  39) 


(540) 


(541) 


Figure  108  Resistance  curves  for 

Q>/±40  Gr/Ep.  Fracture  stress  is 
determined  from  K = Kj^  and 
8K/ ?a^dK^/da.  For  the  com- 
posite shown,  the  terminal  point 
of  corresponds  to  the  final 
fracture.  [29] 
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TABLE  67  ECDZ  SIZE  AS  DETERMINED  BY  EQ.  ( 545)  [24] 


Material 

Laminate 

Gr/Ep 

T300/5208 

[0/±l*5]s 
[0/±l*5]os 
[0/90/ ±45 ]s 
[ 0/ ±u  5/90  ] s 
[0/90]^ 

Gr/Ep 

T300/931* 

[±i45/o/90]s 

[90/0/*U5]S 

B/Ep 

I l02/ih5]s 
[0/+l*5/0/-l*5]s 

1 [ ±U5/ O2 ] s 
l [0/+l*5/0/-lt5]3s 

r 0/ ±l  5/90  ] s 

[0/±h5/90]Us 

[0/90/±l*5]Us 

Gl/Ep 

Scctchply 

[0/±45/90]2s 
[ 0/ 90 ] h s 

B/Al 

606l-F 
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3.279  52 
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Figure  109  Variation  of  c with  crack 
half-length.  [29] 


(4)  Angle  Crack 
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Figure  110  Geometry  of  specimen  with  Figure  111  Normalized  fracture  strength  vs 

an  angle  crack.  normalized  equivalent  crack 


Figure  112  Comparison  between  normal  CTack  and  angle  crack,  Qr/Ep. 


FATIGUE  AND  LIFE  PREDICTION 


UNNOTCHED  FATIGUE  BEHAVIOR 


Unidirectional  Laminae 


(1)  Longitudinal  Fatigue 


Cycle*  to  Failure 


Figure  113  Constant  amplitude  0°  fatigue  of  B/Ep  at  RTD.  R=0.  1,  -1,  10 


X=1 33 1 MPa,  X 


Failure  mode 


Frequent  tab  debond 


Failure  initiation  in  the  form  of  Longitudinal  cracking 
Residual  strength  degradation  minimal  until  immediately  before  final  failure 


No  change  of  modulus  in  B/Ep,  Gr/Ep,  Gl/Ep 


curve  and  hence,  large  scatter 
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(2)  Transverse  Fatigue 


Figure  114  Constant  amplitude  90°  fatigue  ol  B/Fp  at  R 1 jj. 
R =0 . 1 , Y=6 1 MPa,  Y'=276  MPa.  [31] 
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Figure 

i 115  Constant 

amplitude  *45°  fatigue  of  B/Ep  at  RTD. 

R=0.  1,  5 

C =133  MPa,  S=X 

45  45 

n.  [31] 

b.  Multidirectional  Laminates 

Fiber  controlled  laminates 

Character  i a tic  A of  longitudinal  fatieme 
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(1)  S-N  Relationships 


TENSION-  TENSION  FATIGUE 


0 2 4 « • 

LOG  N 


Figure  116  S-N  relationships  of  B/Ep  laminates  having  various 

fraction  of  0 0 plies,  R=0.  1.  Maximum  fatigue  stress 
is  normalized  with  respect  to  static  strength.  [ 31] 


1234567  01234567 

LOG  N LOG  N 

Figure  117  Effect  of  compression  in  fatigue  of  Gr/Ep  laminate  .[  32] 


Figure  118  Constant  life 

diagrams  for  B/Ep 
laminate  and  A1  - 
7075-T6. 

0 0.2  0.4  06  08  1.0 


(2)  1 

Failure  Mode 

. 

Li  1 

ITailtirp  cpnupnf 

:e  i 

- -T T-  ~ r r-  7 

n ro/±45/90l  laminate 

* w — 1 

i . 

Failure  of  90°  ply— i failure  of  ±45  ° plies-^failure  of  0°  ply 
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Figure  119  Typical  failure  mode  observed  along  edge  of  [0/90] 
laminate  • 
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Figure  120  Change  of  secant  modulus  in  fatigue  of  [0/±45/90]g 


Gl/Ep,  Eq  = 20.3  GPa.  [33] 
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increase  in  fatigue  of 
Gl/Ep  Laminate.  The 
increase  is  proportional 
toS^maxf,  where  f is 
the  frequency.  [33] 
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2.  NOTCHED  FATIGUE 

a.  Failure  Mode 


Longitudinal  cracking  at  notch  tip  more  extensive  . han  in  static  tension 
More  effective  relief  of  stress  concentration  at  crack  tip 


(a)  (b) 

Figure  122  Typical  crack  tip  damage  in  [0/±45l  Gr/Ep:  (a)  static  tension; 

o 

(b)  fatigue.  Dotted  lines  represent  the  boundary  of  delamination. 


b.  Notch  Sensitivity 

Residual  Strength  increases. 

Compliance  based  on  notch  opening  displacement  increases. 
No  penalty  for  notch. 


| 


Time 

Figure  123  Increase  of  residual  strength  and  compliance  in  fatigue 
of  composite  laminate  (schematic). 


3.  APPROXIMATE  PREDICTION  FOR  FATIGUE  LIFE 


a. 


Off-Axis  Fatigue  Strength 

^L’  ST’  SS’  1 lonS’-tudinal»  transverse,  shear  fatigue  strength,  respectively. 

Sg  can  be  determined  by  applying  any  static  failure  criterion. 

For  example,  the  maximum  stress  criterion  of  Section  II.  2.  a.  yields 


Sg=min 


S^/m2,  ST/n2,  Sg/(mn) 


(547) 


Figure  125  S-N  relationships  in  off-axis  fatigue  of  [o]  Gl/Ep.  [34] 


Assume  30°  off-axis  fatigue  is  controlled  by  shear  failure  and  60°  off-axis  fatigue  by 
transverse  failure.  The  data  can  be  fit  by 
Sg  = -afl  log  N+b0 


0 

e 

ae 

be 

as 

bs 

aT 

bT 

MPa 

MPa 

MPa 

MPa 

MPa 

MPa 

30 

7.  14 

67.97 

3.09 

29.  43 

— 

— 

60 

2,0  29 

35.  85 

1.  72 

26.89 

(548) 


For  0 = 30°,  m=cos  30°,  n = sin  30°, 


SS=mn  S30= 

-ag  log  N + bg 

Vmna30' 

b =mn  b„„ 

S 30 

(549) 
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For  8=60 

O 

» 

m=cos  6 

0 1 

i=sin  60  ° 

. . 

' 

S = 

2 

:n“ 

s,  = 

-a 

log  N + . 

b 

T 

60 

1 

r ‘ i 

T 

; 

(550) 

aT  = 

£ 

:n 

a60* 

■”\o 

1 

Pr^di  ctcd 

s- 

<1  = 5 

to 0 

A / » A 

1 _ . 

\ 1 

S20 

y.  o i4*  iog  in  t y 1.0  r 

si 

iea 

r failt 

ire 

= - 

0 

C'- 

’ N + 229 

. 87 

MPa  : 

transverse 

failure 

(552) 

. 

The  firat  Eq.  gives  lower  fatigue  strength  and  hence  represents  the  S-N  relation. 


S-N  equations  based  on  shear  failure 


O 

8 

a«,  MPa 

b„.  MPa 
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35.  59 

338.96 
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172.  10 
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Figure  126  Comparison  between  data  and  prediction  for  off-axis 

" * t — - I J 4 — "l  j 'I 

S-N  relationship  of  [0  ] Gl/Ep.  [ 34} 
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b.  Angle  - Ply  f±0] 

Laminates 

r~FT  i — [ i — n i — r • — i — r- 

Assume  elastic  moduli  do  not 

-<  1 

1 1 1 

cnange  in  fatigue. 

l 1 ' * i i 

Application  of  the 

maximum  s 

itiess  criterion  to  fatigue  strengths  leads  to 

— f~ -• -f-— — 4 - 4- 

. 

Sg  = min  ' 

(j 

Sr  3 £ 

Vks) 

— 

(553) 

where 

r 

/tl  4-  -> 

A 

k - 1 

1 

+ sec  20  - 

'1  ‘ 

i-v/  \*xa  C.V 

1 554) 

L 2 

u2  + tan^ 

20 

ZlTj 

5 . 

k = — 

( 

Uj  + sec 

20)  tan  20 

T 2 

1 

- sec  20  + 

" 

u2  + tan' 

1 20 

- 

(555) 

— 

- 

S 2 


4-  TO 

w 

w 

1 

U1 

= • 

L T 

i , -»v 

/ I-l 

f 1 

1 1 

TTj 

. 

U2  = 


s6  = 'ae log  N + be 

a0  = ai/ki  1 = L-  T,  or  s 


bfl  = b./k. 
o i i 


(556) 


(557) 


V°i.t 

1 + 2V.  _ + E. 

. /E_ 

(558) 

(559) 

(560) 

(561) 


-1-  J - . i i 4 i 
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MPa 


For  Gl/Ep 


Type  of  failure 


STRESS  CYCLE  n./8  v* 
STRAIN  CYCLE  n-l 9 c/S 
STRESS  CYCLE  0*11.5  VS 
STRESS  CYCLE  n-rJ  c/% 
STRESS  Ct  CLE  n=  5 fi  r/S 


N (cycles) 


Figure  127  Experimental  and  predicted  S 


Fiber -Controlled  Laminates 


Fatigue  strength  of  laminate 


(tens  ion -tens  ion) 


x(fraction 


Ik  1 

b- 

I S' 

e 

6 

0.271 

24.76 

235.7 

0.5 

13.42 

127.8 

0.  428 

3.29 

51. 31 

el 

E^ 

T 

glt 

‘'lt 

U1 

U2 

GPa 

GPa 

GPa 

54.  7 

17.  75 

6.01 

0.285 

-0.  4475 

1. 957 

4.  LIFE  PR.EDI CTI ON  AND  ANALYSI  $ QF  SCATTER 


a.  Application  of  Failure  Potential 

: ' T T,  [ t~'j  ! I 11  1 

(1)  General  Formulation 

Probability  of  surviving  time  t under  a loading  history 


R(t)  = exp  C -0(Tj  ] . 

l r t a 

T = — j K(s(£))  d£  , t has  the  dimension  of  t. 
4 o 

Failure  potential 

. a 

«T)  = T 

Breakdown  rules  (36] 

Power  law:  K(s)  = (s/C  4 

1 1 

Exponential  law:  K(s)  = exp  (s/C  ) 

C2  3 


(2)  Stress  Rupture 
s(£)  = s const. 

Power  law  breakdown  rule 


T = “J*  (s /Cj)Pd4  = ( s/Cj)^  t/t' 

4 o 

R(t)  = exp  (s/C^^ft/?)01!  = exp  |~-(t/t0)a  1 

A O ^ J 


(to/tj  (s/Cj)^  = 1 


(568) 


or  log  (t0/t)+  (Slog  s = Slog  C.:  Power  law  (log  stress  - log  time)  representa- 
tion of  S-N  relationship 


Exponential. 

breakdown  rule 

R(t)  = exp  | - 1 ~~  1 

exp 

(s/c3)  1 

* 1 = exp  T -(t/t0)al  (572) 

t L 2 

(t0/£)  exp  (s  /C3)  = 

A /,  \ 

C2 

J > L J 

: Exponential  (log  time  - log  stress)  representation 
of  S-N  rela.tionshin 

ios(«0/f,+  (lea*.) . 

= : 

Log  C2 

(573) 

v 1 

_|  fj  H 

1 — 7 j h 4 — -4 * 4 -4  *4 4-  -4  -4 f 1 

I 

1 L.  J I l 1 
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In  the  failure  range  exp(x/C^)  « 1, 


R(x)  = exp/- 


;; exp(ttx/C  ) 

\tLC2  / 


— 3 r A , 

x = [aXnftLC^/C^)  - yj  y:  Euler  constant  (=0.  5772) 


TABLE  68  PARAMETERS  FOR  STRESS-RUPTURE  DATA  [37] 


Strand  Type 


Kevlar  49/Ep 
Gr/Ep 
S-Gl/Ep 
Be/Ep 


] 

Power  Law 

Exponential  Law 

p 

C1 

a 

C3 

C2 

MN/m 

MN/m 

42 

2239 

0.  87 

50.  5 

19 

1. 86  x 10  y 

78 

922 

0 30 

10.  8 

8.  32  x 1036 

30 

2109 

0.  75 

62.9 

14 

1.91 x 10 

26 

733 

75 

2 5 1 

and  Nonsimilarity 

( 1 ) Definitions  [ 38] 

Objects  that  follow  the  same  law  of  behavior  are  similar. 

Objects  that  follow  different  laws  of  behavior  are  nonsimilar. 

Two  objects  are  identical  if  there  is  an  exact  duplication  of  the  atoms  of  the  first  one 
in  the  second. 


DISTRIBUTION  OF 
INITIAL  STRENGTH 


LOG  STRENGTH 
LOG  LEVEL 


DISTRIBUTION  OF  LIFE 


I 1 J 

(a) 

(b) 

Figure  1 

32  Life  < 

iistribution  of  similar  objects:  (a)  similar  deterioration 

(b)  nonsimiiar  deterioration. 
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Figure  133  ! 

Life  distribution  of  nonsimilaT  objects: 

fb)  nonsimilar  deterioration. 
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similar  deterioratior 
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(2)  Relationship  Between  Static  Strength  and  Life  of  Similar  Objects  With  Similar 
Deterioration 


(a)  Pairing  of  static  strength  and  life 

Static  Strength  Data 

6 £xi5.x2£‘‘‘1  xn  > 3 : fatigue  stress 

Life  Data 

t.  < t.  < < t 

1 — 2 — — n 

• 

LINEAR  LEAST  S0UARES 

/ 

fit 

3.  —60  (3) 

EO  ( 7) 

•/s'' 

r 

IS 

glass /epoxy  ( o,) 

^ s 

APfUEO  STRESS  758  MNW* 

/* 

TIME  TO  FAILURE.  MIN 


100  (000 


Figure  134  Plot  of  (t.,x.)  for  folGl/Ep  subjected  to  constant  stress*  [39} 

(b)  Similarizing  operation  (proof  testing) 

Proof  load  to  x^.  Then  the  minimum  guaranteed  life  corresponding  to  x.  is  t. 


SIMILARIZED 

STRENGTH 

DISTRIBUTION 


LIFE  DISTRIBUTION 


Figure  135  Effect  of  similarizing  operation  on  life  distribution. 
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(c)  Mathematical  formulation 

Determine  strength  and  life  reliability  functions 

R (x  ) and  R (t) 
s l 

Relationship  between  x and  t is  obtained  from 


R ( x | s ) = R (t) 
s X 

where  R ( x|  s ) = R (x  ) / R (s) 
s s s 

In  particular,  if 

Rg(x ) = eXp[-(x/xo)as] 

R/t)  = exp  £-(t/t0)ai  ]] 


then 


= t0  |j 


(*/x0r8  - (s/xo) 


■] 


1 / a, 


(586) 

(587) 

(588) 

(589) 

(590) 


In  terms  of  residual  strength  at  t 


R ( xl s)  = R (x  ) 
s r r 


(591) 


: 


I 


where 

R_(s)  = R (t) 

l 592) 

r V 

Figure  136  Relation  between  static  ftrength  and  fatigue  life  of  [o]  Gr/Ep: 

Qtg  = 7.77,  xq  = 1.44  GPa,  0 = 0.74,  tQ  = 1.6xl06  cycles, 

Smax^xo  = Results  of  proof  testing  are  also  bhown.  [40] 


i 
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c.  Strength  Degradation  Model  [39,41] 
(1)  General  Formulation 


x : ideal  static  strength 
s 

x^:  residual  strength  at  time  t under  a loading  history  s(t) 
T:  material  age 


Introduce  normalized  variables  (nondimensional) 

x =x/x,  F = t/*?  , x=x/£ 

r r s s 


dx  - 1_“r 

r - - x . or  > 1 

r 


dT 


— a — a 
x r-5T  r - Or  T 

sr  r 


dr  = K(s,  t ) dt 
Assume 
K(s 


r — 


£ or.  - 1 


(\H 


(593) 

(5941 

(5951 

(596) 

(597) 


Assume  a Weibull  distribution  for  static  strength 


R (x  ) = exp 
s s 


•(i  r 


R (x  ) = R ( x ) = exp 
r r s s r 


Failure  occurs  when  x (t)  = s(t). 

r 


( + «r*  \V“r 

'I  ^ J 


(598) 


(599) 


(2)  Tension  With  Constant  Loading  Rate 


s ( t ) = Lt  = (Lt)  t 


, -S'  (gy  ^ 

Jo  \C\)  9 fi  + or^ 

V*/c/r««, 

• X r = X - ■ ^ t 

'•  r s P+a  t 


( Lt  / C i )P  - P + 


(600) 

(601) 

(602) 


At  failure,  x^  = s(t)/  x = (Lt/  x)  t 


(603) 
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-a,.  -ttr  020+0/ 

X = X - T 


x/  * “1 

0 + ax) 


(604) 


xf-*x8  as  L — «-x:  Xg  is  the  strength  under  a very  high  .loading  rate. 


R (x  ) = exp 
r r 


J_  /_or  + 5'+M 

“s  \ r (fl  + at)  c/(Lt)a*  r / 


V“r., 


(605) 


If  the  second  term  is  negligible  compared  with  the  first  term,  which  can  happen 
when,  e.g.,  Lis  large  (fast  loading  rate),  then 


X ex  X 


.’.R  (*  ) = exp 
r r 


/x  \ a 

i-f-) 

' OS' 


(606) 

(607) 


If  the  first  term  is  negligible,  which  can  happen  when,  e.g.,  L is  small  (slow 
loading  rate),  then 


1 

[ / x N 

(0  + o )a  /a  ' 

(x  ) = exp 

- (rH 

i 9 r 

(608) 

r 

[ V or/ 

-i 

1 

t . a , - > 0 x ®r 


(tL)  (fi  + ut)  QjT* 


°s_ 


1/0+  a^) 


(609) 


Note  that,  if  OL  = 0^  s 1 , then  Eq.  (608)  is  equivalent  to  Eq.  (579). 
(3)  Stress  Rupture 


s(t)  = s const. 

(610) 

T = (s/C  j)^T 

(611) 

R (x  ) = exp 
r r 


,|xr°(r  -t  (Of /Of)  (s/C;)^  t 011  je‘s/<*r 


(612) 


L ' 

X T ’ , 

08  J 

. A ft  8 - a , a8/ar' 

^ (s/x  )*'  + (Oy/a^)  (s/Cj)  t i ^ 

A 

x ) = exp 

V - 04.  / 

(613) 


In  the  fatigue  failure  region,  the  first  term  can  be  neglected. 
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' 

63Cp 

f f 

- 

p/a 

S'  1 

= 1, 

, tl 

t.O  (S/Cl>  = (x  08  r) 


1/al 


Note  that,  if  a = tt  si,  then  Eq.  (614)  is  equivalent  to  Eq.  (569). 


(4)  Constant  Amplitude  Fatigue 

x : realistic  static  strength 

{■■■  -1 — | — . -j — , — . — J — J — J — L.  | L — 1 

s(n)  = f(Smax,  R,  n) 

. , Smax  ,0  Of. - 1 

K(s,  n)  = ( — ) n ■* 


“r-x0^  = - — ( ^ 

r “e  Ci 


R (x  ) * R (x  | S ) •»  exp 
r r s s I max' 


/Smax\“sl 

*(—  j 

' xo  > 


. l^max^l  )^nat)Qs  /0tr 


/ N 

. \ 

<W“rl 

kjw)  » 

exp 

'(-I 

rJ 

= (*< 

VQ r 

(a)  « = « 

_ x 8 + (a  /a  )(S  /C  .fn*1  fs  \“s 

W " “P  [-  ' (^f ) 

o 


-(f)’  - (w-)“‘ ' (W* 

' O/  I T . — + J 
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Rx(N)  = exp 


P/af  a 1/a 

N (S  /C,)  = (x  s a /a  ) 

o max  1 ole 


RJxJs  = exP 

r r]  max 


CL 

x \ s 
r 


(623) 


(624) 


(625) 


li  uic  iiiak 

R (S  ) = R 

r max 

■£(N) 

(b)  a = 1 

X/ 

(x“r  + a (S  /< 

\ r r max 

a \a  /a 

~ A s r 

D / v 1 — ovn 

-1>  "/ 

X\  IX  ) — cau 

r r 

- <*s 

X 

_ o 

= expf  - 


/-  \a 

' x \ r 


R^(N)  - exp 


~ N i 

x / 'o' 

o/ 

a /a 

W* 


a /a  /_  \a 

s r /S  \ s , 

I max] 


N (S  /C.)p  = x /a 
o max  1 or 


Strength  and  fatigue  properties  of  Q)/45/90/ -45  /90/45/cT]  Gr/Ep. 
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24.  12 


x , MPa 
o 


487 

487 


24.  12 
22.  16 


a. 


1.089 

1 


P 


21.68 

19.92 


C t , MPa 


6.  153 
6.231 


(626) 


(627) 


(628) 


(629) 


x,  MPa 
6.8948 


(5)  Cumulative  Damage  Model 
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Figure  137  Weibull  plot  of  normalized  life  data.  [32] 


Residual  Strength.  MPa 


Figure  138  Experimental  correlation  of  analytical  models  in  terms 


of  residual  strength 


(c)  Experimental  correlation 


Strength  and  fatigue  properties  off  0_/±45]  B/E 
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Figure  \¥)  Distribution  of  N : theory  and  data.  [31] 


SECTION  XII 


ENVIRONMENTAL  EFFECTS  ON  COMPOSITES 


1.  HEAT  CONDUCTION 


a.  Formulation 

Heat  Transfer  (temperature  effects) 
q = - KVT 

One-dimensional  steady 

T -T 

.,12 

q = + K 

x x L 


Composite  Materials 

2 2 

= KjjCos  a + K^^sin  a (640) 
Approximations  (v^.  < 0.785) 


K11  = ;i-vf,Km+ vfKf 


* 


S K 

X \ 22  >K, 


K22  * KnS(1  ' 2VV^>  + F" 


BtC  $ Z\*i  ’ 1 


4 -i  V-(VfM 

' tan  — — — - 

h - B^Vf/TT  1 + Vs^Vf/* 


Input:  K , K a v , a 
mfi 

°Utput:  K 1 1 • K22’  Kx 


Thermal  conductivities  of  selected  materials 

K - BTU/hr/ft^/°F/in;  room  temperature 
Resins 


Filaments 


1 

Ph 

ien< 

jlic 

1.- 

4 

Carbon  fabric 

40.0 

^nenoiic,  »liuub 

1. ! 

5 

Glass  fabric 

7.5 

I 

Polyethylene 

i.' 

7 

Graphite  fabric 

1000.0 

Quartz  fabric 

1 

7.  5 

Silica  fabric 

7.  5 
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>inch  thick  fiber  reinfor 


epoxy  plate  is  exposed  to  a 120  F e 
. The  fiber  volume  fraction  is  68% 
heat  transfer  through  the  plate  for  the  following  fiber  orientations:  a = 0,  30, 


2.  MOISTURE  DIFFUSION 


a.  Formulation 


Moisture  Effects 


Steady  Bound 


Cond 


One  Dimensional 


Temp  Constant 


Uniform  Initial  Distributions 


Distribution 


Moisture 


<2j  + l) 


(weight  gain): 


Moisture  Content 


(2j+l)  TT 


Time  to  reach  99.9%  saturation 


Output:  C , c 


Maiximum  Moisture  Content 


I — 

1 — 

■ ' 

c 

1 

= a 

,Hh 

(649! 

i 

Fibe 

r — 

f y 

/ / * 

1 >X 

Diffuaivity: 

z'y 

D 

_ 

D, 

2 

cos 

» 

‘a  + D 

H 

. 2 
[in 

a l 

65C 

)) 

S y 

D,„ 

X 

11 

22 

r 1 

y y 

Approximations  (i 

V 

C 0. 

785) 

1 

i 

• 

__ 

A,- 

m'-m 

Dzz  ■ D 


- zVa 


For  Dr  « D 
f m 


— D +-T-  * 

ir  / m 


4 -1  - W" 

~ tan  

BDvf 


3X  = Dx[](I-vf)cos2 3 * *a  +(1  - 2^/vf/ir)sin2<f] 


Input:  D , D,,  v,„  a 
m f f 

Output:  Du,  D22,  Dx 


b.  Examples 

(1)  Both  sides  of  a 12.  5 mm  thick  Graphite  T-300  Fiberite  plate  are  exposed  to  air  at  350°K 

and  90  percent  humidity.  The  initial  moisture  concentration  is  uniform  inside  the  plate.  The 

initial  moisture  content  of  the  plate  is  0.  5 percent.  The  diffusivities:  D and  D are  given 

11  22  6 

in  the  accompanying  figure  (vf=0.  68).  The  constants  a and  b are  0.0014  and  2,  respectively. 

a)  Estimate  the  time  required  to  reach  one  percent  moisture  content 

b)  Estimate  the  time  required  to  reach  at  least  99.9%  of  the  maximum  possible  moisture 
content 

c)  Estimate  the  maximum  possible  moisture  content  inside  the  material 

d)  Estimate  the  moisture  content  in  5 years 

e)  Draw  the  moisture  distribution  inside  the  material  after  5 years 
Perform  the  calculations  for  fiber  orientations  a = 0,  30.  60,  and  90  degrees. 


Figure  140  Change  oi  diffusivity  with  temperature 


490  «00  990 

TCMKIUTUNC  T.H 


(2)  The  plate  described  in  the  previous  example  is  exposed  to  air  on  one  side  only.  The 
other  side  of  the  plate  is  insulated.  Repeat  the  calculations  for  this  plate. 

(3)  The  initial  moisture  distribution  ip  uniform  inside  a 12.  5 mm  thick  Graphite  T-300- 

Fiberite  1034  plate.  The  initial  moisture  content  of  the  plate  is  one  percent.  The  plate  is 

then  exposed  on  both  sides  to  humid  air  at  333°K  and  10%  relative  humidity.  Estimate  the 
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moisture  content  of  the  plate  after  10  days.  The  fiber  volume  fraction  is  v = 0.68,  the  con- 

1 r 

stants  a and  b are  0.0014  and  2,  respectively.  The  diffusivity  of  the  resin  is  D = 1.6x10 

2-1  m 
mm  s . Perform  the  calculations  for  fiber  orientations  a = 0,  30,  60, and  60  degrees. 


c.  Tests  to  Determine  cx  and  D 


Edge  effects  negligible  [42] 


for  a)  accelerated  test  procedures 
b)  edge  effects  corrections 


(655) 


Observations  (approximations) 

D ~ f(T)  (but  not  M?  ) 

cx  = f(H)  (but  not  T?  ) 
cx  = aHb 

Present  information  b = 1~2 


(656) 


T(t) 

fa 

i 

T(t) 

(tonal 

HL(t) 

\ 1 

A 

HR<t) 

Temp  same 

on  both  sides,  or 

\\ 

,/  / 

A 

yy 

One  side  insulatei 

d 

. 

r 9 

s 

/ 

A 

_Ap 

s 

/ y 

A 

y 

Numerical  Solution  ("W8GAIN") 

Input:  T(t),  HL(t),  HR(t) 

Initial  moisture  distribution 

D(T),  c_  (H),  K(T),  p,  h (for  each  layer) 

~T  pry  [ ' 

Output:  Total  Weight  Gain 

(As  a Weight  Gain  of  Each  Layer 

Fn  of  Time)  Moisture  Distribution  in  Each  Layet 


^ 4 A1 
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TABLE 

70 

SUMMARY  OF  EXPERIMENTAL  DATA  ON  THE  EFFECTS 

OF  MOISTURE  AND 

TEMPERATURE  ON 

THE  ULTIMATE  ' 
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(Continued) 


Laminate  Lay-Up  Orientation 


Composite 


Moist  Tern 


Moist  Temi 


Moist  Tem| 


Browning 


Polyimide 


Boron/AVCO 


Hofer 


Boron/Narmco 


Kaminski 


5505 


Browning 


Negligible  effect  (b)  L = Little  effect  (<30%)  (c)  S = Strong  effect  (>30%) 


Moisture  Content 


Temperature  - 
'500K 


300  350  <100  450  0 

TEMPERATURE,  K 


05  1.0  1.5 

MOISTURE  CONTENT, % 


iFigure  141  Ejects  of  temperature  and  moisture  on  strength 


of  T300/ 


3.  TRANSIENT  ANALYSIS 


Figure  142  Input  to  transient  analysis. 


Figure  143  Weight  gain  over  20  years  for  1/2" 
hybrid . 
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ACTUAL  AMBIENT 

BEFORE  FLIGHT 

-—AFTER  FLIGHT 


CONSTANT  AMBIENT 
TSF.82%  RH 


300  CYCLES 
1000  DAYS 


JBflS# 


10  CYCLES 
GO  DAYS 


1000  CYCLES! 
6000  DAYS  ' 


“400  CYCLES 
2400  DAYS 


CYCLES 


200  CYCLES! 
1200  OAY 


003  02S  O 

POSITION,  inches 


Figure  144  Moisture  distribution  at  boundary  layer 


(1-4  plies)  and  inside  1/2"  hybrid  over 


20  years  - note  flat  top  level  at  0.  42% 


after  10  years 


ACTUAL  AMBIENT 


CONSTANT  AMBIENT 


PERCENT  OF  MAX  MOISTURE  CONTENT 


0.999  Mm 


0 99  Mm 


TIME  , YEARS 


TIME,  y*or$ 


rates  of  constant 


146  Moisture 


Figure  145  Time  required  to  reach  flat  top 


and  transient  ambients  (a  guide  to 


moisture  level  under  constant 


ambient  (irtseftsitive  to  relative 


humidity). 


SECTION  XIII 


TIME  DEPENDENCE  IN  POLYMER 


DEFORMATION 


1.  LINEARLY  VISCOELASTIC  CHARACTERIZATION 


a.  Uniaxial,  Homogeneous  Deformation;  Isotropic 


Solid 


Relaxation  behavior  (prescribe 


step  strain) 


A cross-linked 
material 


Results  in 


B uncros  8 -linked 
material 


Relaxation  Modulus 


ijcally  10-20  decades 


typically 

2-3 


Creep  behavior  (prescribe  step  stress) 


Results  in 
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Example 
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Approximate 
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answer 
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(5)  Homogeneous  shear  deformation  (totally  analogous) 
Relaxation  modulus  in  shear  ^ (t) 

Creep  compliance  in  shear  J (t) 


Shear  stress  t( t)  = Zt  p (t) 


Shear  strain  2e  (t)  = 


(6)  Cyclic  deformation 
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Then 
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b.  3-D  Stress -Strain  Law  for  Isotropic  Solid  (Neat  Resin) 
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2.  EFFECT  OF  TEMPERATURE  ON  TIME  DEPENDENCE 

a.  WLF  Equation 


Log  E^(t)  ~ Absolute  temperature  towards  rubbery  domain 
Short  time  not  well  known 


WLF  Equation  Log  a 


-ci<t-tr> 


T C2  1 (T-Tr) 


b.  Glass  Transition  Temperature 


Cj  i 8.86 

C2  i 102 

T * T + 50 
R g 


AV 

V 


hard 

f 

glass 

rubbery 

leathery 

soft 


Figure  147  Determination  of  glass  transition  temperature. 
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Behavior 


LOO  tk  (tb  IN  MOUftJl 


Relationship  between  stress  and  bi 
HMDA-crosslinked  epoxy  polymer. 


i oft  14  o«  too  ftrt^Cift  MHoufto) 


Figure;  149  Master  curves  of  ff-log  t plots  of  the  HMDA 


d.  Creep  Compliance  of  Shell  58-6BR  Epoxy 


(Example  problem  in  time-temperature  shifting) 


LOG  TIME  (minutes) 


LOG  TIME  (minutes) 


Creep  compliance  for  Shell 


58-68R  epoxy  as  a function 


Figure  151  Net  creep  compliance,  AD 
for  Shell  58-68R  epoxy  at 
different  temperatures. 


TE'tPERAtURE  (*F) 


Figure  152 


Temperature  dependence  of  initial  compliance 


» i 

LJ 

. 

I i 

LOG  t/aT  (minutes) 


Figure  153  Master  curve  for  net  creep  compliance, 
AD,  for  Shell  58-68R  epoxy. 


3.  PRACTICAL  LIMITS  ON  TIME  DOMAI N OF  MATERIAL 


a.  Time  Scale  of  Structural  Problem  Determined 

b.  Time  Scale  of  Material  that  Needs  to  be  Known 

c.  Principle  of  Recent  Memory 
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d.  Strength  and  Life  Distributions 

X : strength  » F(x)  : strength  distribution 

F(x)  = P [strengths]  (700) 

X : life  ■ » F(x)  : life  distribution 

F(x)  = P [life^x]  (70I) 

Also 

F(x)  : probability  of  failure 
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2.  WEIBULL  AND  NORMAL  DISTRIBUTIONS 


a.  Weibull  Distribution  (Two-Parameter) 
(1)  Properties 
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(d)  Sample  problem 

Static  strength  data  for  [ 0 /45/90 /45/0] 


T300/934  laminate  are 


as  followb 


X(MPa) 


x(MPa) 

Fix.) 


Find  x,  s,  C.  V. , ft  , and  x . Use  the  median  rank  for  F 


23  MPa,  C.V 


From  method  (a) 


From  method  (c) 


the  strength  data,.  [ 32] 


Figure  156  Weibull  plot 
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b.  Normal  Distribution 
(1)  Properties 

r ■ :&£  "ph 

Probability  density  function 


/ \21 

M 1* 

f = 


r\f2» 


exp 


(712) 


(713) 


Mean  : 

U 

, 

Standard  deviation 

(T 

(2) 

Estimation  oi 

w 

anc 

cr 

A 

(a)  U 

= X 

» 

A 

cr 

= 

s 

(714) 

(b) 

Vlaximum  likelihood  estimate 

5ame  as  aoove 

(c 

) 

Linear  least  squares  method  ( graphical  method 

Define  - 

£ - 

jj. 

= 

y 

- 

t - 

AL 

= 

V 

(715) 

- V 
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lv 

(716) 
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(d)  Sample  problem 

For  the  data  of  the  previous  sample  problem  2.a.(2).(d), 
determine  u and  <r  . 

From  method  (a),  M = x = 478  MPa,  a = s = 23  MPa. 
For  method  (b),  first  determine 


i 

l 

2 3 4 

5 6 7 8 

9 10  11  12  13 

y 

- 

l.S 

11  -1.50  -1.25  -1.05 

-0.90  -0.76  -0.63  -0.52 

! -0.40  -0.30  -0.20  -0.10  0 

i 

14 

5 

16 

17 

o 

(\J 

o 

00 

t— t 

21  22  23  24  25 

y 

O 

O 

o 

20 

C 

. 3( 

0.40 

0.52  0.63  0.76 

0.90  1.05  1.25  1.50  1.91 

The  linear  least  squares  method  yields 
U = 477  MPa,  cr  = 24.3  MPa. 


Figure  157  Strength  data  plotted  on  normal  probability  paper. 
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3.  RELIABILITY 


a.  Definition 

The  probability  of  a successful  operation  of  the  device  in  the  manner  and  under  the 
of  intended  use. 


**’)  : distribution  function  with  parameters  0^ , 0^,  ••• 
xj  - x2  : limits  defying  the  event  : Success. 


Reliability  function  R 

"v.. 

R = P(xj<X<x2)  = F(x2;0l,02>  •••)-  F(xl;ei>02,  •••) 
= R(01.e2,  •••  ; *J.  x2) 


b.  Reliability  for  Strength  and  Life 
X : strength 


R(x)  = P [ strength  >x]  = 1 - F(x) 

X 

: life 

R(x)  = P[  life  > xj  = 1 - F(x) 


c.  Estimation  of  R 


H = *<W  -IV  x2, 


An  estimator 

A , 

* 9 is 

unbiased  if 

EfSl 

e 

? 

J he  unbiased 

l estimator  of  < 

T 

is 

s'2 

= 

L 

S 

n/(n-  1) 

An  estimator 

A , 

9 IS  consistent-  if 

A 

line  ] 

p n 

9 - 

0|  > e] 

~ 

0 for  arbitrary  < 

E . 

n 0 

An  estimator 

should  have  ai 

9 small  a 

variance  as 

possible  for  any  sample  size  n 

• t - > - +■ — H 

, 

1 

1 . 

conditions 


(718) 


(719) 


(720) 


(721) 
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4.  DESIGN  ALLOWABLES 

a.  Definition 

"A"  allowable  x^.  Probability  is  95%  that  at  least  99%  of  the  distribution  will  be  contained 


within  the  interval  (x^,  x ). 


’B"  allowable  xg  : Probability  is  95%  that  at  least  90%  of  the  distribution  will  be  contained 


within  the  interval  (x  ,»). 

B 


b.  Weibull  Distribution 
Suppose  & is  known. 


-i  Aft 
2n  x 


has  the  V distribution  with  2n  degrees  of  freedom: 


- A a 
2nx 


< y 


1 r fn-l  -t/2 

“ | t e dt 

2n  I ' ( n)  Jo 


(722) 


Given  the  confidence  level  Y , let  XZ  be  defined  bv 

2n;y  ’ 

P [xl  <xl  ] = Y 

L 2n  2n;  y J 


.‘.P 


2nx  , 

2-  < x2 

a “ 2nr  y 


= Y 


= p 


2n£ft 

° n 

—? - x 

V2  ° 

L*2n;  y 


(723) 


(724) 


(725) 


O 

y 

/ 7 

,1/a 

A 

X 

y 

-hr 

“ * 

I O 

\X*n;y  / 

Since 

Xol 

< x 2 implies  R < R , 

the  lower  confidence  limit  for  R,  Rv  , is 

A 

1 / x\  " 

K 

y 

= exp 

[-(r) 

. V J 

(72b) 


(727) 
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For  "A"  allowable  x 


For  'B"  allowable  x 


n • NO.  OF  SAMPLES 


COMPOSITE 


SHAPE  PARAMETER  a 


Figure  158  Relation  between  "A"  or 


allowable  and  shape  parameter 


Normal  Distribution 


K is  chosen  in  such  a way  that  the  probability  is  y that  at  least 
distribution  will  be 


a fraction  R of  the 


contained  within  the  interval  between  x 


’A"  allowable 


For  "B"  allowable 


i C 

>•95 

(731) 

En 

=*  l 

■ I 

- 2nf  nRy  /x  y 

= [ 

1 V * i i/a 

T 2-  x-  J 

n i J 

i=l 

1 r- 

Figure  159  One-sided  tolerance  factors  for  normal  distribution,  and  K^. 
d.  Sample  problem 

For  the  data  of  sample  problem  2.a(2)(d),  determine  the  "A"  and  "B"  allowables. 
For  the  Weibull  distribution  assume  3 = 2 3.4 


Solution 

• 

For  the  Weibull  distribution 

j , 

2 5 

>*> 

II 

, 1 23# 4 ,1/23. 

[ Tg  L xi  ] 

4 

488 

SJ 

25  i=l  1 

’ j ; ■ p | 

| 

1 L J L 

LJJ 

±±±] 
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r~ 

— 

TT 

From 

x2 

distribution  table 

— 

■ — 

X 

2 

50; 

0.  <5 

= 6" 

.5 

.* 

XA 

= 

t- 

50  1 n 

o.« 

59/67. 

5] 

1/2 

3 .4 

488 

= 

396  MPa 

X1 

= 

[- 

50  /n 

0.' 

50/67. 

5] 

1 / £ 5. 

488 

= 

438 

MPa 

For  the  normal  distribution 


x = 478  , s'  = 23.5 


From  the  figure 

K.  = 3.15  K 

A i 

= 

1.1 

35 

.*.  X.  = 478-  3.  15  x 23.5 

_ 

4 

i04 

MPa 

A 

K = 478-  1.85  x 23.5 

£5 

= 4 

■35  MPa 

5.  STATISTICAL  INTERPRETATION  OF  FAILURE  PROCESS 

I — — L H 

a.  Typical  Failure  Process 

Failure  rate  at  time  t : the  probability  of  failure  per  unit  time  after  having  survived  to 
time  t. 


Time 

Figure  160  Typical  failure  rate  vs.  time. 


Initial  failure  period  : Break-in  period 

Failure  due  to  initial  defect  or  weakness. 


Chance  failure  period  : Failure  due  to  unusually  severe,  unpredictable,  and/or 

unavoidable  environmental  conditions. 


Wear-out  failure 

period:  Failure  due  to  fatigue,  aging. 

b.  Mathematical  Representation  of  Failure  Rate 

Aft)  : failure  rate,  hazard  function,  risk  function 

f(t)dt:  proportion  of  the  initial  population,  which  will  fail  in  the  time  interval 

(t,  t 

-Hit) 

R(t)  = l-F(t)  : 

proportion  of  the  initial  population,  which  have  survived  to  time  t. 

A (t)  = 

f(t)  _ 1 dP  _1_  <*R 

R(t)  = 1-F  dt  R dt 

~ = - Aft)  dt  fnR  = - /Ad, 


R(t) 

exp  [ - f A (f)dfl 

R(o) 

r 1 l » 1 

1 1 

i L I 1 1 

' : 

. 

1 — ) — 1 — p ■ j — j ; — * — * — t ~ 

(732) 

(733) 

(734) 
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For  Weibull  distribution 


1 61  Failure  rate 


and  probability- 


density 


Weibull 


distribution 


normal  distribution 


Figure  162' Fail  i 


distribution 


c.  Concept  of  Failure  Potential 


l 

) 

’ 

I 

I 


A(t)  - — \jj  : failure  potential 

Introduce 

r : material"age" 

r depends  on  the  loading  history  s(£)  through 


r = f K(s(f))df 

Jo 


K : breakdown  rule 


A (t)  = 


dip  d r 
dr  dt 


(7391 


(740) 


(741) 


R(t)  = exp  [-'I'  (r)  ] 


(742) 
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6.  SIZI 

I EFFECT 

1 ! 1 

: i it: 

] 

r~rr 

T~ 

a. 

Simple  Cha 

Ln  or  £ 

series 

VI  ode 

Consider 

a 

chain 

consisting  of  N identica 

. links 

RJ»)« 

reliability 

of  link 

, i.  e 

, probability  of  strength  exc« 

reding  tr. 

R(°)  : 

r 

eliability 

of  chain 

The  chain  tails  if 

at  least  one 

link  fail 

6. 

R 

- 

R 

N 

o 

*.  F 

c 

1 - R 

= 1 

-rn 

= l-U 

. i 

(743) 

If 

R 

s 

eKp 

W,ro)flt 

o 

J 

l '****/ 

R(<r) 

a exp  l-NJoVcr  )<* 

= 

exp 

-((T 

. Q 

/<roN> 

. 

(745) 

1 

r 

1 

j 

a 

= a 

-l/« 

N 

/ 

oN 

O 

i n 

i 

/ 

y 

£ 

/ 

/ 

I 

/ 

/ 

/ 

g 0.5  - 

/ 

“ 

i 

y 

'i  / 

y 

Strength 

i — 

Figur 

e 163  . 

_u 

Comparison  in  probability 

of 

failure  between  single  link 

and  20  - 

link  chain. 
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b.  General  Chain  Model 

Homogeneous  material  : uniformly 
distributed  flaws 


cr  is  uniform  in  each  element. 

Av. 

Ri  • R<M 


N 

lAYi 

R 

ri 

R 1 

;<ri> 

^ 

5 

i=i 

l 

o 

Jc  \ 1 

r 

\ 

N r t 

fnR  = T AV.  fn  R (<r.) 

i=l  1 L° 
fnR  = / fn  [RJa)]  dV 

• R = exp  j ' f fn  [RQ(a)]  dvj 


For  Weibull  distribution 

Ro(<r)  = exp  [^-(cr/o-^J 


fn  [ro(<t)J 

R = exp 


* -(<r/<r  ) 


[•  / 


R is  the  probability  of  surviving  the  stress  distribution  <r . 
Let 

o-  a X f(x)  x : position  vector  within  the  body. 

X:  reference  stress 


R(X)  = exp  ^(X/Xo)°  J 

1 1 4 tkf  t t i~ 


i_^-L  14--LI 
M M 
-LU  J--M-  - 


; charictferistid  strength 

1 . I 


: : 

I — — r 


I ri  i 

i i 


-J L 4_ 

TTTT 
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(1)  Sample  Problem 

From  tensile  coupon  tests  X and  a are  found  to  be 

X = 134  MPa  , ft  = 4.126. 

o 

What  are  the  expected  values  of  XQ  and  ft  in  3-pt  and  4-pt  flexure  tests? 
Specimen  dimensions  are: 


W 

L 

H 

1 

mm 

i 

nm 

mm 

tension 

13 

50 

4 

3- pt  flexure 

4- pt  flexure 

25 

54 

4 

4 

Use  the  volume  model  in  Table  71 

Solution 


< 

rf- 

II 

tv 

2 

600  mm 

i 

. 

L 

Vf3  " 

5400 

2 

mm 

= 5 

2 

400  mm 

I4* 

(X  ) = 

O f3 

(D* 

/ V 

0>f4 

774 

\/r 

i<7  J 

iro 

» 

V A 

ivira 

c.  Size  Effect  in  Fatigue 

For  the  renresentative  element 

Ro(t 

'V 

: I 

Probability  of  surviving  t whei 

si 

lbjected  to  a 

loading 

characterized  by  <r.. 

For  the 

entire  body 

R(t| 

fV  1 — 

l Ji r 

/ r-  ft  \ 

<T) 

| L o ’J  uv 

If 

is  uniform 

throughout  the  body. 

and 

Ro<‘ 

la-) 

R (t)  = exp 

-(t/t/ 

(759) 

j 

then 

= t 

ixn 

[ 

-V  (t/t  )^| 

/ 7 4n\ 

r L 
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’T"  f — 
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( 1)  Sample  Problem 

The  life  distribution  of  a tensile  coupon  of  unit  volume  subjected  to  a constant  stress 
a is  given  by 


RQ(tl(r)  = exp 


J 


Determine  the  life  distribution  under  a three-point  bend  stress  rupture.  Use  the 


volume  model  in  Table  71. 

Solution 

in  [R 

Jt  iO 

-amr 

o 

R(t  | s)  = exp 


From 

Table 

71, 

J 

u 

(T 

Y*  dv  * 

: /- 

s_  ' 

♦ 

s 

• ..  r 

2 I17* 

2(3  tl)  / v 

V V 

(T 

O 

/ V 

3 

O / 

^ L 

3/a  . 

/ 1 V 

* / 

\0 

. .xmu 

1 S) 

- exp 

VJ  ■ 

s i 

o / 

F 

- 

d.  Parallel  Model 
Definition 

The  system  fails  when  and  only  when 


all  subsystems  fail 

• 

■ 

Probability 

of  failure  of  the  system 

F s 

fn 

(761 

) 

O 

_ a.  _ - - t] 

oystem  reuaDiniy 

R = 

l-i 

r 

= 

i-i 

fn 

(762 

/ 

o 

(1)  Sample  Problem 


Determine  the  reliability  of  a parallel  system  of  10  bars  under  a strain  controlled 
test.  Assume  that  the  probability  of  failure  of  each  element  is 

= 1 - exp  - (e/0.6)3  , where  e is  in  %. 
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I 

i 


Solution 

Since  e is  the  same  in  all  bars,  and  since  e is  controlled,  failure  of  one  bar 
does  not  affect  the  strain  in  the  other  bars.  Therefore, 


l 

F 


N 

F 

o 


exp 


[-  (e/0.6)]3 


I 


t 

. 
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7.  DATA  AVERAGING 


a.  Compliances 

Invariants  : independent  of  fiber  orientation 


Measure:  S^,  S^.  S^,  , S^. 


Calculate:  1^,  I^.  Rj.  R2 


*1  = <s'11+S;2+2S;2)/4 

i2  = (sn+s;2-2s'l2+s;6)  /s 


Rl=  [('Sn  + S22)2+  (S'i6  + S26)2] 1/2  / 2 

I i I i o » '7X17. 

R2=  [(S11+S22-2S12-S66>  +4(S26-S16>  J /* 


Calculate  average  values  : 1^,  I2>  R^,  R2 


Average  compliances  : 


S11 

= Il+I2  ~ R1  ~ R2 

*22 

= Xl+  X2  + R1  ' R2 

®12 

= !j  - \ + r2 

*66 

= 4T2+4R2 

(1)  Sample  Problem 

Test  data  for  compliance  of  unidirectional  composite 
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1 767) 


n rr 

[ I _ i 1 4 

| |1  \\  ' | 

1 . 1 . .1  1 I I . 

! 

] { 1-  - 1 1 - _L  ! 4-1  1_  L : i i 1 

N 

1*1  f I1} 

f|  [ 1 T f ||  IT  J li  lt 

J i 1 X 1 . 1 

r * t . 

M 

: 6 x n matrix 

: 1x6  matrix 

n : no.  of  tests 

|i} 

: 1 x n matrix 

T 

M 1*1  * 1M 

■■■frl  = (M  [*]>  M 

T N 

Minimum  no.  of  tests  required 

= no.  of  components  of  F 

Distribution  of  f 

/ ; v / i \ 

r(i)2+  F F x F 

1 + F12  1 2 + F22  2 F66  6 

w 

- r a 1 + F (T  1 + F a 

■ 11  Z Z 11 

7 = 

i If. 

n . i 

(768) 
(7  69) 


(770) 

(771) 


ri  -2n1/Z 

9 ] <?72> 

F(f)  = 1 - exp  [-(f/fo)a]  (773) 


If  no  coupling  is  assumed  between  ffj  and  (o-^,  a ^),  then  one  can  define 


= F^+F^ 

C 

1 • 

(774) 

•> 

2 2 

= F_a  + F„j,  + F.  ,<r,  . 

(775) 

CL  22  2 00  0 

(1)  Sample  problem 

Failure  occurred  at  the  following  values  of  cr,  and  <r, . 

c.  b 


Determine  F^,  F > F^,  and  f. 


0" 

9 ' 

MPa 

« 

1.99 

12.78 

2 

5.  ( 

)3 

1 

11.05 

-18.62 

12.82 

, 

a 

6’ 

MPa 

31 

'.30 

-45.26 

5.  ( 

>2 

1 

51.34 

69.50 

0 

1 

1 

1 

1 

Solution 


M 


1522.78  15766.25 

591264.27 


symmetric 


-20203.75 
249450.29 
36546209.5  . 


t17 

JT  {1}  = 

f53.05 

\ 

f2  = 

= 3 

. 374 

x 1C'2  (MPa)'1 

< 

j 1522 

.78 

► 

F22  : 

= 4 

.688 

x 10'4  (MPa)'2 

(10951. 29 J 

F66  = 

= 2 

. 872 

x 1 

0‘  (MPa)' 

f : 

0.78.  1.10 

1.15, 

1. 

19, 

0 

.92, 

0.  51 
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SECTION  XV 


STRUCTURAL  ELEMENTS 
1.  LAMINATED  COMPOSITE  BEAMS  - STATIC  BEHAVIOR 

a.  Formulation 

The  classical  beam  bending  equation  is  derived  by  combining  the  appropriate  moment- 
curvature  relationship  and  equilibrium  equation. 


k = 
x 


,2 

-d  w 


dx 


M 

El 


Then 


d2M  . , v « 

— r-  + q(x)  = 0 

dx 


El  = q(x) 

dx 


(776) 


(777) 


(778) 


For  laminated  beams  it  is  convenient  to  replace  El  by  an  "effective  bending  stiffness". 
The  suitability  of  this  approach  has  been  investigated  by  deriving  a beam  theory  from 
classical  plate  theory.  The  plate  constitutive  equations  - inverted  form 


f o 

1 

~ — 

e 

A'  B' 

N j 

! 

1 

f 

s 

< 

| 1 

* ' 

1 

C'  D' 

l 

\ 1 

i 

(779) 


A solution  for  the  deflection  of  a wide  beam,  simply  supported,  loaded  by  a uniform 
moment  is  available.  The  apparent  bending  stiffness  is 


(El) 


bR2 

A fDll  1,2  * D14  R * °12> 


(780) 


where  R-depth  to  width  ratio. 

The  effect  of  bending  anisotropy  can  be  significant. 


Figure  165  Unidirectional  beam  under 


pure  bending. 

Fibers 

oriented  w to  x axis. 

a « Aft 

beam  at  supports. 

ft 

LX1 
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T 


Figure  166  (El)  /(b/D^)  as  a function  of  6:  as  a function  of  R for  a 10°  off  axis 
idir  - * 


unidirectional  beam. 


The  D'  matrix  is  defined  as  [D-BA  B]  and  involves  the  B-coupling  matrix  for  un- 
symmetrical  laminates. 

Effects  of  bending  anisotropy  on  (EI>A  depend  on  E^/ E 22>  h,  R,  and  stacking  sequence. 
Effects  of  B-coupling  on  (EI)a  depend  on  Ejj/E^.  h and  stacking  sequence. 

R large  or  small 


<EI)a>  = FT; 

(781) 

ii 

D16  = 0 


(EI)a  = 


A Dj j ~(D11  * B11  /All)b 


(782) 


(Ejj  "®n  ^A-jj)b  is  a "reduced  bending  stiffness"  and  is  comparable  to  the 

transformed  area  method  of  computing  El. 

Mid-plane  symmetry  and  no  bending  anisotropy,  B *D  =0 

ij  16 


<EI)A~Dllb 


(783) 


occasionally,  one  finds 


(EI)a  = F Exk  Ik 
k 


(784) 
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where  1^  is  computed  about  the  beam  mid- surface 


and 


cos^fl 

E11 


1 2^12  2 2 

+ (7: +7; ) sin  p cob  P + 

C12  E22 


• 4 

sin  Q 

E_  _ 

22 


(785) 


Significant  error  can  result  if  Eq.  (784)  is  used  when  B coupling  and  bending  anisotropy 
are  present. 


b/D 


11 


E1 

X 

kbl 

h3 

12 

Figure  167  Ratio  of  bending  stiffness  computed  by  Eqs.  (781)  and  (784-785) 

The  assumed  linear  strain  in  a beam  results  from  considering  bending  deformation  only. 
Shearing  deformation  can  be  significant  in  composite  beams. 

Example.  Simply  supported  beam,  point  loaded  at  center,  deflection  computed  for  the 
center  (by  virtual  work). 


8c  = 4EI  t 12  + > nr  1 = 4EI^Cb  + 6S1 

t. 


l— ti 


-shear  term 


(786) 


where 

2 

l 

= beam  length  r a I/A 

n 

= E /G  ) = form  factor 

X xz 

consider 

a unidirectional  composite,  fibers  oriented  to  x-axis 

i 

= s 

(,  4,  2 in 

E 

b 

:n 

= E = 25  x 10  psi,  G,  , = G = 0.  5 x 10^  psi 

x r 1 3 xz  r 

1 in,  t = 0.  25  in.  1 = 0.  001 3 in4 

. ..  , _ . . . . . 2 

d 

- (rectangular  Deam)  a = u.  L's  in 

2 

T = 0.0052  in 
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8b  *=  lZ/12 
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b.  Problems 


(1)  Assess  the  effect  of  B-coupling  on  (EI)^  for 
a>°°24/90°24 


b)<r12/9<r12/o°12/90°12 

O (0-2/90Vi2 

where  Ejj  = 20  x 10^  psi 


Gj2  = 0.  5 x 10  psi 


E22  = 1 x I0V  psi 


,12  = 0.3 


t = 0.  005  in. 


(2)  Use  the  equilibrium  equations  given  below  and  Eq.  (779)  to  derive  a governing 
equation  for  the  beam-column. 


dN 
x 

dx 


d2M 


= 0 


dx 


X , d2W  ^ „ 

— + N — + q(x)  = 0 

x , Z 
dx 


. 

1 — 1 1 — * — 4-  H 1—  - — A 
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2.  LONG  CYLINDRICAL  TUBES  (AXIAL  LOAD,  TORSION,  PRESSURE,  BEAM  BENDING) 


a.  Formulation 

Membrane  theory  with  anisotropy  and  B- coupling  included. 
Stress  fields  more  than  one  diameter  from  ends. 


Define 


A11 

A12- 

bi2/r 

A16  + B16/R 

A12 

A22- 

B22/R 

A26  + B26/R 

(794) 

A16 

A26  ' 

B26/R 

A66  + B66/R_ 

and 


b..  = 
11 

a. . 
ij 

(795) 

nov 

r 

o 

.iL 

i 

h 

j 

and  stres 

ts  in  a 

layer 

* 

C7. 

[ 

b., 

q -Lk 

_ . 

z 

(Q;a  1 

56i  ' 5i2  b2i]  Ni  (797) 

N 

Lxial  load 

L . 

X 

Ne 

= I 

>R 

- pressure  load 

• 

1ST 

xe  - 

torsion  load 

Axial 

compliance 

b. . 

i 

€ 

X 

:> 

1 1 

h 

N 

t 

S = b 

ll/h 

(798) 

Effective 

axial 

modulus 

E 

1/S 

= h/b. 

(7991 

X 

X 

J 

L 1 

Cylindrical 

I 

tube  in  beam  bending 

(approx. 

) 

,4 

E 

I 

X 

dx4 

= q 

(x) 

(800) 

Tubes  with  low 

\ 

, terms 

(all  0* 

. o 

1/90), 

include  effect  of  shear  deformation.  Form 

factor 

, * 

= 2 



4 
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b.  Problem 


Design  a laminated  tube,  0°./±45°.  ply  configuration,  and  calculate  the  stress  in  the 
innermost  and  outermost  45°  and  0°  layers. 

R . = 2 in. 

1 6 I 

En=  20x10  psi  ej  = +0.066,  -0.005  Nx  = 2300  lbs/iu. 

E = 1 x 106  psi  e * = +0.003,  -0.006  N „ = 500  lbs/in. 

£ - . c x A 


°12 

= 0.  5 

X ] 

[0  psi  e,  - 

o 

o 

-H 

. 01 

5 

V12 

= U.  3 

t = C 

1.  005  in. 
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306 


< 


3.  COLUMNS  (COMPRESSION  LOADED  1-D  MEMBERS) 


a.  Formulation 

Laminated  column  = 0 or  D ^ ^ / D ^ j very  small 

Bending  stiffness 

IEI,A  = 'Dii-Bu2/An»b 


Column  equation 


(EI,  ^ .pA  * " 
A dx  dx 


(801) 


(802) 


Euler  buckling  load,  simply  supported  columns 

nZ  (Ei)A 


Pe  = 7T 


(803) 


For  general  end  restraints  = M.  /a.  @ x = o 

dx  11 


^ w M /ft  @ x = L 

dx  m2,a2u 


(804) 


(El), 


(El), 


defining  X , a — 

1 a 

.rv  , 

jL  ’ *2-  <*2L 

i = kL 

, k2  = p/(EI)  . 

the  characteristic  eauation  is 

i 1 *2.  A 

2 

U-A 

1 ‘ A 2 ‘ ' 

V 1 X2  * ) * 

sin  # + (2  + X ^ # ) cos  #-2  = 0 

For  special  boundary  conditions,  the  effective  length  approach  can  be  used  with 


Leff  = KL 


ff2(ED 

3 2. 

cr  2 
Leff 


Buckled  shape  of 
column  is  shown  by 
dashed  line 


Theoretical  K value 


Recommended  design 
value  when  ideal  con- 
ditions are  approxi- 
mated 


4 — )- 


End  condition  code 


0.  5 


0.65 


T 


0.7 


0.80 


1.0 


1.2 


Y 


1.0 


1.  0 


2.  0 


2.  10 


*Y  Rotation  fixed  and  tranaintion  fixed 

Y Rotation  free  and  translation  fixed 

hf3  Rotation  fixed  and  translation  free 

Rotation  free  and  translation  free 


(805) 


(806) 


2.  0 


2.  0 
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Role  of  Bjj  coupling 

(1)  Reduced  bending  stiffness  - lower  P 

cr 

(2)  The  reduced  bending  stiffness  can  also  be  calculated  by  the  transformed  area 

method.  This  illustrates  the  additional  effect  of  shifting  the  columns  effective 
centroid  (neutral  axis  for  bending) 

/ b / 

On 

I b 

Figure  169  A 'transiormed  section1,  all  0 , b^E^/E  ^ )b. 

Therefore,  axial  loads  applied  coincident  to  the  mid  surface  produce  bending 

deformations  with  possible  failure  due  to  excessive  deformation  with 

P < P (reduced) 
cr 

Role  of  low  interlaminar  shear  modulus  G 

xz 

Reduces  resistance  (total  effective  stiffness)  to  out-of-plane 
deformation  (buckling) 

P = k P 
cr  cr 

where 

k = , ■ - form  factor 

' cr 
AG 

Laminated  circular  tube  - compression 
no  anisotropy 

(E!)  A = EXJ  (809) 

E defined  by 
x 

Secondary  column  failure  modes 

For  cylindrical  tubes,  high  R/h,  high  compressive  load,  there  can  be  cylindrical 
buckling  - many  analytical  tools  available. 

Anisotropy  and  heterogeneous  - Cheng  and  Ho 
heterogeneous  - Tsai 

orthotropic  - NASA  and  Air  Force  Design  Guides,  Handbooks 


(807) 

(808) 
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Homogeneous,  orthotropic  approximation 


O = KE,  (h/D) 
cr  1 

K is  evaluated  for  various  laminates  and  materials 


Boron  Epoxy 
U.  H.M.  Graphite /Epoxy 
H.  M.  Graphite/Epoxy 
H.S.  Graphite /Epoxy 
Fiberglas 


N. 

Ns 


10  20  30  40  30  60  70  00  00 

LAMINATE  TYPE  1*01,  On 


Boron/ Epoxy 

— •• — U.H.M.  Graphite /Epoxy 

H.M.  Graphite /Epoxy 

H.S.  Graphite /Epoicy 


Figure  170  Buckling  parameter: 

(a)  [±  9]  laminate; 

(b)  [0/90]  laminate. 


80  70  60  SO  40  30  20  10 

PC8CENT  0*  F0«  0-90"  LAMINATE  TYPE 


b.  Problems 

(1)  Develop  an  escpression  for  the  layer  stresses  in  a laminated  beam-column. 

(2)  Siae  a laminated  tube  for  the  given  compressive  load  and  geometry 


p = 

s Z0,  000  lbs 

K ! 

c 2 in. 

i = I D It. 

K = 

« 1 

l 

E1 

= 40  x 10  psi 

C12  = °- 

5x10  psi 

t = 0.  005  in. 

E * 1 x 10  psi 

V,- 

= 0.! 

1 

B*  1 

( 3)  Assess 

b a ■ 

the  effect  of  low  shear  s 

itiffness  on  P 

. 

(4)  Re-evaluate  the  design  of 

(2)  for  local 

cr 

. instability. 
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4.  COMPOSITES  FOR  SELECTIVE  REINFORCEMENT  OF  STRUCTURAL  ELEMENTS 


a.  Formulation 

Standard  structural  forms  can  be  reinforced  by  unidirectional  (0°)  composites  for 
enhanced  specific  stiffness  and  strength. 

Examples 


YUTTT1 

TTTT7 

•>  ■■ — - =; 

=55F  ? 

Ez 

lOJ  l u 

(LI  U 

Beams 


Panel  Stringers 


Mid  plane  symmetric  no  symmetry 

(EI)A  = X(EI)  (811) 


(El)  A = D - B /A 
A x x x 


or  transformed  area  method 


Specific  stiffness 

/ - £(EI)  - 

( I)A//1T  E(PV)  " ^PV1 


bnear  stress  in  adnesive 

- IQ. 


(812) 


(813) 


' ib 

(814) 


I - moment  of  inertia,  transformed  area,  all  material  (I) 
b - width 
V - shear  force  at  x 

Q - first  moment  of  area,  composite  reinforcement 
about  centroid  axis 

Secondary  failure  mode 

local  buckling  - compare  to  sandwich  beam  face  sheet. 


V 
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b.  Problem 

rtr 

LI  1 II  1 II  1 

T — 1 — 1 — 1 — I — f — ! — T" 

4 — J 1 1 - | 1 i 1 

specific  stiffness. 

j i I * 

An  aluminum  stri 

L i j j T } — i - i • r 

neer  reiniorced  panel  is  to  be  upgraded  for  higher 

specific  i 

l _ I i - l j 
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1 T ’ - 
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r~ 
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t h 

P 
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—¥ 
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, 4 1 

*-  t 

. -L  h 
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s 
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f 1 — ' 

i» 

zzz 

n± 

= 1.  0 in 

6 

b = 4.  0 in. 

1 s c 

x = 

A 

", . H 

• E 

A /b 

b 

- stringer  spacing 

X 

-• 

B S 

B = 
X 

B11  ' 

► z 

1 

E A /b 

B s s 

subscript  s - relates  to  stringer 

D = 

E 

• E 

I /b  + Z 

2E  A /b 

X 

11 

1 

B S 

s s s 

E 

11 

= 35  x 10 

psi  Gj2  = 1 X 10  I 

>si 

E 

6 

= 3.  5 x 10  v..  = 0.25 

id 

p 

AL 

= 0.1  lbs/in  p = 0.075  lbs/in 
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, 1 1 
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5. 

THERMAL  STRESSE 

:s 

. 

i r~ 

T~ 

T 

* 

b 

f- 

it 

, , zr- 

* 

- 

Jn 

j 

(7 

Qij 

1 

Nk 

- e 

iJn 

RhlV 

Ji6 

3, . 

* 

^i2b2j 

h/2 

r 

where  N. 

= 

Qi 

. e 

. dz 

L 

l/2 

e. 

_ 

O'. 

AT 

1 

or. 

— 

lamina  coef 

of  thermal 

expansion  in  x . fi  cnnrrf. 

i 

L 

(a 

X*’  a n 

from  a, . aJ 

U X 

1 C. 
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. 
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- _* L 

i — L 

l 
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(815) 

(816) 

(817) 
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6.  SIZING  FOR  STIFFNESS 

a.  Formulation 

Stiffness  critical  designs  : bending  - deflection 


compression  - stability- 
dynamic  - frequency  control 

In  1-D  structures  (beams,  columns),  unidirectional  laminates  will  predominate  in  designs 

In  2-D  structures  (plates,  shells),  angle-plied  laminates  become  more  efficient. 

Plate  stability  (assume  quasi-homogeneous,  orthotropic  : B. . =0,  D ,,  D , - negligible) 

ij  16  26  =>  o / 

simply  supported,  compression  loaded  on  X edges 


n ‘ 

2 . 

| 

C.  IT  ] 

N = r- 

xcr  ,2 

b ' 

!<D11D22>  +D122D66| 

; (8i8) 

or 

’xcr  = K(h/b)  (819) 

Numerous  sources  exist  which  give  K for  different  materials  and  laminate  configurations. 
Examples  are  given. 




1 1 1 1 — 1 1 r ■ 

1-  i J - 1 1 

I 1 1 

l_J 

, 

1 TT ? 

L I.  1 1 11  i 1 II  _j_!  j__L 
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0°/90°/±45°  laminates  PLATE  STABILITY  OF  0°,  90°,  AND  ±45°  FAMILY  OF  PLATES 

example 

Loading 


Thickness  Fraction  of  Ply 
Orientation,  % 

K x 10'6  psi 

cr 

Symmetric 

Unbalanced 

. _o  o 

45  0 

90° 

135° 

Plvine  Order 

Plvine  Order 

Plvine  Order 

50  0 

0 

50 

54.5 

45.8 

27.  1 

40  20 

0 

40 

50.0 

42.7 

25.2 

30  40 

0 

30 

45.  5 

40.7 

25.8 

20  60 

0 

20 

41.1 

39.  1 

29.0 

10  80 

0 

10 

36.  5 

36.8 

33.5 

0 100 

0 

0 

32.0 

32.0 

. 0 

40  0 

20 

40 

50.0 

42.8 

c-j  .Z 

30  20 

20 

30 

45.  5 

40.8 

27.7 

20  40 

20 

20 

41.1 

39. 1 

29.6 

10  60 

20 

10 

36.  5 

36.8 

30.4 

0 80 

20 

0 

32.0 

32.0 

24.  3 

30  0 

40 

30 

45.  5 

40.5 

25.8 

20  20 

40 

20 

41.1 

39.2 

29.  6 

10  40 

40 

10 

36.  5 

36.8 

28.2 

0 60 

40 

0 

32.0 

32.0 

18.5 

20  0 

60 

20 

41.1 

38.9 

29.0 

10  20 

60’ 

10 

36.  5 

36.8 

30.4 

0 40 

60 

0 

32.0 

32.0 

18.  5 

10  0 

80 

10 

31.2 

32.  6 

33.2 

0 20 

80 

0 

32.0 

32.0 

24.  3 

Q 0 

RHU 

0 

20.8 

20.8 

20.8 

Quasi -Isotropic 

25  25 

25 

25 

43.  3 

40.0 

28.  6 

1.  Quasi- homogeneous  plying  orders  contain  a sufficient  number  of  consistently  repeated 

sequfnces  of  lamina  to  make  negligible  the  effects  of  resulting  bending  couples.  The  basis 
for  these  sequences  is  45°,  0°,  90°,  and  135°, 


2.  Symmetric  plying  orders  contain  balanced  ply  orientations  where  the  sequence  of  lamina 
directions  are  exactly  mirror  image  about  the  laminate  midplane.  The  bases  for  these 
are  45°.  0°,  90°,  135°,  90°,  0°,  and  45°. 

3.  Unbalanced  plying  orders  contain  only  one  sequence  of  orientations  based  on  45°,  0°,  90°, 


and  135' 

• - — | 

1 

t 
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cr 

= K (z) 

t 

Thickness 


It  is  evident  that  for  Nx  loading  the  ±45°  plate  gives  the  maximum  resistance  to  buckling. 
This  results  from  the  dominant  2D^  term  in  Eq.  (818).  Since  strength  must  also 

be  adequate,  0°  plies  may  be  added  or  ± 6 laminates  at  lower  angles  may  be  considered. 

Hole  of  stacking  sequence  - example 

± 45°  plies  on  the  outer  surfaces  produces 
the  maximum  D,  ,. 


TABLE  73  EFFECTS  OF  ORIENTATION  SEQUENCES  ON  STABILITY  OF 
SYMMETRICAL  (BALANCED)  QUASI- ISOTROPIC  PLATES 


Orientation  Sequenc< 

of  Laye 

K x 

10 

-s 

>si 

cr 

4 

5° 

l 

35° 

0° 

90° 

90 

o 

o 

° 

135° 

45 

O 

50 

.6 

13 

5° 

0° 

45° 

90° 

90 

o 

45° 

0° 

135 

O 

44 

.0 

0° 

on° 

90 

o 

45° 

90° 

1 

o 

O 

4n 

A 

rn 

/vQ 

0° 

o 

no 

Q 

5 

U 

9 

1 33“ 

1 33 

0 

**3 

. u 

0° 

135 

9 

0° 

4 

5~ 

45 

90° 

1 35u 

0° 

37 

.7 

0° 

90° 

— 

4 

5° 

135° 

135 

o 

45 

0 

90° 

0 

O 

34 

.7 

1. 

All  plat< 

iS 

ire 

equal 

in 

thicknes 

1. 

2. 

Loadine 

and  laminate  olane  coordinate 

58 

ire 

the  sam 

e as  ir 

Tabl 

e 7 

2 

Consider 

compression  strength 

= 

P 

Jj 

h 

(820) 

Requirin 

8 

r 

= 

and  defining  as  a 

ontimum 

i 

3 

cr 

3 

• 

2 

<T 

_ 

K(h/b 

)2 

N 

/h) 

2 

(821) 

1 

OrT 

xc 

r > 

X 

eliminates  h ai 

< N 

2/ 

3 

I 

anDT 

a 

<I/3 1 



■ 

(822) 

/ N 

x\ 

2/ 

3 

J 

' 

i 

- 

- 

structural 

index 

— . 1 

A 

D / 

ZD 

— 

J 



_ 

_ 

±j 

zd 

1 

J 

1 

Combining  with  a laminate  failure  criterion,  the  optimum  laminate  can  be  designed  for  the 


specific  class  being  considered  (±6,  ±9  /0°,  0°/90°/±45°) 

m n 


Figure  172  Examples  of  Uniaxial  strength  allowables  (compression). 


b.  Problems 

(1)  Calculate  the  individual  bending  stiffnesses  (D's)  of  Eq.  (818)  for  the  following 
laminates  and  material 

[+45,  -45,  +45,  -45] 

[ 0.  90,  0,  90] 


[+4! 

5,  -45,  0,  90] 

J s 

L. 

i 

=11 

- ' 

SO  x 10  psi  C 

"12 

= ( 

1.5 

X ] 

10 

psi 

t 

= 0.005 

—4 f- — l — -i- — 

E 

^22 

6 

l . 5 x 10  psi  t 

; 

12 

= t 

1.3 

(2)  Using  the  charts  and  tables  available  for  strength  and  K for  (0°/±45°/90O)  laminate 
class,  determine  a minimum  weight  design  for  a compression  panel  where 

N = 5,000  lbs /in.  N = N =0 

x y xy  . 

a = b = 10  inches 

Plot  WT/b  vs  N / b for  stability  controlled  and  strength  controlled  designs 
i X 2/5 

plot  "opt  v,(Vb) 


-+ — I — — f- 


tiJ  I LI 


JJ.  L 
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DESIGN  PROBLEM  - STRINGER  REINFORCED  COMPRESSION  PANEL 


SECTION  XVI 


SPECIMEN  CONFIGURATIONS  AND  LOADINGS 
1.  TESTS  FOR  LAMINA  PROPERTIES 

a.  Properties  to  Be  Measured 

1.  TENSION  - El,  Et,  Vlt,  Fl\  Ff 
COUPON,  SANDWICH  BEAM,  TUBE 

2.  COMPRESSION  - EL , ET,  VLT,  FLC,  f' 

COUPON,  SANDV.'ICH  BEAM,  TUBE,  SHORT  BAR 

3.  INPLANE  SHEAR  - GLr,  F*T 

[t45]  COUPON,  CROSS  SANDWICh  BEAM,  TUBE 

4.  FLEXURE-  Ef,  Ff 
BEAM 

5.  INTERLAMINAR  SHEAR  - G,s,  F,s 
SHORT  BEAM 

6.  GENERAL  STATE  OF  STRESS 

OFF-AXIS  COUPON,  TUBE,  CROSS  SANDWICH  BEAM 


b.  Tension  and  Compression 


*SiQE  SUPPORT 


323 


gitudinal 


and  Interlaminar 


(±45  SPECIMEN) 


L * 06,  W * 0 25.  T«  0075-0  090 
S»0.4 

T * 3P/4WT 

FIBERS  PARALLEL  TO  LENGTH 


Specimen  configurations:  (a) 
(b)  interlaminar  shear. 


longitudinal 


CENTER  RAIL 
SLIDES  THROJGH 
GUIDE 


STRAIN 

GAGES 


STRAIN  PACE 
LOCATION 


e 175  Rail  shear 


double, 


tests 


single 


T i | r | | I i 

-1-4-  4-4  — f 

Flexure  Tests  I _J 

X-OT  Cl  CVI  IOC 


3-PT  FLEXURE 
(Longit  Jmo!  Flexure) 
,P 


hwi_i 


4-PT.  FLEXURE 
(Trcnsverse  Flexure) 

P/2  P/2 

£s/2-*| 


Ls4,  W*  0 5,  1 = 0,075-0.090,  S=25 
<7  - 3PS/2WT2 
E,  = (S5/4\VT3)(fVU) 

U DEFLECTION 


04 


L=  3,0,  W=  0.5,  T=  0075-0090,  S=2.0 
a « 3PS/4WT2 
E(=  (IIS5/64WT5)(P/U) 

It  DEFLECTION 


(a)  (b) 

Figure  176  Flexure  test  configurations:  (a).  3-pt.  ; (b)  4-pt. 
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Figure  1 

10  15  20 

L/h 


shear  stress  in  flexure  tests, 
e.  Sandwich  Besm  Test  i I 


SANDWICH  BEAM  DETAILS 


H — t — ( — H 


ORIENTATION 

No  OF  PLIES 

SECONO  FACE 
MATERIAL 

CORE  ALUMINUM 
HONEYCOMB  DENSITY 
(pcf) 

0* 

6 

STAINLESS  STEEL 
( 125"  THICK) 

23 

90* 

8 

GLASS/EPOXY 
( 08"  THICK) 

4.5 

CT  • 4P/t[l.5  » (TtTf  ) /z] 


1 

. 

u 

. Figure.  178 
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nfieuration. 
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, 

LLJ 
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t rxttrJ 

1 1 J 

1—1 1 J 1 1 J L-  1 „ i ..L  . 




HYDROSTATIC  SEAL  ( 06) 


(DiNSlOE  SUPPORT 
(2)  INSiOE  Tab 


Tube  configuration, 


Comparison  Among  Different  Test  Methods 


LONGITUDINAL  TENSION 


TRANSVERSE  TENSION 
B/Ep 


COUPON  30.1  21 

FLEXURE  28  7 - 

SAND  BEAM  32  2 .50 


COUPON  2.73 
FLEXURE  2.76 
SAND  BEAM  3 56 


04 

STRAIN,  % 


STRAIN,  % 


longitudinal;  (b)  transverse.  [43] 


e 181  Shear  stress-strain  r< 


elation*.  [44J 


e Test 

TYPE  OF 
SPECIMEN 


TEST 

TEMP.  (°F) 


ORIENTATION 


(MSI) 


COUPON 


S.  BEAM 


COUPON 


COUPON 


S.  BEAM 


T300/5208 


Ic  • 0.5"  FOR  COUPON 


.34 

218 

.39 

247 

.31 

206 

.50 

214 

TABLE 

74  COMPRESSION  PROPERTIES  [45] 
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problem^  Resulting  from  ani  sSotRopy  and  nonhomogeneity 


Coupling  Effect* 


Shear  coupling  * S1A’  S?a  ^ 0 


Extension  - twisting  coupling 


SMtAR  COURUNO 


CONVERSELY 


Figure  182  Effect  of  shea 


r coupling 


Figure  183  Stress  concentration  at  free  edge  of  [45/-45]  angle 
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Figure  184 
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APPENDIX  A 


SYSTEME  INTERNATIONALE 


1.  METRIC  PREFIXES 


Prefix 


Abbreviation 


Multiplier 


tera- 

T 

10 

giga- 

G 

10 

mega  - 

M 

10 

kilo- 

k 

10 

hecto- 

h 

10 

deca  - 

da 

10 

deci- 

d 

10 

centi- 

c 

10 

milli- 

m 

10 

micro- 

H- 

10 

nano- 

n 

10 

pico- 

P 

10 

feroto- 

f 

10' 

atto- 

a 

10' 

2.  CONVERSION  EQUATIONS 

l kg  = 2. 20  lb 

1 kgm  3 = 3.61x10  3 lb /in  3 

IN  = .225  lb-f 

IN  = .102  kg-f 
1 Pa  = 1 . 45  x 10  4 psi 


1 lb  = . 454  kg 
I lb/in3  = 2.  77  x 104  kgm'3 
1 lb-f  = 4.45  N 
1 kg-f  = 9.  81  N 
1 psi  = 6.  895  Pa 


1 

MPa 

= 

. 145  ksi 

1 

ksi  = 6 

»• 

1 

GPa 

= 

. 145  1 

0 ps 

i 

1 

„6 

0 psi  = 

1 

Pa  = 

. 

102  kg 

-f  m" 

2 

1 

1 r -2 

kgf  m 

1 

MPa 

= 

. 10Z  ke-f  mm 

.2 

1 

kg-f  mix 

x 

Nm"  J 

l = 

, 00571 

iKf /i, 

n 

1 

1 hf  / in 

; 

1 

Mm'1 

l_ 

1 Al  L 

I 

l 1 l/l  f All 

a " a i 

1 

NTT) 

- ( 

, IU^  Kgl  / TTi 

1 

kgf/m  : 

s 

L— J 

I 

Nm 

= 8.98  ii 

i-lbf 

1 

inlbf  i 

= 

i 

L—JJ 

-1 


APPENDIX  B 


TRIGONOMETRIC  FUNCTIONS 

. 


C0S  20 

CBS  40 

SIN  20 

SIN  40 

THETA 

(DEG) 

-90 

-1.00000 

1 .00000 

-.00000 

.00000 

-90 

-85 

-.98481 

.93969 

-.17365 

.34202 

-85 

-80 

-.93969 

.76604 

-.34202 

.64279 

-80 

-75 

-.86603 

.50000 

-.50000 

.86603 

-75 

-70 

-.76604 

.17365 

-.64279 

.98481 

-70 

-6  5 

-.64279 

-.17365 

-.76604 

.98481 

-65 

-60 

-.50000 

-.50000 

-.86603 

.86603 

-60 

-55 

-.34202 

-.76604 

-.93969 

.64279 

-55 

-50 

-.17365 

-.93969 

-.98481 

.34202 

-50 

-45 

.00000 

-1.00000 

-1.00000 

-.00000 

-45 

-40 

.17365 

-.93969 

-.98481 

-.34202 

-40 

-55 

.34202 

-.76604 

-.93969 

-.64279 

-35 

-30 

.50000 
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-.86603 

-.86603 

-30 

-25 

.64279 

-.17365 

-.76604 

-.98481 

-25 

-20 

.76604 

.17365 

-.64279 

-.93481 

-20 

-15 

.86603 

.50000 

-.50000 

-.86603 

-15 

-10 

.93969 

.76604 

-.34202 

-.64279 

-10 

-5 

.98431 

.93969 

-.17365 

-.34202 

-5 

0 

1,00000 

1 .00000 

.00000 

.00000 

0 

5 

.98481 

.93969 

.17365 

.34202 

5 

10 

.93969 

.76604 

I 

.34202 

.64279 

10 

15 

.36603 

.50000 

.50000 

.36603 

15 

20 

. 76604 

. 1 73  65 

.64279 

.98481 

20 

25 

.64279 
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25 

30 

.50000 
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.86603 
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30 

35 
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.64279 

35 
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.98481 
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40 
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.00000 
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45 

50 

-.17365 
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50 

55 
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-.64279 

55 

CO 

-.50000 
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60 

65 
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65 

70 
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.17365 

.64279 
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70 

75 
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75 

80 
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80 

85 
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.93969 
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85 

90 
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APPENDIX  C 


NORMAL  DISTRIBUTION 
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1 02 

8461 

1539 

1 52 

9357 

06*3 

03 

3120 

4880 

53 

7019 

2981 

1 03 

8485 

1515 

1 53 

9370 

0630 

04 

5160 

4840 

54 

7054 

2946 

1 04 

8508 

1492 

1 54 

9382 
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05 

5199 

4801 

» 

7 086 

2912 

1 05 

8531 

1469 

1 55 

9394 

.W,06 

06 

5239 

4761 

56 

7123 

2877 

1 06. 

8554 

1446 

I 56 

9406 

0594 

0? 

5279 

47. 

57 

7 157 

2843 

1 07 

8577 

1423 

1 5* 

94  18 

••-.82 

08 

5319 

4681 

58 

7190 

.'**10 

1 0« 

8599 

i 40 1 

1 S-s 

9429 

0571 

09 

5359 

4641 

59 

7224 

27  7*. 

1 0 , 
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1379 

1 59 

••441 

0559 

10 

5398 

4602 

60 

7257 

27  4 1 

1 10 

8643 

1357 

1 60 

9152 

O'.  48 

II 

5438 

4562 

61 

7291 

2709 

I 11 

8665 

1335 

1 61 

9463 

0537 

12 

54  78 

4522 

62 

7324 

26  7 1. 

1 12 

8686 

1314 

1 62 

9474 
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13 

5517 

4183 

63 

7357 

2443 

1 13 

8708 

1292 

1 03 

9484 

0516 

14 

5557 

4443 

** 

7389 

2611 

1 14 

8729 

1271 

l 64 

9495 

0505 

13 

5596 

4404 

65 

7122 

2578 

1 IS 

8749 

1251 

1 65 

9505 

OHS 

1C 

5636 

4364 

66 

7451 

2 *-4  6 

1 1<> 

8770 

1230 

1 6o 

9515 

.485 

IT 

5675 

4325 

67 

7iv 

.11 

1 1* 

8790 

1210 

1 6' 

9525 

"•75 

18 

5714 

4 286 

68 

7 517 

. 483 

1 16 

6810 

1190 

1 68 

9.535 

0<65 

19 

5753 

4247 

69 

7549 

21.51 

1 19 

8830 

1170 

1 69 

954  5 

0 4 55 

20 

5793 

4207 

70 

7 5s- 

- 120 
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8849 

1151 

1 7u 

9554 
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21 
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4168 

71 

7611 

2 l“9 
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8869 

1131 

1 71 
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OJ  <6 

22 
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72 

7642 

. < .s 
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23 

5910 
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73 

2-.27 

1 23 

8907 

1093 

1 73 

a 

. .18 

24 

5948 

4052 

74 

7704 
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1 24 

8925 

1075 

1 74 

9591 

0409 

25 

5987 

4013 

75 

7734 

2266 

1 25 

8944 

1056 

I 75 

9599 

0401 

26 

6026 

3974 

76 

7761 

2236 

1 26 

8962 

1038 

1 76 

960S 

0J92 

27 

6004 

3936 

77 

7791 

2206 

1 27 

8980 

1020 

1 77 

9616 

0384 

28 

6103 

3897 

78 

7823 

2177 

1 28 

8997 

1003 

1 78 

0.175 

29 

6141 

3859 

79 

7852 

2148 

1 29 

9015 

0985 

1 79 

9633 

OJ67 

30 

6179 

3821 

80 

7881 

2119 

1 30 

9032 

0968 

1 80 

M41 

0359 

31 

6217 

3783 

81 

7910 

2090 

1 31 

9049 

0951 

1 81 

9649 

0351 

32 

6255 

3745 

82 

7939 

2061 

1 32 

9066 

0934 

1 82 

9656 

UJ  4 4 

33 

6293 

3707 

83 

7967 

2033 

1 33 

9082 

0918 

1 83 

*664 

0316 

34 

6331 

3669 

8' 

7995 

2005 

1 34 

9099 

0901 

1 84 

9671 

0329 

6368 

3 132 

85 

802J 

1977 

1 35 

9115 

0885 

1 85 

9678 

0322 

36 

6406 

3594 

86 

805 1 

1949 

1 36 

9131 

0869 

1 86 

9686 

0314 

37 

6443 

3557 

87 

8079 

1921 

1 37 

9147 

0853 

1 87 

9693 

0307 

38 

6180 

3520 

R8 

8106 

1894 

1 38 

9162 

0838 

1 88 

9699 

6301 

39 

6517 

3483 

89 

8133 

1867 

1 39 

9177 

0823 

1 89 

9706 

0294 

40 

6554 

3446 

90 

8159 

1841 

1 40 

9192 

0808 

1 90 

9713 

0287 

41 

6591 

3409 

91 

8186 

1814 

1 41 

9207 

0793 

1 91 

9719 

0281 

42 

6628 

3372 

92 

6212 

1788 

l 42 

9222 

0778 

1 92 

9726 

0274 

43 

6664 

3336 

93 

8238 

1762 

1 43 

9236 

0764 

1 93 

9732 

0268 

44 

6700 

3300 

94 

8264 

1736 

1 44 

9251 

0749 

1 94 

9738 

0262 

45 

6736 

3264 

95 

8289 

1711 

1 45 

9265 

0735 

1 95 

9744 

0256 

4fi 

6772 

.1228 

96 

8.115 

; 85 

I 46 

92  79 

0721 

1 96 

9750 

025 0 

47 

6808 

3192 

97 

8340 

1660 

1 47 

9292 

0708 

1 97 

9756 

0244 

48 

6844 

3156 

98 

8365 

1635 

I 48 

9306 

0694 

1 98 

9761 

0239 

49 

6879 

3121 

99 

8389 

1611 

1 49 

9119 

0681 

1 99 

9767 

0233 

50 

6915 

3085 

1 00 

8413 

1587 

1 50 

9332 

0668 

2 00 

9773 

0227 

APPENDIX  D 

CHI -SQUARE  DISTRIBUTION 
PERCENTAGE  POINTS 


F(x'> 


r 

Jo 


(IY 


w-2 
2 „ 


APPENDIX  E 


GAMMA  FUNCTION 


Value»  of  I'  (n)  - [ t >dx.  (n  + I)  - nHn) 

Jo 


" 

1 (n)  ^ 

n 

I 1 H) 

I 

rw 

n 

r in) 

1 00 

1 00000 

1 2 ft 

90640 

I .ftO 

88623 

1 7ft 

91906 

1 01 

133 

1 26 

90440 

1 51 

88659 

1 76 

92137 

1 02 

98884 

1 27 

90250 

1 ft  2 

88701 

1 77 

92376 

1 03 

983.7ft 

1 28 

90072 

1 53 

H87.ft7 

1 78 

92623 

1 04 

9784  t 

1 29 

8990 1 

1 54 

88818 

1 79 

92877 

1 Oft 

97  JftO 

1 30 

89747 

1 55 

88887 

1 80 

93138 

1 Oft 

90874 

1 31 

891, 00 

1 56 

88964 

1 81 

93408 

1 07 

9041ft 

1 32 

89464 

1 ft7 

8904*1 

1 82 

1 08 

9.ft973 

1 33 

89  <38 

1 58 

89  142 

1 83 

9.3969 

1 09 

95.716 

1 34 

89222 

1 59 

89243 

1 84 

94261 

1 10 

95135 

1 3ft 

89115 

1 60 

89352 

1 8 ft 

94561 

1 II 

9 4739 

1 36 

89018 

1 61 

89468 

1 86 

94869 

1 12 

943ft* 

1 37 

88*131 

1 02 

89592 

1 87 

95184 

1 It 

93993 

1 38 

R88ft4 

1 63 

89724 

I 88 

95507 

1 1 i 

93642 

1 39 

8878ft 

1 64 

89  HO  4 

1 89 

9 .838 

1 1ft 

93304 

1 40 

88726 

1 6 ft 

90012 

} 90 

96177 

1 >6 

92980 

1 41 

881.76 

1 66 

90167 

1 91 

96523 

1 IT 

92670 

1 42 

886  16 

1 67 

90330 

1 ‘*2 

96878 

1 18 

92373 

1 43 

8860  4 

1 68 

90500 

1 93 

97240 

1 19 

92088 

1 44 

88  ft  80 

1 69 

'•0678 

l.*4 

97610 

1 20 

91817 

1 4 ft 

88ft«ft 

1 70 

90864 

1 95 

97988 

l 21 

91  ft.ftH 

1 46 

88  ft  60 

1 71 

91057 

1 96 

98374 

1 22 

91  til 

1 17 

88563 

1 72 

91258 

1 97 

98768 

1 23 

91075 

1 48 

88  ft  7. ft 

1 71 

9 1 466 

1 98 

99171 

1 24 

908  ft  2 

1 4*1 

88ft*ift 

1 74 

91683 

1 99 

99581 

2 00 

1 00(8)0 

• Kyr  larxe  positive  value*  of  i,  Hx)  appruionatea  the  asymptotic  aeriea  x*r  » \^*  T I -f 
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